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Introduction 



From an etymological point of view, physics means Science of Nature. However, 
over the course of time the word 'science' had his definition changed. In early 
times, such as Ancient Greece, 'science' meant the knowledge one has about 
one subject. Later, in the early XIX]3century, 'science' meant the collection of 
'mathematical science', physics and 'natural science', where 'science' in its old 
meaning is used to define the modern concept. 

On the other hand, the definition of physics itself has changed all over the 
centuries. For instance, physics as the science of nature was one of the three 
parts of philosophy, according to Plato (about 400 B.C.), together with ethics 
and logic. Note that Aristotle (about 350 B.C.) also included physics, or the 
science of nature in his definition of philosophy, which distinguished theoretical, 
practical and poetic philosophy. 

It is indeed interesting to note that in the Xlljjcentury, physics had a double 
meaning; namely 'medicine' and 'what is related to nature'. Note that the word 
which designates someone who practices medicine is 'physician'. It is obvious 
that the words physicist and physician have the same etymologic origin. In 
the XVJjcentury, 'physics' was defined as the knowledge of natural phenomena 
and was also called natural philosophy. The title of the celebrated book by 
Newton (around 1700) is indeed PhilosophiceNaturalis Principia Mathematical 
Latin for 'Mathematical Principles of Natural Philosophy\ 

It is only from the period of Galileo and Descartes (around 1600) on that 
physics took its modern meaning of classical Newtonian physics. As an illustra- 
tion, the first edition of the French Academy Dictionary (1694) defines physics 
as Science which goal is the knowledge of natural things. In its eighth edition 
(1935), the Dictionary defines 'physics' as Science which observes and classify 
material world's phenomena in order to extract the underlying laws; which is 
closer to the modern definition of physics: The scientific study of forces such 
as heat, light, sound, etc., of relationships between them, and how they affect 
objects., according to the Oxford Dictionary. 

Anyhow, modern physics is composed by various domains, such as con- 
densed matter physics, astrophysics, particle physics, among many others. 
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Moreover, there has always been an intrinsic link between theory, which for- 
malizes abstract mathematical concepts and leads to observable predictions, 
and experiments, which try to measure phenomenon to which theory gives a 
conceptual representation. The experiments are guided by theories, while theo- 
ries are developed when experiments reveal new phenomena, this is a constant 
interplay. As an example, in the '60s, Gell-Mann found a symmetric struc- 
ture in hadronic particles, which led him to postulate the existence of quarks. 
These have been observed from late '60s to mid '90s. This is when theory 
influences experiments, but the converse is true as well: experiments on black 
body radiation led Planck to postulate the quantization of the energy of light. 

On another hand, since the beginning of the XXJjcentury, there has been 
a distinction between theoretical physics and experimental physics, unlike in 
chemistry where many brilliant theorists are also brilliant experimentalists. 
This can be a consequence of the deep mathematisation of physics, which relies 
on advanced concepts of mathematics such as group representation theory, dif- 
ferential geometry, algebraic geometry, topology, etc. Such a mathematisation 
has opened the area of theoretical and mathematical physics, at the boundary 
of mathematics and physics. This is an inter-dependent link between the two 
disciplines. 

As an example, Newton and Leibniz invented differential calculus in order 
to formalize the laws describing motion of bodies, and Einstein relied on works 
by Riemann shedding the foundation of differential geometry around 1850 in 
order to develop his theory of General Relativity (1915). 

General Relativity is precisely the framework in which this thesis has been 
realized. But what is General Relativity? There are entire books dealing with 
this question, but we will try to give a very brief idea of the main concept of 
General Relativity. 

The theory of General Relativity states that there is a intrinsic connection 
between the spacetime and its content; the energy and matter content of the 
spacetime curves it while the spacetime tells the matter how to move. 

This interconnection is of geometric nature: the energy /matter content of 
the universe is closely related to the geometry of the spacetime. This is encoded 
in the Einstein equations: 

Gab^^Tab, a, 6 = 0,1, 2, 3, (1) 

where Gab is the Einstein tensor, a geometrical object related to the space- 
time curvature while Tab is the stress tensor, representing the matter/energy 
content of the spacetime; G Newton's constant and c is the speed of light. In 
particular, this confers a dynamical aspect to the spacetime, in opposition to 
the Newtonian theory where the spacetime is fixed and thus not dynamic. 

Einstein found that his equations predicted an expanding or contracting 
universe; but as he believed in a static universe anyway, he refused this option 
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and added a term in his equations: the cosmological constant term. By prop- 
erly tuning the value of this cosmological constant, he could recover a static 
universe. But a decade later, Hubble pointed out from the observations of the 
relative motions of stars that our universe should be expanding. Einstein finally 
abandoned the assumption of the cosmological constant. Note that nowadays, 
cosmological observations predict a small and positive cosmological constant. 
We will come back to this question later. 

The theory of General Relativity works quite well at the scale of the solar 
system and has passed numerous experimental tests, such as the displacement 
of Mercury's perihelion or the deviation of light rays. At larger scales how- 
ever, the observation of angular velocities of stellar objects in our galaxy is 
not compatible with the predictions of General relativity. Physicists generally 
postulate the existence of some dark matter, i.e. matter that we don't observe 
but that has an effect on the spacetime curvature. 

Black holes, which are also one of the main aspects of this thesis, are proba- 
bly among the most fundamental solutions to the Einstein field equations. The 
first example of black hole has been provided by Karl Schwarzschild p] in 1916 
on the Russian front during world war I. The Schwarzschild solution describes a 
vacuum spacetime geometry with a spherical symmetry. The Schwarzschild so- 
lution depends on an arbitrary parameter Af , asymptotically, the gravitational 
field behaves like the Newtonian gravitational field of a punctual object of mass 
M located at the origin. Note that the Schwarzschild solution also describes 
the region outside a neutral non-punctual static spherical massive source. 

The Schwarzschild solution is further characterized by a particular radius, 
the Schwarzschild radius. If a spherical massive object has a radius smaller 
than the Schwarzschild radius, it forms a black hole, i.e. an object so massive 
that massive or massless particles cannot escape from the gravitational field of 
this object, once they are too close. 

Black holes are characterized by an event horizon coinciding with the Sch- 
warzschild radius, it is the boundary of the region where the gravitational field 
is so strong that nothing can escape from it. In the coordinate system used 
by Schwarzschild, the event horizon corresponds to an apparent singularity in 
the solution; it is however possible to formulate the solution in a different co- 
ordinate system where the apparent singularity is removed. Black holes have a 
region of very strong gravitational attraction, which can be inside the horizon. 
They are therefore interesting from a theoretical point of view since they can 
reach the validity limit of General Relativity. Black holes are indeed sometimes 
called theoretical laboratories for high energy physics. We will come back on 
the question of black holes with more details in the end of the introduction. 

From a mathematical point of view, the number of dimensions of spacetime 
is just a parameter of the theory, although physically, we observe four dimen- 
sions. In the twenties, Kaluza and Klein [3] noticed that there was an elegant 
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way of formulating gravity and electromagnetism, assuming the existence of a 
fifth dimension. But since we don't see five dimensions, they assumed that this 
fifth dimension was very small, of the order of the Planck length (of the order of 
10~^^m). Unfortunately, their theory suffered from numerous problems, such 
as the existence of a bosonic fundamental field which did not correspond to 
anything physical at those times. 

A few years after Kaluza and Klein's theory, quantum mechanics was devel- 
oped and the idea of extradimensions was not developed further. The theory 
of Kaluza and Klein was then completely forgotten until recently when String 
Theories started to develop. 

What's in this thesis? 

This thesis is organized in five chapters (and an additional concluding chapter) . 
The first three chapters are independent of the last two. Here, we briefly 
summarize the content of the chapters and present some theoretical aspects of 
the models addressed. 

The first chapter is devoted to braneworld models, we review the main 
ideas of Kaluza and Klein regarding extradimensions. In particular, we intro- 
duce the Kaluza-Klein model in 5 dimensions. Next, we present the Kaluza- 
Klein reduction mechanism and turn to the Arkani-Hamed - Antoniadis - Dvali 
- Dimopoulos [H|3] model. Then, we introduce the first warped extradimension 
model, namely the Randall-Sundrum model [6l [7] . Finally, we present other 
types of branes available in String Theories, namely p-branes and _D-branes. 

Someone said extradimensions? 

As mentioned earlier, the idea of extradimensions is not new. It was pos- 
tulated by Kaluza and Klein in the '20s and forgotten until the early '80s, 
when String theories started to develop. Originally, String theory was 
designed to describe hadrons and strong interaction, but it turned out 
that the theory contains a spin-2 particle, resembling the graviton. It was 
therefore realized that String theory can provide a good candidate for a 
theory of quantum gravity. The main idea of string theory is to consider 
string-like instead of point-like fundamental objects. Whereas the formu- 
lation of the classical theories considering point-like objects minimizes the 
proper length of the worldlines of particles. String theory minimizes the 
area of world sheets, i.e. of the surface defined by the motion of a string 
along the proper time direction. 

However, the quantization of String theory is consistent only if the 
number of spacetime dimension is 26 (this is actually true for the bosonic 
String theory); or 10 for supersymmetric String theories [5]. Anyway, 26 
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(and 10) is much larger than 4; so the idea of extradimension is nowadays 
necessary. 

As we just said, the spectrum of String theories contains a spin-2 ob- 
ject, reminding of the graviton. Moreover, the low energy description of 
String theories includes gravitation. It then makes sense to study higher 
dimensional gravity, inspired by String theories. On the other hand, we 
also mentioned the fact that the theory of General Relativity can be formu- 
lated in more than four dimensions. Mathematically speaking, the number 
of dimensions is then just a free parameter. However, at our energy scales, 
we do actually observe only four dimensions. So if one deals with extradi- 
mensional theories, one has to ensure that the low energy description of 
his theory implies the observation of four dimensions. 

This is indeed the case for braneworld models, where our four dimen- 
sional universe is embedded in a higher dimensional spacetime. Either the 
extradimensions are small, so that they are not observable at our energy 
scale, either there exists a confining mechanism that sticks our observable 
world on the four dimensional slices of the entire spacetime. 

The second chapter is devoted to braneworlds with an extension along the 
extradimensions. The brane models considered in chapter two are described by 
the Einstein gravity lagrangian (with a cosmological constant term) extended 
by matter lagrangian describing bosonic fields. The matter lagrangian are 
chosen as usual classical fields theory known to admit 'soliton' type solutions. 
Generally speaking, a soliton refers to classical localized solution to nonlinear 
equations. It is well known that some field theories admit soliton solutions 
(cosmic strings in the case of the Maxwell-Higgs model, monopoles in the case 
of Yang-Mills- Higgs models). 

The philosophy of the approach of chapter two is to make use of these 
solitons by imposing the direction of localization of the solution in the extradi- 
mensions of the full model. 

This chapter is based on original results (SUTOlIII]; first we will present brane 
solutions to the Einstein-Maxwell-Higgs model with inflating four dimensional 
slices and a bulk cosmological constant (i.e. a cosmological constant for the 
entire spacetime) |10| . Then, we present solutions to the Einstein non abelian 
Higgs and Einstein- Yang-Mills-Higgs theory, again with inflating branes |10) . 
Finally, we consider the six dimensional inflating baby-Skyrme model jllj . The 
baby-Skyrme is a toy model for the Skyrme model somehow describing nucleic 
matter. The baby Skyrme brane model can therefore be viewed as a toy model 
mimicking a brane composed by nucleic matter. 
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What's the cosmological constant? 

Roughly speaking, the cosmological constant is an extra term in the Ein- 
stein theory of gravity, which leaves the equations consistent. An impor- 
tant issue of this thesis is the study of the effect of a cosmological constant 
on some classes of classical solution to Einstein equations. The Einstein 
equations with a cosmological constant read 

Gab — ^4 Tab — ^Sab, (2) 

where Gab and Tab are the same objects than in equation ([T]), A is the 
cosmological constant and gab is the metric. The metric is the essential 
unknown quantity of equation ([2|). It encodes the mathematical tool al- 
lowing to compute the distance between objects in spacetime. 

The cosmological constant has witnessed a strange history. It was first 
introduced by Einstein when he realized that his theory predicted a non 
static universe; either expanding, either contracting. The cosmological 
constant term was required in order to have a quasi-static universe, com- 
patible with the apparent staticity of the stars [12]. However, in the late 
'20s it was pointed out by Hubble that stars were actually not static, in- 
stead, they were getting further and further away from each other |13j . 
When Einstein heard about Hubble's discovery, he qualified the introduc- 
tion of the cosmological constant as 'the biggest blunder of [his] life'. 

It should be mentioned that the cosmological constant is a strange 
object by nature. The cosmological term can be viewed as part of the 
stress tensor, it leads to a constant energy density all over the universe 
and a constant pressure of opposite sign. In particular, the energy density 
p and the pressure P associated to the cosmological term are p = A and 
P = —A, leading to the equation of state p = —P. 

As already stated, current astrophysical observations point to a posi- 
tive cosmological constant. It is then interesting to see how the inclusion of 
a cosmological constant deforms solutions to Einstein equations obtained 
in the absence of a cosmological constant. On the other hand. String the- 
ories predict a correspondence between gravity solutions with a negative 
cosmological constant {AdS, see appendix [B]) and (conformal) field theory 
(GET), namely the AdS /C FT correspondence [T3]. There are then moti- 
vations to study the effect of both signs of the cosmological constant. 



The third chapter presents a simplified model for the coupling of fermions 
to topological defects (see [15]). This can be seen as a toy- model for the local- 
ization of fermionic fields on branes. This chapter is based on original results 
published in [16j . 
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Chapter four is a (non-exhaustive) overview of various black holes avail- 
able in higher dimensional general relativity. In higher dimensions, it turns out 
that the allowed horizon topology is much reacher than in four dimensions. In 
particular, there exists black holes with horizon topology Sd-2 as well as black 
strings with topology Sds x 5*1 (here Sp refers to a p-dimensional sphere) 
among many others. It should be mentioned that apart from their geometrical 
meaning, black holes can further be characterized by physical quantities which 
governed by laws strikingly resembling the laws of thermodynamics. A short 
overview of this is presented in the next section entitled black holes. As a 
consequence, thermodynamical properties of black objects can be emphasized 
and provide numerous informations about the physics of theses objects. 

In this thesis, we focus on the two simpler cases: black holes and black 
strings. First we present known solutions generalizing the electrically charged 
or rotating black hole solution in four dimensions to higher dimensions. It is in 
particular pointed out that {{A)dS) electrically charged and rotating black hole 
solutions are not available in higher dimensions, at least in an analytical form. 
We then present the numerical construction of such solutions, following the 
original results found in |17[ I18j . Afterwards, we study black string solutions. 
Black strings are solutions to the higher dimensional Einstein equations having 
the particularity that the horizon topology is not spherical but cylindrical. In 
the case of black strings, one of the spatial dimension of the spacetimes, say z, 
plays a particular role: it is assumed to be compact (or periodic); the axis of 
the axial symmetry coincides with the particular direction. As far as the other 
dimensions are concerned, hyperspherical symmetry is assumed. 

We review the basic asymptotically locally flat black string solution and its 
properties, in particular the fact that these objects are unstable [TO]. Then we 
study the influence of a positive cosmological constant on the equations of black 
strings. In particular, we provide numerical evidences for the non-existence of 
an asymptotically de Sitter black string (following the original results in jl7|). 
but instead of the existence of an asymptotically singular geometry. 

Next, chapter five is devoted to the study of asymptotically locally AdS 
black strings. In opposition to the case of a positive cosmological constant, an 
asymptotically locally AdS black string solution is available numerically and 
has been constructed by [50], extending the work of [3T] (for d ~ b). We first 
review the solution of [20] and its properties. In particular, we reconsider some 
of the thermodynamical aspects of the AdS black string and present a new 
phase of the latter, characterized by a negative tension [22]. Black strings can 
indeed be characterized by a tension, in addition to the mass. The black string 
solutions mentioned here are so-called uniform, in the sense that they don't 
depend on z. 

As mentioned above, asymptotically locally flat black strings are unstable. 
Since the AdS black strings investigated in this chapter depend naturally on an 



X 



INTRODUCTION 



supplementary parameter, namely the cosmological constant A, it is natural to 
investigate the stability properties of the latter as a function of A. The stability 
of these solutions can be studied both through the thermodynamical aspect as 
well as through the dynamical aspect, in the spirit of [H]. For the families of 
solutions considered, we were able to show that they have a dynamically stable 
phase and a dynamically unstable phase, agreeing with the thermodynamically 
stable and unstable phases, along with a conjecture due to Gubser and Mitra 
[23]. 

Then, we perform in detail the construction of the perturbative AdS non- 
uniform black strings and predict a new phase of thermodynamically stable 
non uniform black strings in AdS. By non uniform, we refer here to solution 
non trivially depending on the coordinate z. 

All the results accumulated in the perturbative approach are useful but 
nevertheless are approximations of the real non uniform solution, if exists. In 
the last section of chapter five, we finally consider the construction of the full 
(non-perturbative) solution. We provide strong numerical evidences that the 
non uniform AdS black string solution indeed exist and we foreseen that this 
family of solutions ends-up in an AdS localized black hole, whose construction 
was out of the scope of this thesis. In the asymptotically locally flat case, the 
order of the phase transition between uniform and non uniform black strings 
depends on the number of spacetime dimensions. We argue that in the presence 
of a cosmological constant, the counterpart of this critical dimension should 
depend on the cosmological constant. This chapter presents original results 
published in [Ml [IH [22] . 

Black holes 

Black holes are probably among the most puzzling objects in general rel- 
ativity. They are solutions to Einstein equations and they are governed 
by laws strikingly resembling the laws of thermodynamics. We will try to 
sketch the portrait of these amazing objects here. 

First, what is a black hole? Formally speaking, it is a solution to 
Einstein's field equations admitting a trapped spacelike region, i.e. a region 
of spacetime such that the future light cone of this region does not extend 
to infinity |26| . It formally describes the vacuum spacetime with a point- 
like massive source. In that sense, it is the general relativistic counterpart 
of the electric field of a point charged source in Maxwell theory. The 
solution is however valid outside a static spherical body of given mass M; 
for example, the Schwarzschild solution accurately describes the spacetime 
outside a (nearly) spherical object, such as the sun. 

More intuitively, a black hole is a massive object, so dense that the 
gravitational attraction prevents any objects which are too close to escape 
to infinity. They can be imagined considering Newton's law of gravity; the 



velocity required to escape the gravitational field of a body of mass Af at 
a distance R from this object is given by 



G being the Newton constant. It can be easily seen that if the distance 
between the observer and the body is smaller than 2GM/c?. c being the 
speed of light, the escape velocity is greater than c; i.e. the light itself 
cannot escape the gravitational field of the massive body. 

In the theory of general relativity, the simplest black hole of mass M 
is described by the Schwarzschild metric: 



where c?f2| is the square line element on a unit 2-sphere. 

There is clearly something special around the radial coordinate r = 
Rs = the sign of the time and radial metric coefficients change. 

Note that the surface of equation r = Rs is called the event horizon 
of the black hole. Roughly speaking, the radial coordinate becomes a 
time coordinate and conversely, for radial coordinates smaller than the 
Schwarzschild radius Rs- It follows that the future of an observer crossing 
the Schwarzschild radius is located at r = 0, since just before the crossing, 
the radial velocity is negative, dr/dt < 0, but just after the crossing, time 
and radius interchange and the new time, i.e. the radial coordinate r, 
flows towards r = 0. Note that a massive object is not infinitely dense. It 
follows that there is a critical density in order to form a black hole. From 
([3]) and (jl]), it is obvious that if the radius of the massive object is larger 
than Rs, it is not a black hole, since ((31) is a vacuum solution - valid only 
outside the massive body. 

Another question is: how do black holes form? There are essentially 
three types of black holes: stellar black holes, primordial black holes and 
micro black holes. 

Stellar black holes result from the collapse of massive stars. A star is 
formed by the accretion of dust and hydrogen under their own gravitational 
fleld. As the hydrogen collapses, the heat increases, until thermonuclear 
processes start consuming the hydrogen into helium. During this reaction, 
the radiation pressure of the thermonuclear processes counterbalances the 
self-gravitation of the star, preventing it from further collapse. But at 
some stage, all the hydrogen is burnt; the star then encounters a second 
stage of collapse, until thermonuclear processes start again, this time using 
helium as a fuel for the reaction. 
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These stages of radiation-collapse can continue until a stable element 
is reached, namely Iron. Then, there are essentially three possibilities; if 
the star is not massive enough, nothing special occurs any more and we 
are left with a white dwarf. The second possibility occurs if the star is 
more massive: the star collapses until the degeneracy pressure of the iron 
nucleus, due to the exclusion principle applied on the constituent of the 
star, stops the process. This is a neutron star. Finally, if the star was very 
massive, the collapse continues until the density reaches the critical density 
required to form a black hole. Stars leading to the formation of black hole 
are expected to have masses of the order 1.5 — 3 solar masses and greater, 
at the end of their life. This is known as the Tolman-Oppenheimer-Volkoff 
limit [211 Eg. 

The second type of black holes, primordial black holes, are believed to 
be formed in very early stage of the universe. Shortly after the Big-Bang, 
the density in the spacetime was extremely high. A small perturbation 
in the density might have resulted in the formation of a black hole, after, 
the medium collapsed under its own weight. The important point here is 
the presence of density variation; there are many models, many of them 
predict primordial black holes. 

Finally, micro-black hole are believed to be formed during high energy 
collisions. Consider two particles colliding, if the collision is such that 
these two particles become closer to each other than the value of the Sch- 
warzschild radius for the total mass of the particles, a micro-black hole is 
formed. Note that this would require ultra high energy. These processes 
are expected to occur for energies in the center of mass larger than the 
Planck mass where quantum effect are expected to break completely Gen- 
eral Relativity down. The Planck mass is about 10^^ GeV, micro black 
hole formation is then unlikely to appear, except in the case of extradi- 
mensional models where the fundamental Planck mass can be of the order 
of 1 TeV (see chapter [1]). 

Note that black holes are not just theorist's fantasy: there are strong 
observational evidences for a black hole to reside in the center of our galaxy 

m- 

Another amazing aspect of black holes is their link to thermodynamics. 
The laws governing black holes can be formulated according to three laws 
|301 [3TI [32] (we present them in units where G = c ~ h ~ 1): 

Olp: The horizon has constant surface gravity for a stationary black hole, the 
surface gravity being the gravitational attraction applied by the black 
hole on an hypothetical observer located at the horizon. 

1**: dM = -^dA + ildJ + ^dQ, where M is the mass of the black hole, k the 
surface gravity, A the horizon area, il the angular velocity, J the angular 
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momentum, $ the electrostatic potential and Q the total electric charge 
of the black hole. 

2"'^: dA > 0. 

The zeroth law reminds the zeroth law of thermodynamics stating that 
temperature is constant inside a body in thermal equilibrium, suggesting 
a correspondence between surface gravity and temperature. 

The first law is similar to the first law of thermodynamics, expressing 
the conservation of energy. In General Relativity, mass and energy are the 
same thing. The first law then states that the energy variation of a black 
hole comes from the variation of rotational energy, electrostatic energy and 
variation of the horizon area. Usual thermodynamics include the variation 
of entropy, TdS in the first law. It is then tempting to relate the horizon 
area and the entropy of a black hole. This is comforted by the second law, 
analogue to the second law of thermodynamics stating that the entropy 
can only increase during physical processes. 

These ideas have actually been derived more formally by Hawking, 
Bekenstein and Carter, leading to Th = k/{2tt), S = A/A, Th being the 
black hole temperature, S being its entropy. 

The sixth chapter is the final chapter of this thesis where we review the 
main results presented in the text and give some possible outlook and perspec- 
tives of our work. 

Note that natural units will be used throughout this thesis: h = c = 1. 
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Chapter 1 

Braneworld model 



In this chapter, we wih present some brane models appearing in theoretical 
physics. Braneworld find their existence in extradiniensional theories (say with 
d dimensions); what is called the brane is then a four dimensional subspace of 
the d dimensional spacetime, describing the four dimensional spacetime where 
we live. The idea of extradimension is not new and was first introduced by 
Kaluza and Klein in the beginning of the last century [2l[3]. However, the idea 
of Kaluza and Klein was motivated by unifying electromagnetism and gravity 
as we will see later, but has been forgotten for decades, until String theories 
started developing. 

String theories are consistent in more than four spacetime dimensions, de- 
pending on the model and furthermore include general relativity as a low energy 
description. This is part of the motivation for considering general relativity in 
more than four dimensions. However, in order to explain the fact that we in- 
deed observe four dimensions, one has to find a mechanism that leads ordinary 
particles to evolve in a four dimensional spacetime, just like an ink drop would 
live in the two dimensions of a sheet of paper, although the sheet itself lives in 
three spatial dimensions. 

Among other, one of the first viable brane scenario was introduced by An- 
toniadis, Arkani-Hamed, Dimopoulos and Dvali in 1998 |11[5] (AADD), where 
the authors consider a (71 + 4)-dimensional spacetime with Ricci flat extradi- 
mensions. They were able to give an explanation to the hierarchy problem, i.e. 
the huge difference between the order of magnitude of the Planck mass and the 
electroweak uniflcation mass scale. This is indeed one of the features of brane 
models. Their model excluded n ~ 1 and put strong constraints on the size of 
the extradimensions. 

Later, in 1999, another mechanism has been proposed by L. Randall and 
R. Sundrum [5J [7], where the authors consider a flve dimensional spacetime 
and a thin brane, described by a Dirac delta. The Randall-Sundrum model 
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includes two branes: one is the brane on which we are hving, the other is a 
'mirror', phantom brane. The length of the extradimension can be infinitely 
large thanks to some warp factor. 

This chapter is organized as follows: first, we will review the Kaluza-Klein 
model and the so-called Kaluza-Klein reduction mechanism in the first sec- 
tion. Then, we briefly review the main features of the AADD brane model 
before turning to the Randall-Sundrum models. Finally, we will briefly present 
p-branes, D-branes and black p-branes in (low energy) String theory for com- 
pleteness, but also for the fact that they somehow link branes and black objects 
presented in chapters 2] and [S] of this thesis. 



1.1 Kaluza Klein model 

The original motivation of Kaluza [5] and Klein [31 was to give a unified formu- 
lation of gravity and electromagnetism. The construction is the following (see 
[33] for a review): consider a five dimensional metric 

ds^ = gaidx^-dx", a, 6 = 0, 1, 2, 3, 4. (1.1) 

This five dimensional metric can be written in a factorized form, without loss 
of generality: 

ds^ ^ e-^/^ (g(,t^ + K'^e-^'l'A^A^dx^'dx'' + 2^-^"^ A^dx^'dy + e^^-^V^ , 
/i,i/ = 0,l,2,3; (1.2) 

where y denotes the fifth dimension, y g X C M. The metric component g\^} 
is a tensor, A^ is a vector and is a scalar from the four dimensional point of 
view, K is a real constant. 

Assuming that the fields g)S , A^^, do not depend on the extra coordinate y, 
the five dimensional Einstein-Hilbert action reduces to [Mj 

S = / (1.3) 

IbTTG Ja^xI 

where is the four dimensional submanifold of the 5-dimensional spacetime, 
Fp,u = ^^J,Ay — dyA^j, is the field strength of A^ and Ge// ~ G j l^^dy is the 
effective Newton constant, G being the 5-dimensional Newton constant, i?^*) 
is the Ricci scalar computed with g^^y . 

Setting (/) = Const., the action reduces to the Einstein-Maxwell model, 
providing a unified model of gravity and electromagnetism. Unfortunately, this 
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attempt of unification has been forgotten because of the success of the quantum 
electrodynamics theory. The idea regained interest only in the '80s when String 
Theories brought the question of extradimensions up to date. Anyway, it is 
remarkable that compactifying one dimension, reducing the model from five 
dimensions to four leads to the occurrence of a gauge invariant theory, the 
gauge group being U{1). It should be stressed however that (j) = Const, is not 
a solution of the 5 dimensional model because of the non trivial coupling with 
the Maxwell field. It is anyway remarkable that 5 dimensional gravity somehow 
contains Einstein-Maxwell (-dilaton) theory 

This procedure is known the Kaluza-Klein compactification; note that com- 
pactification from more than five dimensions to four dimensions induce other 
gauge groups, usually non abelian, depending on which manifold the extradi- 
mensions are compactified. 

1.1.1 Kaluza-Klein reduction 

Another mechanism enters the game once one deals with compact extradimen- 
sions: the Kaluza Klein reduction. Consider a scalar field F{x^,y) in five 
dimensions, where the extradimension is compact and Ricci-flat: 

ds^ = gl^Jdx'^dx-' + dy\ y G [0,i], (1.4) 

for some real L. 

The five dimensional massless Klein- Gordon equation leads to 

□5F(a;^ y) = □4^^(x^ y) + d'^F{x^, y) = 0. (1.5) 

where Ds is the five dimensional d'Alembertian while ^4 is the four dimen- 
sional d'Alembertian. In order to derive a four dimensional equation, we 
can perform a Fourier series decomposition in the y direction according to 
F{x^,y) = E„/n(a;'')e^''™2'/'^, n&'L. Equation then reads 

(□4-M2)/„(a;'^) = 0, (1.6) 

which is the equation of a four dimensional scalar field with mass M„ = 
2'iT\n\/L. In other words, a massless scalar field in five dimensions appears 
as a tower of massive scalar fields in four dimensions. 

This procedure can also be applied to vector fields or fermionic fields but is 
more involved, so we won't present it here (see for example |35j). 

The Kaluza-Klein reduction furthermore generalizes to higher number of 
dimensions, the result being again that massless fields in higher dimensions 
appear as a tower of massive fields from the lower dimensional point of view. 
Note that this is also the case for massive fields in higher dimensions, except 
that the mass tower is shifted by the higher dimensional mass of the fields. 
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1.2 A ADD braneworld 

The Einstein-Hilbert action in n + 4 dimensions has the form 

S = M;+2 J ^R„+^d^xd^y, (1.7) 

where Rn+i is the n + 4 dimensional Ricci scalar, Alp the d-dimensional Planck 
mass and g is the determinant of the metric. 

Now, consider a n + 4 dimensional spacetime of the form 

ds^ = r^^^ydx^'dx" + S.jdy'dy^, tj = diag(-l, 1, 1, 1), (1.8) 

where = 0, 1, 2, 3, i, j = 1, . . . , n, (5,;j is the Kronecker symbol and where we 
assume the extra coordinates to be in [0,ii]. The metric (jl.8[) is a vacuum 
solution to the n + 4 dimensional Einstein equations. 

The action S can be reduced to a four dimensional effective action. The 
intuitive procedure is to integrate (jl.7p over the n extra coordinates, assuming 
that the Ricci scalar does not depend on these extra coordinates and that the 
geometry of the extradimensions is given by (|1.8|) , leading to 




i?4 being the four dimensional Ricci scalar, 34 the determinant of the metric 
on the four dimensional slices and the volume of the extradimensions. This 
procedure implicitly assume that the four dimensional part of the metric (jl.8p 
can fluctuate, but not the extradimensional part. 

Identifying the coefficient of the integral (|1.9p as the square of the four 
dimensional Planck mass A^p, we have the following relation between the n + 4 
dimensional and four dimensional Planck mass: 

Ml = m;+V„. (1.10) 

In other words, the four dimensional Planck mass appears as the prod- 
uct of a power of the 71 + 4 dimensional Planck mass and the volume of the 
extradimensions . 

As a consequence, it is possible to obtain a very large 4 dimensional Planck 
mass with a fundamental Planck mass of the order of the TeV (which is the 
electroweak unification energy scale), provided that the volume of the extradi- 
mensions is large enough. 

A straightforward consequence of this class of models is the deviation from 
Newton's gravitational law. Assuming the length Li to be all of the same order 
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of magnitude R, Gauss' law in A + n dimensions implies that the gravitational 
potential felt by two test masses mi, m2 and separated by a distance r is 



Note that this is consistent with (|1.10p : the effective Planck mass is the fun- 
damental one times the volume of the extradimensions. 

As a consequence, one can assume the fundamental A-\-n dimensional Planck 
mass to be of the order of the TeV, along with the electroweak unification scale 
mEw- The large discrepancy between the Planck scale and the electroweak 
scale is then a geometrical effect, due to the volume of the extradimensions. 
Assuming so puts constraints on the radius R [?]: 



Note that 71 = 1 is directly ruled out, implying deviation from the Newton's 
law at astrophysical scale (for rt = 1, i? « lO^^m ). However, from n = 2, the 
deviations would appear at the millimetric scale and beyond, compatible with 
experimental tests of Newton's gravity law 137j . 

This class of models provides an elegant solution to the hierarchy problem. 
However, the geometry of the extradimensions is completely separated from the 
geometry of the 4 dimensional branes. Moreover, these models rely on compact 
extradimensions; there is no natural reason that only three spatial dimensions 
are infinite and all others are compact. These geometries, where the four 
dimensional spacetime is completely separated by the extradimensional part of 
the spacetime are called factorisable geometry. Releasing the assumption of 
factorisability of the spacetime leads to other classes of models where the size 
of the extradimensions can be much less constrained. It is indeed the case of 
the Randall-Sundrum model [6]. 

1.3 Randall-Sundrum braneworld 

Although the AADD scenario sketched above might solve the hierarchy prob- 
lem, another mechanism has been introduced in 1999 by Randall and Sundrum 
in order to solve the hierarchy problem [6l[7], without requiring compactifica- 
tion. Their approach is based on the use of warped spacetimes: the advantage 
of such spacetimes is that the warp factor can allow a large volume of the ex- 
tradimensions without constraining too much the range of the extradimensional 
coordinates. 




(1.11) 




(1.12) 
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We will briefly review the solution and proposal. The model considered in 
[H [7] is described by the following action: 

Sgr = I d'^x I d(f>y/^ {2M^R - A) , 



Shid = J d xy/~ghid {Lhid - Vhid) ■ (1-13) 

where g is the determinant of a five dimensional metric while g^s , ghid are the 
determinant of the induced metric on a visible (resp. hidden) brane, i.e. four 
dimensional subspaces where matter fields live, embedded in the five dimen- 
sional space (the bulk). M is the five dimensional Planck mass, L^is (resp. 
Lhid) is the lagrangian describing the matter fields on the visible (resp. hid- 
den) brane, Vyis (resp. Vhid) is the vacuum energy in the visible (resp. hidden) 
brane. 

The metric ansatz is given by 

ds^ = e'^^'^'ry^^dx^da;'' + r^d^^, g [-tt, tt], (1.14) 

supplemented by a Z2 symmetry. This is typically a warped geometry and 
denotes some compactification radius, which can be infinite. The visible brane 
is placed in = while the hidden one resides at = tt. 

In such a setup, the solution to the Einstein equations is given by 



Note that although this function is continuous, this is not the case for its 
derivative. This is due to the fact that the branes don't have an extension. 
We will see in chapter [2] that considering branes with an extension regularizes 
the metric functions. Furthermore, the bulk cosmological constant A has to be 
negative. 

The result of Randall and Sundrum is such that the parameters of the 
model, namely the vacuum energy of the visible and hidden brane Vvis , Vhid 
and the 5-dimensional cosmological constant have to obey the following relation 

K»s = -VMd = 2AM^k, A = -24M3fc2, (1.16) 

for a given value of fc € M. 

The effective four dimensional Planck scale in this model is given by 

M'^^^{l-e-^'^^-), (1.17) 
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allowing Tc to be large, the control on the four dimensional Planck mass being 
essentially provided by the value of k. 

Let us review how Randall and Sundrum found this result: consider a mass- 
less fluctuation of the metric: 



where is a massless tensor fluctuation (there is no dependence in the ex- 
tradimension, see section fl.ip while T(x) is a massless scalar fluctuation. The 
authors argued that there shouldn't be off-diagonal vector fluctuation in the 
low energy models, since these vector modes would be massive [HIT]. The ef- 
fective four dimensional action resulting from this model contains a term of the 



where R is the Ricci scalar constructed with p'^j^ and doesn't depend on (p. 

This leads to the effective four dimensional Planck mass presented in equa- 
tion (fTTTl . since the only term depending on the extra-coordinate in (jl.l9p 
is g-'2''r<:\'l>\ ^ which can be explicitly integrated. This provides a possible so- 
lution to the hierarchy problem without putting constraint on the length of 
the extradimension: the effective Planck mass is still well defined in the limit 
where the hidden brane is pushed to infinity. Note that the construction of the 
effective action is formally equivalent to integrating out the extradimensional 
dependence of the 5-dimensional action. 

This model is quite elegant and opened a new research area since it was the 
first model to consider warped branes. However, it contains some criticizable 
points: first, as already mentioned, the derivative of the metric function is not 
continuous, second, there is a fine tuning between the parameters of the model 
{VvisjVhid and A). Third, it requires the existence of a mirror hidden brane. 
Note that the third point has been however discussed in [7]. 

1.4 Branes in String theory 

For completeness, we introduce the notion of p-branes solutions and D-branes 
in this section. These two objects are different in their nature, but are conjec- 
tured to be actually equivalent, leading to the basis of the so-called AdS/CFT 
duality, which relates gravity in AdS and conformal field theory defined in the 
background of the conformal boundary of the AdS spacetime. We will not 
enter the details of this duality, neither detail the construction of the effective 
models where p-branes live. We will actually not even give much details about 
how the various object entering the discussion appear theoretically; instead we 
will present the model and p-brane solution and explain briefly what are the 



form 
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D-branes. We refer the reader to |381 [39l |40] between many others for more 
details on the construction of the theory. 

In superstring theory, the fields describing the strings are bosonic and 
fermionic. The bosonic strings as well as fermionic strings can be closed or 
opened; accordingly suitable boundary conditions are imposed; in type II string 
theories (which is relevant for our purpose), the strings are indeed closed. 

Nevertheless, fermionic strings admit two different kinds of boundary con- 
ditions: periodic (Ramond: R) or anti-periodic (Neveu-Schwarz: NS) (see for 
example [S]). 

Depending on the type of boundary conditions imposed the spectrum of 
the string is slightly different. Type // string theories consider closed string; 
for example, the NS — NS sector provides a fundamental 2-form B while the 
R — R sector provides {p + l)-forms C^"*"^, with p odd or even, depending if 
the string theory considered is of type IIA or type IIB (this is related to the 
chirality of the fermionic sector of the strings). 

The idea here is to consider the bosonic part of low energy type II super- 
string theory. 

So finally, the bosonic part of low energy type II string theories contains 
general relativity, dilaton and p-forms action: 



where gs is the determinant of the metric gs,ab, 4> is the dilaton, B is the 
NS — NS 2-form, C^+^ are the R — R (p + l)-forms; d is the exterior derivative 
and the squares are to be understood as contractions on the indices of the form 
components. This action is the bosonic part of the 10-dimensional supergravity 
[41) . which is a good approximation of type II string theories at moderate 
energies. 

Note that gs,ab is the metric in the so-called 'string frame'; it is possible 
to re-express this action in the 'Einstein frame', defining ge,ab ~ e~'^^^5s.ah! 
leading to 



In the Einstein frame, (|1.20|) reduces to the Einstein-Hilbert with the Klein- 
Gordon action for the dilaton, a 'Maxwell-like' actions for the forms and cou- 
pling between the dilaton and the various forms. 




(1.20) 




(1.21) 
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1.4.1 p-branes 

In order to construct the j>-brane solution, we will consider a truncated action, 
taking into account only one (p + l)-form and setting B to zero: 

7l0„ 



Sii,p = I d'''x^^,{[Rig,)~Ada(l)d^ 
1 1 



2(p + 2)! 

where Up = (3— p)/2. The equations of motions resulting from the action (|1.22p 
are given by 

Rab = ^da(t)db(j) + Sab, 

o _ ^ ^".0 ( p p m2...mp+2 _ _P±J_„ . p2^ 

bab " 2(p+l)! |^^-"2...m,+2n 8(^+2)^''"' y' 

Q — |^gap0^am2...mp+2~J ^ 

□ ^ ^ °P ^, e°'"^F^ (1.23) 

^ 2(p + 2)! ' ^ ^ 

where F = dC'P+^ . 

There exists solutions to the equations (jl.23p having the properties that 
they extend in p directions. These solutions are known as p-branes. For con- 
venience, we will split the 10 coordinate in two parts: p + 1 longitudinal 
coordinates a;^, p, = 0, . . . ,p and 9—p transverse coordinates y\ i = p+1, . . . , 9. 
Then, the solution is given by 

ds^ = H^P-^'^'S^..dx''dx'' + H^P+^^'^5,jdy'dy', 

20 _ ff{3-p)/2 
c Hp , 

CIS!:j - H-'-l, (1.24) 
The function Hp is a function of r = \/^. and is given by 

Hp = l + ^, (1.25) 

where Qp is the charge associated to the {p + l)-forni (see [41] and references 
therein) . 

Note that in the string frame, the line element takes the simpler form 

ds^ = Hp^/^r^p^dx^'dx" + H^^S.^dfidy^ (1.26) 

It is interesting to note that low energy description of String theory provides 
extended solutions. In other words, these theories predict brane solutions. 
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1.4.2 D-branes 

On another hand, String Theory predicts also another type of branes, namely 
D-branes. D-branes consist actually on hypersurfaces where opened strings 
end. These open strings are subject to Dirichlet boundary conditions, giving 
the D to D-branes. These objects have many properties, and can be put in 
relation with the p-brane solution presented above, but this is much too far 
from the object of this thesis. 

1.4.3 Black p-branes 

Let us finally remark that there exists extended objects in higher dimensional 
gravity which are also solutions to supergravity described by (|1.22[) , referred to 
as black p-branes |42]. These black p-branes can carry some charge associated 
to the p-forms present in (|1.22|) . However, we will exhibit simpler solutions, 
uncharged with respect to the p-forms: 

ds^ = -f{r)dt' + — + r^dnl_j, + S,,dy'dy^ , j = 10 - p, . . . 10, 

= Const., CP+^ = 0. (1.27) 

where = 1 — (-7^) ^, ^ M+. This is simply a (10— p)-dimensional black 
hole with p transverse flat direction. The construction of the solution with non 
trivial p-form and dilaton can be found in j42) . 

1.5 Concluding remarks 

There are many ways to approach branes; for instance, we can consider pure 
gravity with extradimensions and add some branes to the model, we can con- 
sider low energy String theories, full String theory, etc. Actually, once one 
wants to deal with extradimensions, one has to introduce the concept of branes 
in order to have a chance to explain why we actually don't see the extradimen- 
sions. It seems clear that branes play a very important role in actual theoretical 
physics. Note that extended objects also play an important role in other fields 
of physics, such as surface physics in condensed matter |43j as an only example. 

Another important remark is that in the Randall Sundrum model, the brane 
doesn't have an extension in the transverse space while p-branes are somehow 
extended in the direction transverse to the brane. Recall that the metric func- 
tions were not regular, in the sense that their derivative was not continuous, but 
the metric functions somehow regularize in the p-brane model, the extension 
of the p-brane removing the discontinuous behaviour of the metric functions. 

Note also that the p-brane solutions (|1.27p are actually black objects. In 
that sense, they are a nice junction between the considerations of chapter [5] 
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and [3] of this thesis, deahng with extended brane models and the rest of this 
thesis, deahng with higher dimensional black holes and strings. 
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Chapter 2 



Topological braneworld 
models with extended branes 

In this chapter, we will present four brane models where the brane has an 
extension in the extradimensions. The idea for providing an extension to the 
branes in the extradimensions is to use localized soliton solutions available in 
usual field theory and to interpret the dimensions where the soliton lives as 
the extradimensions of the spacetime. We will consider models containing two 
cosmological constant: one is a cosmological constant for the entire spacetime 
(the bulk cosmological constant) while the second is a cosmological constant in 
the four dimensional subspace of the entire spacetime, i.e. on the brane. The 
second cosmological constant will be modelled by an inflating four dimensional 
subspace, i.e. an inflating brane. 

The aim of this chapter is to study the influence of these two cosmologi- 
cal constant on the pattern of some solitonic brane models considered before, 
without cosmological constant and inflation (for instance [44 l [45 t l46]). 

The first model we will consider is the Einstein Abelian Higgs model [9] in 
4 + 2 dimensions, the second model is the Einstein non abelian Higgs model 
|10| in an arbitrary number of dimensions, the third model is the Einstein- 
Yang-Mills non abelian Higgs model in 4 + 3 dimensions and the fourth model 
is the gravitating baby Skyrme model [47] in 4 + 2 extradimensions [11]. In all 
these models, the brane is seen as a topological soliton extending in n = d — 4 
extradimensions, d being the total number of dimensions. 

In all cases, the matter fields are localized in the extradimensions, defining 
the region where the brane is located. The four models we consider all have 
solutions with non trivial topological properties. This is the reason why we call 
these brane models topological brane models with an extension. 

In addition, as mentioned before, the four dimensional slices are chosen to 
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be inflating, modelling a four dimensional positive cosmological constant. This 
is motivated by the fact that the actual observations points to a four dimen- 
sional positive cosmological constant [351 1321 • Note that this is also relevant 
for inflationary models with extradimensions. 

The general action for all the models we will consider has the following form 

S S gravity ^brane^ (^■-^) 



where Sgrav is the Einstein-Hilbert action: 



^—jd''x^g{R-2Kd), (2.2) 



S gravity - 



Ad being the bulk cosmological constant, Gd is the d-dimensional Newton con- 
stant, related to the d-dimensional Planck mass Md by Gd = 1/M^~'^ and g 
the determinant of the c?-dimensional metric. 



The term Strane in (|2.ip is the action of the matter fields where the brane 
resides and will depend on the model under consideration; we furthermore 
define Lbrane as the lagrangian density such that Sbrane = J \/—gLbraned'^x. 



2.1 The metric ansatz and Einstein equations 

We consider a d-dimensional spacetime, with 4 special dimensions describing 
our universe. There are then n — d — 4 extradimensions (or codimension). 

The ansatz for the d-dimensional metric reads: 
ds^ = AP{p) [~dt^ +e^"*5^Jdx'dx'] + dp^ + l'^{p)dnl_^, i,j = 1,2,3, (2.3) 



where p and 0a-, a — 1, . . . ,n — 1 are the coordinates associated with the extra 
dimensions, drt^_r^ = J22=i Y[b<a ^^'^^ ^bdO^ is the square line element on the 
d — 5 sphere, di e [0,TT],da^i S [0,27r], i,j = 1,2,3 and x^,x'^ and x^ are the 
3-dimensional spatial coordinates. iJ > is the Hubble parameter related to 
the (positive) 4-diniensional cosmological constant. 



2.2. VACUUM SOLUTION 



15 



In the general case, the Einstein equations read 



2 \P P 

= SnGdTi: - Ad, ^ = 0, . . . , 3 

" 2 A/2 + ^ M2 + Ml ^ 2 ^ /2 p J 

(n-2)(n-3) (l^ 
^ 2 lv"^"/2 

= %^GdTll - Ad, i = 1, . . . , n - 1 (2.4) 

where n = d — 4 is the number of extradimensions, Gab is the Einstein tensor 
and Tab is the stress tensor (see Appendix [A|) of the matter fields described by 
Sbrane and i?'"*^ is the scalar curvature of the four dimensional slices; in our 
case, it is 12H^. 

Note that the brane appears only thought i?'^' in the equations; it follows 
that we can replace the four dimensional subspace by any spacetime with a 
constant positive curvature. 



2.2 Vacuum solution 

In order to consider vacuum solutions, we set Tab = in (|2.4p . According 
to the sign of the 4 + n-dimensional cosmological constant, vacuum solutions 
have constant positive, null or negative scalar curvature. The result is very 
similar to the 4-dimensional case where the geometry of the spacetime can be 
opened, closed or flat according to the sign of the cosmological constant. In 
the present case we emphasize a more general situation where 4-dimensional 
slices of the space have a de Sitter geometry, which can induce angular deficits 
in the n-dimensional subspace, as we will see later. Because the equations do 
not explicitly depend on the radial variable p the solutions given below can be 
arbitrarily translated in p. 

In this section we will consider solutions for n > 2 for reason that will 
become clearer later. We will come back on special solutions for n = 2 in 
section 12.3.31 

In the case where H ^ 0, the Einstein equations above possess explicit 
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solutions, depending on the sign of the bulk cosmological constant A^: 



M{p) = — i —^sinujp, l{p) = -ih^^—^smujp, foi Ad > 0, 
oj y n + 2 uj y n + 2 



In — 2 



M{p) = — W— ^sinhwp, = -W^^—^sinhwp, for Ad < 0, 

u y n + 2 u) y n + 2 

2 _ 2|A,| 

" 12 + (n-l)(n + 6)' ^^'^^ 

Note that these solution depends crucially on the two cosmological constants 
and from the form of the / function, it is clear that n = 2 is a special case. 
Indeed, in the 6-dimensional case where n = 2. the equation for M decouples 
[9] and leads to I cx M', which is not compatible with the solution above. The 
case n = 1 is also special since for 5 dimensions, the function l{r) is not defined. 

Note also that the solution (|2.5p are not regular since the Kretschmann 
invariant K = RabcdR"'^'"^ 

24i/4 48H'^M'^ 4(n-l)Z"2 16{n-l)PM'^ 2{n - 2){n - 1)1"^ 
^ ^ ~Af^ M4 + p + + jpfi 

4(n - 2)(n - l)^'^ 2{n -2){n- 1) 16M"^ 24M'^ 

¥ ^ ^4 + M2 + 

evaluated with the above solution gives 



(2.6) 



•J \ [n — 2)sm(h) [ujp) I 

where we define sin(h)(a;) = sin(a;) for A^ > 0, — sinh(x) for A^ > and x for 
Ad = 0. The invariant (|2.7p is obviously singular at the origin. Moreover, in 
the case A^; > 0, it is also singular for r — kn juj {k integer). 
The Ricci scalar, given by 



AP AP J \ AI ^ ' All 

is constant for these solutions (it is proportional to A^, using the equations of 
motion). Let us finally emphasize that the warp factor Al is directly propor- 
tional to the parameter i7, both are related to the four dimensional inflating 
subspace... 
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The solution p.5p has a natural geometric interpretation: it describes the 
surface of an n + 1 dimensional manifold of constant curvature (sphere, hyper- 
boloid or plane) in the extra dimensions and presents an angular deficit relative 
to the angles 9i. 

The angular deficit is computed in the following way: the extradimensional 
part of the line element with the solutions (|2.5p reads 



n — 2 1 

dsl = dp' + —-—sm{hf-iojp)dnl_,. (2.9) 
n + Z ui'^ 



Setting O — top, we find 



ds„ = 7^ ( de' + '-^—^ sin{hfQdnl_^ ) . (2.10) 



From this expression, it is clear that the factor induces an angular deficit 
1 — in the angular direction. Note that the angular deficit vanishes in the 
limit where n oo. 

In addition, the radius of the constant curvature manifold is found to be the 
parameter \/oj defined in p.Sp . In the case of Ad > 0, where the extradimen- 
sions present a closed geometry, it defines naturally a compactification radius 
for the extradimensions. 



2.2.1 Melvin-type universe in d dimensions 

In this section, we look for solutions of the form: 

M{p) = Moip - poT, l{p) = h{p - Por, (2.11) 

where Mq, Iq, pq, p,, a are constants to be determined. 

In the following, we will distinguish the cases where the solutions develop 
the above behaviour for p — >■ oo (asymptotic solution) and the case where the 
solutions develop a singularity in the neighbourhood of p = po for some real po 
(see section 12.51 and reference [SD] for the n = 3 case) , respectively. 



Asymptotic solution 

Because of the occurrence of non-homogeneous terms (e.g. 1/P and H'^ /AP) in 
the Einstein equations, power-like solutions of the form above cannot be exact 
for generic values of H, Ad, n. 

However, we will see that solution of the form (|2.1ip as asymptotic solutions 
and appear as 'critical' solutions when matter fields are supplemented in the 
form of global and local monopoles (see Section and 1^75)) . 



18 



CHAPTER 2. EXTENDED BRANES 



Let us for a moment neglect the non-homogeneous terms in the Einstein 
equations. Inserting the power law above in the Einstein equations leads to the 
following conditions for the exponents fj,, a: 

in — l)fn — 

3//(Ai-l) + 3Ai^+3(n-l)^a + (n-l)a(a-l) + -^ '-^ —a^ =0 (2.12) 

6/^2 + 4(71 - l)//a + ^^^^^^^^^a^ ^0 (2.13) 
V(Ai-l)+6Ai' + 4(n-2)Aia+(n-2)a(a-l)+-^^^— ^fc^a^ ^ q ^2.14) 

The solutions then read: 

2± V(n + 2)(n-l) (n - 1) T 2v/(n + 2)(n - 1) 

u = — ; r , a = -. ——. ^ 2.15) 

' 2(n + 3) ' (n-l)(n + 3) ^ ' 

Note, however, that with these exponents, it is not justified to neglect the 
inhomogeneous terms (x l/l"^ and oc 1/M^ except in the particular case of the 
critical solutions (see Section [2. 4p . where such terms vanish. 



Singular solutions 

If we want to interpret the functions (|2.1ip . (|2.15p as the dominant terms of a 
solution of the vacuum Einstein equations which is singular in the limit p — > po, 
the exponents should be such that fi < 1, a < 1 (along with the assumption 
that the inhomogeneous terms are sub-dominant). It turns out that these 
conditions are fulfilled for both values of the sign ± in (|2.15p . 

Note the relation between fi and a: Afi + (n — l)a = 1. This reminds the 
Kasner conditions [5T] : 

4/i2 + {n- l)a^ = 1, Afi + {n- l)a = 1, (2.16) 

except that only the linear relation is fulfilled. We will refer to this type of 
solution as Kasner type. 



2.3 Einstein abelian Higgs Model 

In this section, we consider a six dimensional spacetime with matter fields 
described by the Abelian-Higgs model. 

The action Sbrane for the Einstein- Abelian-Higgs (EAH) string is given in 
analogy to the 4-dimensional case [52l [53] by: 

S,rane = J <f (^-^Fm^F^^^ - \ (Dm^) {D'' <f)* " ^(0> " v'f^ 

(2.17) 
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where Dm = Va/ — IcAm is the gauge covariant derivative, Fmn = QmAn — 
On Am, M,N = 0, . . . ,5 the field strength of the U(l) gauge potential Am, e 
is the gauge coupling, v the vacuum expectation value of the complex valued 
Higgs field (p and A the self-coupling constant of the Higgs field. 



2.3.1 The ansatz 

The ansatz for the 6-dimensional metric is given by p.3p with n = 2. Note 
that here, 6i G [0,2tt]. We will denote the angular coordinate 6 since there is 
no possible confusion with other angular variables here. 

The ansatz for the non vanishing gauge and Higgs field reads |52| : 

e) = vf{p)e'^\ Ae{p, e) - -{N - P{p)), (2.18) 

e 

along with the ansatz of the Nielsen- Olesen cosmic string |52], and where TV is 
the vorticity of the string, which throughout this section will be set to iV = 1. 

Note that this ansatz is symmetric under rotations in the two extradimen- 
sions. 



2.3.2 Equations of motion 

Introducing the following dimensionless coordinate r and the dimensionless 
function L: 

r = VAwp , L{r) = VXvl{p), (2.19) 

the set of equations depends only on the following dimensionless coupling con- 
stants: 

a=T ' 7'-8^G6«2 , , (2.20) 

A Au^ Au^ 

The gravitational equations then read 



M" L" L' M' M'2 . K ^ f (l-P? 



2L2 2aL2 



(2.21) 



M'2 L' M' . K ■? f (l- P? PP^ P'^ 



L M AP ' V 2 4 2L2 2aL^ J ' 

6— ^ -h4— + A- 2— = -7^ ' . y J > J \ 



M"^ M A'P V 2 4 2i2 2aL2 
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the equations for the matter fields read 



+ = 0, (2.22) 



- af^P = 0, 



the prime denoting tlie derivative with respect to r. For later use, we note that 
these equations are solved by / = f , P = {). We will refer to this configuration 
as the vacuum solution. 

The equations ([2:2T|) can be combined to obtain the following two differen- 
tial equations for the two unknown metric functions; 



M^L 2 2 \2aL^ 4^ ' 



(2.23) 

12 



2.3.3 Six dimensional vacuum and asymptotic solutions 

Explicit solutions to equations p.23p can be constructed for /(r) = f and 
P(r) = 0. These are by themselves of interest, of course, but are also interesting 
for the generic solutions since we would expect that the metric fields take the 
form presented below far away from the core of the string. A similar analysis 
has been done for global defects in |54| . 

The solution for M, L can be written in terms of quadratures as follows: 



3"^" 10" 



r — ro — / dAI . / 7 , L — Ln — Ln\i 

(2.24) 

where fo, C are integration constants and Lq is arbitrary. 

Unfortunately, we could not find a general solution to the integral (|2.24[) . 
We will now discuss some particular cases where it is indeed possible to perform 
the integral explicitly. 

Static branes {k = 0) 

For A = the system admits two different types of solutions |551 156) : 

Ms = l , Ls = r-ro (2.25) 



and 



Mm ^Ci{r- ro)2/5, Lm = C2{r ~ ro)-^/^ 



(2.26) 
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where Ci, C2 are constants related to C and Lq. 

By analogy to the case of section 12.2.11 we refer to the second type of 
solution as the 'Melvin' (M) branches. The first type of solution is the flat 
space solution, if we add localised matter fields in the 2 extradimensions, it 
will look like a string, the geometry of the two extradimensions being similar 
to a plane in polar coordinates. This is the reason why we will refer to the first 
type of solution as the 'String' (S) branch. The Melvin solution can be found 
by setting k = A = 0, C 7^ in (|2.24p . The string solution is the Minkowski 
space and corresponds toA = K = C = Oin the integral. 

For A > the explicit solution reads: 

flr ■ r\/ ~ \\2/5 r ^ cos (A(r - fp)) 

M = Ci sm (A(r - ro)) ' , C2 — (2.27) 

sm(A(r-ro)) 

where C'l, C2 are parameters again related to C, Lq and = 5 A/8. 

This solution is periodic in the metric functions. In the 4-dimensional ana- 
logue, i.e. Nielsen-Olesen strings in de Sitter space, periodic solutions also 
appear for trivial matter fields [571 EH] ■ 

It is easy to see that the Melvin solutions (|2.26p can be obtained from this 
solution for special choices of the free parameters and specific limits (for A ^ 0) 
of the trigonometric solution in p.27[) . 

The solutions for A < are given by 

M = A± exp(±crr), L = B± exp(±crr), = -J^ (2-28) 

where A± , B± are constants related to Lq in the case where C = and 

M = sinh (A(r - ,o) f^ , L = ^^^^Hr ^Jo)) ^9) 

sinh (A(r — fo)) 

where Ci, C2 are parameters again related to C, Lq and A^ = — 5A/8. 
Inflating branes (k > 0) 

The general explicit form of M is involved; it depends on elliptic functions. 
If A = 0, in the case k = 0, the integration can be done by an elementary 
change of variable, as we have just seen. However, in some particular cases, it 
is possible to integrate (|2.24p explicitly for k ^ 0. 
In the particular case A = 0, C = we find 

Af = W ^(r - fo), L = Lo = LqJ ^ = constant. (2.30) 
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This latter solution corresponds to a cigar-type solution (in the extra dimen- 
sions; 'Cigar-type' refers to the fact that gee is constant) and we find that for 
C 7^ the solutions are also of this type, however they cannot be given in an 
explicit form. 

Consequently, the circumference of a circle in the two extra dimensions 
becomes independent of the bulk radius p. An analogue solution in 7 dimen- 
sions with monopoles residing in the three extra dimensions has been found 
previously in [501 [55] . 

In the case A > 0, k > and C = 0, we find 
A/=/^sin(^yA(,_,„)^, i.^„eos(^yA(,_,„)y (2.31) 

Again, the metric functions are periodic. The periodicity of string-like solu- 
tions seems to be a generic feature - independent of the number of space-time 
dimensions or of the type of brane present in the case of a positive cosmological 
constant. 

Note that by analytic continuation we find 

^ = (/^^'' " ^"^) ' ^ = (/^^'' ''^) ^^-^^^ 

for A < 0, K > and C = 0. 



2.3.4 Boundary conditions 

We will now consider the system (|2.22p . (|2.23p for general /, P. Before solving 
the equations, we need suitable boundary conditions. We require regularity at 
the origin r = which leads to the following boundary conditions: 

/(O) = 0, P(0) = 1, M{0) = 1, M'\r=o = 0, i(0) = 0, L\r^o = 1- (2.33) 

Along with [SOI [HI |5D1 [SS] , we assume the matter fields to approach the vacuum 
configuration far from the string core (i.e. for r 3> 1): 



/(r — >• cx)) = 1, P{r oo) ^0. 



(2.34) 
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2.3.5 Behaviour around the origin 

Close to the origin r = 0, the functions have the fohowing behaviour: 



f{r « 1) 

P(r< 1) 
M{r < 1) 

L(r < 1) 




where po, fo are real constants. 

Note that the behaviour found in [60] is recovered for k = 0. 

For large values of r, the matter functions reach their asymptotic values 
P — > 0, / — 1. As a consequence, the metric functions asymptotically approach 
the special solutions described in section [2.3.31 

2.3.6 Four dimensional effective action 

In this section, we will apply the ideas presented in the first chapter regarding 
the effective lower dimensional theory. To this end, we will consider the vacuum 
solution, assuming that the radial extradimension ranges in the interval [ri, r2], 
ri < r2. Note that we don't exclude the possibility ri = 0, r2 — >■ oo. In this 
context, the length of the extradimension in the radial direction is r2 — 7'i. 

In the vacuum configuration, where / = 1, P = 0, it follows from the last 
two equations in (|2.2ip that L{x) oc M'{x) unless M{x) is constant. This will 
allow to write explicitly the effective action in four dimensions by formally in- 
tegrating the extradimensional dependence of the fields in the Einstein-Hilbert 
action. 

Due to the ansatz we use for the metric, a dimensional reduction of the 
gravity action ()2.2p can be performed. For this purpose notice that we can 
write the 6-dimensional Ricci scalar R according to: 

rW m" M'2 L'M' L" 

where P^^-* is the 4-dimensional Ricci scalar associated with the metric on the 
brane. 

Using the fact that L{r) = cM'{r), c constant, we can integrate the ex- 
tradimensional dependence, along with the remark following equation (|1.19p of 
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chapter [TJ Doing so, we obtain the following 4-dimensional effective action: 



where 



gravity 



167rG 



eff 



- 2A 



(4) 
e// 



(2.37) 



eff 



5 



(2.38) 



and where 



= Mi = 



eff 



2cn 
3Gi 



-M' 



(2.39) 



where Mp is the four dimensional Planck mass. 

Furthermore, we can evaluate the value of c for the vacuum solution, using 
the boundary conditions: £'(0) = cM"(0) = 1. Using the equation of motion 
and the boundary conditions on M,M', we can compute the value of A/"(0), 
leading to c = 

It should be noted that due to the solitonic nature of the matter fields, we 
expect this low energy effective action to be a good approximation of the case 
where the matter fields are not in a vacuum configuration; the four dimensional 
Planck mass should be well approximates by (|2.39p with c = -^-^ . 



2.3.7 Numerical results 

In absence of explicit solutions for non trivial f,P, we have solved the system 
(|2.22p . p.23p numerically and present the solutions in the next few sections. 
Following the investigations in |50( I59| , here we mainly aim at a classification 
of the generic solutions available in the system. The results are organised 
according to the domain of k and A. 



2.3.8 Static branes {k = 0) 

We solve the system of ordinary differential equations subject to the above 
boundary conditions numerically. The system depends on three independent 
coupling constants a, 7, A. We here fix a = 2, corresponding to the self dual 
case in flat space, and we will in the following describe the pattern of solutions 
in the 7-A plane. Results for the 4-dimensional gravitating string 153] make us 
believe that the pattern of solutions for a 7^ 2 is qualitatively similar. 

For vorticity ri > 1, the situation might change (see [53]), however, we don't 
discuss this case here. As will become evident, the presence of two cosmological 
constants leads to a rather complicated pattern of solutions. The interconnec- 
tion of the different type of solution available is illustrated in Fig |2.9l 
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Zero or negative bulk cosmological constant 

The case A = was studied in detail in [55l [SQ. We review the main results 
here to fit it into the overall pattern of solutions. The pattern of solution in 
the 7-A plane is given in Fig l2.1l 




For A = and 7 < 2 two branches of solutions exist, with an asymptotic 
behaviour of the metric functions given by (|2.25p and (|2.26p . Referring to their 
counterparts in a four-dimensional space-time |53| we denote these two families 
of solutions as the 'string' and 'Melvin' branches, respectively. The terminology 
used e.g. in [S3] will be used throughout the rest of the chapter. 

Specifically for a = 2, we have M{r) = 1 and L{r ^ 1) ^ ar + b, a > 0. 
For 7 = 2 the two solutions coincide and L{r ^ 1) = 1. When the parame- 
ter 7 is increased to values larger than 2 the string and Melvin solutions get 
progressively deformed into closed solutions with zeros of the metric functions. 
For 7 > 2 indeed, the string branch smoothly evolves into the so called in- 
verted string branch (again using the terminology of [53]). The inverted string 
solutions are characterized by the fact that the slope of the function L{r) is 
constant and negative, L{r) therefore crosses zero at some finite value of r, say 
r = ris. 

On the other hand, the Melvin branch evolves to the Kasner branch, a 
configuration for which the function M{r) develops a zero at some finite value of 
r = tk while L{r) becomes infinite for r ^ tk- More precisely, for <C r < rn 
these solutions have the behaviour M cx (r^ — r)^/^, L oc {tk — r)^'^/^ . The 
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transition between Melvin solutions (for 7 < 2) and Kasner solutions (for 7 > 2) 
is illustrated in Fig. 12.21 




Figure 2.2: The pattern of the transition from the Melvin solution (7 < 2) to 
the Kasner solution (7 > 2) for A = if given for the metric functions M and 
L. 



For negative A and 7 < 2, the string and Melvin solutions are still present 
and merge into a single solution at some critical value of A. For A < and 
7 > 2, we have the zero r/s of the inverted string solution increasing with 
decreasing A. It reaches infinity for the exponentially decreasing solution, the 
so-called "warped" solution that localizes gravity on the brane [551 [5S] . If the 
cosmological constant is further decreased the solution becomes of Kasner- type. 



Positive bulk cosmological constant (A > 0) 

Up to now static branes have only been discussed A = or A < 0. Here, we also 
consider the case of static branes with A > 0. First, we examine the evolution 
of the string and Melvin solutions for A > 0. For all solutions constructed with 
A > 0, we were able to recover the behaviour (|2.27p asymptotically. 

This evolution is illustrated in Fig 2.3(a) and Fig 2.3(b) for 7 = 1.6 respec- 
tively for the string and Melvin solutions. Here, we show the metric functions 
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M, L for A = and A = ±0.005. For A > we find a solution with the metric 
function AI{r) possessing a zero at some finite r, say r = ri. At the same time 
L{r ri) +00. These solutions tend to the string solutions in the limit 
A 0. Following the convention used in the A = case, we refer to these 
solutions as of 'Kasner' type. 




(a) The profiles of the metric functions A'l, (b) The profiles of the metric functions M, L 
L of the string solution (for A = 0, A = of the Melvin solution (for A = 0, A = —0.005), 
—0.005), respectively of the Kasner solution respectively of the inverted string solution (for 
(for A = 0.005) for 7 = 1.6. A = 0.005) 7 = 1.6. 

Figure 2.3: Profiles of the metric functions L and M for the string, Kasner and 
Melvin solution. 



The second type of solutions that we find has metric functions behaving for 
< r2 < ra as 

L{r2) = , A/'(r2) = , M(r3) = , lim L{r) ^ -00 . (2.40) 

r— i-ra 

These tend to the Melvin solutions in the limit A and we will refer to them 
as of "inverted string"-type. 

For fixed A > and increasing 7 we find that both Kasner and inverted 
string solutions exist for all values of 7. This is demonstrated in Fig l2.4l where 
the values of the parameters Ci, C2, fo (defined in Ea. (|2.27p ) are plotted as 
functions of 7. 
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T 

Figure 2.4: The values of the parameters Ci,C2,ro defined in (j2.26p for the 
inverted string and Kasner-like solutions as functions of 7 for A = 0.003 



2.3.9 Inflating branes k > 

Here, we discuss inflating branes (k > 0). Again, we fix a = 2. The pattern of 
solutions can largely be characterized by the integration constant C appearing 
in 



Zero bulk cosmological constant 

First, we have constructed solutions corresponding to deformations of the string 
solutions residing in the extra dimensions. We present the profiles in Fig. 12.51 
for K = 0.003 and for comparison for k = 0. Obviously, the presence of an 
inflating brane (k > 0) changes the asymptotic behaviour of A/. L drastically 
in comparison to a Minkowski brane. The function M(r) now behaves linearly 
far from the core of the string, while L(r) tends to a constant. The solutions 
approach asymptotically (|2.30p . The space-time is then cigar- like. 
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Figure 2.5: The profiles of the metric functions Af , L for A = and for k = 
and K ~ 0.003; respectively. 



This can be explained as follows: In the case a = 2 the equations are 
self dual, in particular the equation determining the function Af on the string 
branch is MAf" + (3/2)(Af')2 - k/2 = (the combination of the energy mo- 
mentum tensor vanishes identically for q = 2 and for the string like solution). 

The value of C compatible with the boundary condition turns out to be 
C = — k/4, since we are interested is small values of k, the integral (j2.24p 
can be reasonably approximated by (j2.30p . in complete agreement with our 
numerical results. 



The parameters Lq and C appearing in (|2.24p can be determined numeri- 
cally. It turns out that for the cigar-like solutions we always have C < 0, while 
Lq is positive. For a fixed value of k and varying 7, we find that C — > for a 
critical, K-dependent value of 7. 7cr('t)- At the same time Lq ^ for 7 — > jcr- 
This is shown for three different values of k in figure 2.6(a) and 2.6(b) In 
the limit 7 — > 7c,-, the diameter of the cigar tends to zero and the cigar-like 
solutions cease to exist. 
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^=0-001 

K=0-003 

K=0 005 



1.6 1.65 1.7 1.75 



1.85 1.9 1.95 2 



(a) The parameter Lq. 



(b) The parameter C. 



Figure 2.6: The parameters C and Lq for the cigar-hke solution as functions of 
7 for different values of k and A = 0. 



We have also studied solutions corresponding to deformations of the Melvin 
solutions. It turns out that the k = Melvin solution is smoothly deformed 
for K > 0. C is always positive and tends exponentially to zero as function 
of 7. For 7 > 2 we find that the solutions have a zero of the metric function 
L at some 7 and k dependent value of r. The solutions are thus of inverted 
string- type. 

We present the pattern of solutions in the 7-K plane in Fig |2.7l 
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Figure 2.7: The pattern of inflating brane solutions in the j-k plane. 



2.3. EINSTEIN ABELIAN HIGGS MODEL 



31 



Positive/negative bulk cosmological constant 

We have limited our analysis here to 7 < 1.8. For 7 > 1.8, the numerical 
analysis becomes very unreliable, this is why we don't report our results here. 

Fixing 7. e.g. to 7 = 1.8, the inverted string solutions (available for A > 0, 
K = 0) gets smoothly deformed for k > 0. The constant C is positive for all 
solutions. For fixed k and A ^ the Melvin solution is approached. 

The Kasner solutions (for A > 0) get deformed and the functions M, L 
become periodic in r asymptotically. These are deformations of the periodic 
solution ((OT|l with C ^ 0. This is illustrated in FigHjH 

The value of C for these solutions is negative. The function M oscillates 
around a mean value given by Mmv = {-yT'Y^^ a^nd stays strictly positive, the 
solution is therefore regular. The period of the solution depends weakly on k. 
In the limit, A — >■ the periodic solutions tend to the cigar-like solutions. 

We notice that the values of the constants k and A can be chosen in such 
a way that 2k — A becomes arbitrarily small. This seems compatible with the 
domain of parameters leading to periodic solutions. The corresponding four 
dimensional Planck mass determined through (|2.39p can therefore be made 
arbitrarily large, compatible with observations. It is remarkable that the region 
allowing a realistic Planck mass is naturally compact. 




Figure 2.8: The profiles of the metric functions M, L for A = 0.0008 and for 
K = 0.003 and 7 = 1.8. 
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A sketch of the pattern of the solutions is proposed in figure 12.91 The fact 
that we obtain periodic solutions constitutes a natural framework to define 
finite volume brane world models. With view to the discussion on the effective 
4-dimensional action of section !^. 371 we could imagine to put branes with 'large' 
effective gravitational coupling at r ~ rj-, k = 1,2, 3,..., where corresponds 
to the position of the maxima of M{r). Correspondingly, branes with 'small' 
effective gravitational coupling could be put at r = f^, k = 1,2,3, .., where 
are the positions of the minima of M{r). Though there is no 'warping' in the 
sense of Randall and Sundrum [SlIZ]) nevertheless gravity looks stronger on the 
branes positioned at the maxima of M{r) in contrast to gravity on the branes 
positioned at the minima. Note that at the positions of these branes, the metric 
function L(r) vanishes, i.e. L{rk) = L{fk) = 0. Thus the cylindrical geometry 
of the extra dimensional space shrinks to a point and the whole space-time 
geometry becomes effectively 5-dimensional. 




Figure 2.9: The pattern of brane solutions in the K-7-A domain. The origin of 
the coordinate system here corresponds to k = 0, 7 = 1.8, A = 0. M, S, IS, K, 
Per, C denote Melvin, string, inverted string, Kasner, periodic and cigar-type 
solutions, respectively. 



In order to complete this diagram, we summarize the various type of solu- 
tions with the corresponding domain of parameters for A > in the following 
table: 
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Type 


Form of the asymptotic 


Range 


K 


7 


A 


String 


M oc Cste, i cx Lo(r — ro) 


[0,oo] 


> 


< 2 





Melvin 


M oc (r - ro)^, i oc (r - ro)" 


[0,oo] 





< 2 





Cigar 


M oc (r — ro ) , i oc Cste 


[0,oo] 


> 


< 7c < 2 





Inv. String 


M oc Cste, L oc — (r — tq) 


[0,oo] 


> 


> 2 





or 


M 0, L ^ —GO, for r ^ Tc) 




> 




> 


Kasner 


M oc (r - rcY,L oc {r - re)" 


[0,r,] 


> 




> 


Periodic 


M oc cos, L oc sin 


[0,oo] 


> 




> 



2.4 Global 0{n) monopoles 

In this section, we will consider global 0{n) monopole models in n+4 spacetime 
dimensions, n being the number of extradimensions. In this model, we put 
the the 0{n) symmetry of the extradimensions in relation with the internal 
symmetry group 0{n). 

Global monopole solutions occurs in the Goldstone model describing n 
scalar fields in a field theory globally invariant under the 0{n) transformations. 
The symmetry is spontaneously broken by a potential. In the present context, 
the Goldstone model and the corresponding scalar fields are formulated with 
respect to the n extradimensions: 

Sbrane = J (^{d A^)Hd'' 'S>) - ^(<i>t$ - v'f^ d^+^X (2.41) 

where the n scalar fields $ = {(t?- ^ . . . form a fundamental representation 
of the group 0{n). A is the self-coupling of the potential, v the vacuum ex- 
pectation value of the scalar field. We define the following quantities, which 
should not be confused with the rescaled variables of the previous section: 

7 = Ad, /3 = 87rG<j. (2.42) 

The ansatz for the scalar ri-uplet reads 

0' = MC/P , (2.43) 

where the ^' denote the Cartesian coordinates representing the extra dimen- 
sions. 

Correspondingly, the energy momentum tensor has only diagonal compo- 
nents given by |62) : 

T,' = + ^-^^^^ + ^ic^' ~ (2.44) 
T; = -0'^ + ^^^^ + ^(0^-.Y (2.45) 
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T^:=^"+'-^^^ + ^i<t>'-vr (2.46) 

where (f>' = d(j)/dr. The equation corresponding to the Goldstone field reads: 

/' 4 M' 1 

r + {n- 1)(- + - — )</.' --^cj^ = Xcj^icp' - v') (2.47) 

t n — 1 M 

The appropriate boundary conditions read; 

M(0) = 1 , M'{0) = , 1(0) = , l'{0) = 1 (2.48) 

for the metric functions. In the case when the radial variable can be extended 
to /O = oo, the usual boundary conditions for the scalar field function are 

0(0) = , 0(oo) = V, (2.49) 

for regularity and finiteness of the energy. 

Note that the presence of a cosmological constant can lead to a cosmological 
horizon at p = pc- In such case, an appropriate boundary condition for 4'{pc) 
has to be imposed; this will be discussed in due course. 

The expressions of the energy momentum tensor contain terms of the form 
\/P which also appear in the Einstein tensor. If the gravitational constant /? 
is chosen such that 

^/^v = V^r~2 , (2.50) 

the two inhomogeneous terms cancel. This value determines the so-called 'crit- 
ical monopole'. In this case, (and assuming in addition H = Kd = Q) the 
asymptotic Melvin solutions of section 12.2.11 are acceptable as asymptotic so- 
lutions isniEa- 



2.4.1 Sub-critical monopoles 

The case \fPv < y/n — 2 corresponds to the case of sub-critical monopoles |62) . 
The vacuum solutions for which the functions M , I go asymptotically to infinity 
(i.e. corresponding to A„+4 < 0) are such that the term 0^/^ — >■ 0. 

The asymptotic behaviour of the metric function M{p) remains the same 
irrespectively of the presence of a global monopole while the function l[p) must 
be renormalized according to 

l{pU = , ^^^^^ (2.51) 

Here 1^ corresponds to the function of the vacuum solution while Z,„ corresponds 
to the solution in the presence of the monopole. 
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For 7 > the arguments above do not apply because the terms cannot 
be neglected, since in this case / — >■ for r — >■ cx). However, the profiles of the 
metric functions M, I remain very close to those of the vacuum solution. For 
larger values of p, the metric become singular at some finite value of p. The 
singularity is of Melvin type and is of the same nature as in the case of local 
monopoles (see discussion below). 

2.4.2 Mirror symmetric solutions 

The coupled system of equations possesses several symmetries, namely invari- 
ance under translations of the radial variable p and invariance under the re- 
flections p —p and (p — > —(p. These symmetries suggest that solutions which 
are invariant under suitable combinations of the symmetries could exist. 

In the case of vacuum solutions, the solutions corresponding to 7 > possess 
such a symmetry. The most natural combination of the symmetries suggests 
to look in particular for solutions with the following properties 

Kpo - P) = Kpo + P)^ ^Hpo - p) = AI{pq + p), (j){pQ - p) = e(f){po + p), e = ±1 

(2.52) 

where the reflection point po depends on the various coupling constants. The 
existence of solutions presenting one of the above symmetries can be analysed 
by solving the equations supplemented by the following boundary conditions 
at p = Pa- 

I'ifo) = , M'ipo) = , 0(po) = if e = -1 , cj)'{po) = if e = 1 . (2.53) 

These conditions complete the ones given above for p = and allow for a 
numerical study of the equations (no explicit solution was found for </> 7^ 0). 
Our numerical analysis suggests that 

1. solutions corresponding to (j){po — p) = 4>{po + p) do not seem to exist. In 
fact we were able to construct numerically such configurations but they 
do not persist when increasing the accuracy. 

2. Solutions obeying (f>{pQ — p) = —(j){po + p) do indeed exist for peculiar 
values of the coupling constants. 

The existence of these 'odd' solutions can be related to the fact that, in the 
neighbourhood of the symmetric point po, the Goldstone field equation can be 
put into the form 

0" - Xcf,^ + {Xv' - J^)'t> - (2.54) 

where we used the fact that M'{po) = 1'{pq) = 0. 

This simplified version of the Goldstone equation is identical to the kink 
equation provided \/XI{pq)v > 1 which turns out to hold true in the region 
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of parameters that we have explored. Kink-hke solutions can therefore be 
expected. 

Two solutions of this type are presented in figure 12.101 These solutions are 
similar to the so called 'Bag of Gold' solutions discovered in the 4-dimensional 
Einstein- Yang-Mills equations with a positive cosmological constant [63]. Sim- 
ilar phenomena in 4-dimensional space-time were observed in |64| and more 
recently in |65j . However, to our knowledge, it is the first time that compact 
solutions relative to the spatial extra dimensions are constructed. 




10 20 30 40 10 20 30 40 



(a) The profiles of the metric and scalar (b) The profiles of the metric and scalar 

functions for the mirror symmetric solution functions for the mirror symmetric solution 

corresponding to pc = 40. This is a solution corresponding to pc = 40. This is a solution 

on the first branch. on the second branch. 

Figure 2.10: Typical profiles of the two mirror symmetric branches. 



It has to be stressed that this type of solutions exist only for peculiar values 
of the cosmological constants 7, H once the coupling constants A, /3 are fixed. 
Setting A ~ 0.2, /3 = 0.3, the relations between 7. H and pc allowing for 
mirror symmetric solutions are given in figure 12.111 Our results suggest that 
several branches of solutions could exist. We were able to find two non trivially 
different branches, see figure 12.111 The construction of the second branch is 
however involved numerically. 
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Figure 2.11: The relations between the parameters 7, H, po allowing mirror 
symmetric solutions. 

It is clear that these solutions can be extended symmetrically for p e 
[Pc,2pc] and further for p G [2pc,4pc]- This leads naturally to periodic so- 
lutions on [0,4/9c]- The corresponding spacetimes are deformed hyper-spheres 
with angular deficit. 

2.5 Einstein Yang-Mills Higgs 5*0(3) monopole 

It is well know that monopoles exist in some spontaneously broken non-abelian 
gauge theories. The simplest case is the Georgi-Glashow model with gauge 
group SO (3) and a Higgs triplet. In this context, the matter lagrangian reads 

Lrn = -(L'a*'')^^*") - ^FXbF"^^^ - ^(«'°$" - (2.55) 

where A,B = 0,1,. ..,6, a = 1,2,3 and where the usual definitions for the 
covariant derivative and field strengths are used: 

^M$" = dM^'^+ee-'^W'^j'^^, F^jj^ = dMW^~dNW^j+ee'''''WljW^ (2.56) 

Along with [S3] , we use a spherically symmetric ansatz for the gauge and Higgs 
fields, where the 5*0(3) symmetry is imposed in the three-dimensional space 
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defined by the three extradimensions: 

W^; =0, W^r = (l-i/;(p))eay^, a,z,j- 1,2,3, (2.57) 

$'' = </.(p)^, a = l,2,3, (2.58) 
P 

■C^i C^i being the Cartesian coordinates in the extradimensions such that = 
■Ci + ^2 + ?3- The metric ansatz is given by (|2.3|) with n = 3. 

The resulting energy momentum tensor has the following non-vanishing 
components: 

h 1 f — 

,J2\ , '-^ I ±2 in2 



1^ £.272 \ o;2 ^ ^ J ^ A W ^ J 



T^ = -A£)! + ^ + ^((i^_^'2)^«(^2_i) 

2 P eH^^ 2/2 i+ 41'?' 

TS: - ^^-'^^\(-(' -"')■'■ (^■'"» 

111 the absence of a bulk cosmological constant, this model was studied in 
|66j . Here, the equations are solved with the usual boundary conditions for the 
functions w{p),(j){p): 

(t){Q) = 0, w{0) = 1 (2.60) 
for regularity at the origin and 

4>(p oo) = V, w(p oo) = (2-61) 

outside the monopole core, for finiteness of the energy. 

We will analyse the influence of the bulk cosmological constant on the so- 
lutions obtained in |66j . Again, mirror symmetric solutions will be available, 
specific boundary conditions being imposed at some intermediate value p = pc- 
111 the discussion and the figures, we use the following rescaling of the param- 
eters: 

/3 = STre^^Gd, 7 = Ad/e^v^^H/ev H. (2.62) 

In general, the three types of vacuum solutions presented in section 12.21 are 
deformed by the presence of the local monopole on the brane. The deformation 
of the metric fields for 7 = 0, ±0.01 are sketched in figure [2J2l for H — 0.15. 
We will now discuss the different solutions more qualitatively. 

2.5.1 7 = 

In this case, the solutions have been analysed in details in [S3]. Far from the 
monopole core, the metric functions behave linearly: 

Mip) = ^Hp + mo + ^, lip) = ^p + l,+ ^j. (2.63) 
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Figure 2.12: The profiles of the metric functions I, M corresponding to an 
inflating brane with H ~ 0.15 regularized by a monopole are given for 7 = 
0,±0.01. 
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Figure 2.13: The profiles of the metric functions I, M and the ratio /'//, M' /M 
corresponding to a static brane are shown for 7 = —0.01 



These approach the corresponding vacuum solutions of section [2.21 for p ^ 00. 



2.5.2 7<0 

For 7 < our numerical analysis reveals that there are regular solutions which 
obey asymptotically the form 



M{p) — Mpsinh/xp , l{p) — li^smhpp^ p = 




(2.64) 



The profile of such a solution is presented in figure 12.131 for = and 
7 = —0.01. In particular, we can see the deviation of the functions M, I from 
their asymptotic behaviour. 



2.5.3 7>0 



We studied the solution for 7 > and found qualitatively how the vacuum 
trigonometric solutions are deformed by the matter fields. This is illustrated 
in figure 2.14(a) As expected, we observe that the monopole regularizes the 



singular vacuum solutions such that the spacetime become regular on the brane 
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(a) The profiles of the metric functions of an (b) The details of figure [2! 14(a) | in the region 

inflating brane solution with H = 0.07 and of the singular point p Rj 291.7. 
positive cosmological constant 7 = 0.005 are 
shown. 



Figure 2.14: Profile of metric functions for a singular solution. 



p = 0. The functions M, I reach their maximum value when the matter fields 
have already reached their asymptotic values. The mirror symmetry of the 
vacuum solutions is broken, in particular the numerical analysis reveals that 
the two functions M, I do not reach their maximum exactly at the same value of 
p. It is likely that mirror symmetric solution exist in this case as well for tuned 
values of the constants H , 7. Mirror symmetric solutions will be discussed in 
the next section. 



The challenging question is to understand how these solutions (regularized 
at the origin) behave in the neighbourhood of the first period of the associated 
periodic solution of the vacuum equation. 



This is illustrated in figure 2.14(b) which clearly indicates that the solution 
develops a singularity at some p — pc- For H — 0.07,7 = 0.005 we have 
Pc ~ 291.7. A detailed analysis of the numerical solution for p k, pc reveals 
that the singularity is of Melvin type discussed in section 12.2.11 For generic 
values of the inflating parameter and positive bulk cosmological constant, the 
local monopole solutions present a singularity at some finite value of the radial 
variable relative to the extra dimensions. 
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2.5.4 Mirror symmetric solutions 

In the previous sections we reported solutions obtained for generic values of the 
different coupling constants. The counterpart of the vacuum periodic solutions 
turned out to be irregular. However it could be that regular solutions exist for 
specific values of H, 7. Mirror symmetric solutions, like the ones obtained in 
the context of global monopoles are good candidates. It turns out that also in 
the case of local monopoles, it is possible to tune H and 7 in such a way that 
a mirror symmetric solutions exist. In order to construct such solution, we im- 
posed the constraints (p = 0, w = M' = 0, /' = as supplementary boundary 
conditions at some finite value p = pc] then we tried to determine numerically 
if the values of H, 7 can be adjusted for all constraints to be fulfilled. It turns 
out to be possible. 

For fine tuned values of the two cosmological constants, indeed, our nu- 
merical integration of the equations exhibits solutions such that the function 
4i{p) bends backwards and reaches the value (j> = at p = pc after developing 
a plateau where '--^ 1 for < p < pc- At the same time the function w{p) 
and the derivatives of the metric functions / and M, Z', M' approach zero for 
P Po- 

An example of such a solution is given in figure 12.151 It corresponds to 
Pc = 100 which fixes 7 w —0.003, H w 0.187. We constructed a branch of 
solutions for lower values of pc and obtained , e.g. for pc = 50 the values 
7 Ki 0.018, H ss 0.217. The sign of the bulk cosmological constant 7 leading 
to these kind of solutions seems to be negative for pc > 82 and positive for 
Pc < 82. The value H is positive and does not vary significantly with pc. 

The properly speaking mirror symmetric solutions can then be obtained 
by extending the numerical solution discussed above by mirror symmetry for 
p G [0, 2pc] {M,l,w are symmetric under the refiection p — >■ 2pc — p while </) 
is antisymmetric). This is possible by using the discrete symmetries of the 
equations as discussed above. These mirror symmetric solutions can further be 
extended periodically on the interval p € [0, 2ppc] for p integer. The periodicity 
of the function suggest that the geometry of these solutions in the extra 
dimension consists of the surface of a deformed 3-sphere. We can then interpret 
the radial coordinate p as the "colatitude" and the solution looks like a magnet 
(with non-abelian fields) whose positive and negative poles lay respectively on 
the north (p = 0) and south (p = 2pc) poles of the sphere. 

In the region of the equator (p = pc) the Higgs field form a domain wall 
which separates these two poles. This is illustrated by figure 12.161 where the 
Ricci scalar R (see Eq. (|2.8|) ). the full energy density Tq and the contribution 
to the energy of the Yang-Mills fields T^y^^ of the solution of figure [2T5l are 
represented in the north-pole region and in the equator region. In particular, 
we see that only the Yang-Mills energy is zero in the equator region, so that 
the energy content in this region is due to the Goldstone boson. 
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Figure 2.15: Profiles of metric functions for the periodic monopole solutions. 




Figure 2.16: The Ricci scalar R, the energy momentum tensor Tq and the Yang- 
Mills contribution to the energy Tq corresponding to the periodic monopole 
of figure |2.15(a)[ |2.15(b)[ 
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2.6 Gravitating baby Skyrme model 

In this final section, we consider another type of solitonic matter fields, namely 
the baby-Skyrme model. The standard baby-Skyrme model is formulated in 2 
spatial dimensions and admits a solution: the baby-Skyrmion. 

The original Skyrme model was proposed in the 60's [ST] as an effective 
model to describe the pion-nucleon interaction. It is based on a mapping be- 
tween the 3 + 1 dimensional spacetime and the Lie group SU{2), which is known 
to be equivalent to the manifold 6*3. The model admits a fundamental soliton, 
the Skyrmion which is interpreted as a nucleon. The Skyrme model and the 
corresponding classical equations are notoriously difficult. For instance, the de- 
scription of nuclei of high baryonic number is possible but technically difficult 
to construct since it involves several interacting Skyrmions. 

The baby-Skyrme model was proposed in the 90's [47] in order to provide 
a toy model where the theoretical ideas of the grand brother Skyrme model 
could be tested. The baby-Skyrme model is constructed for a 2 + 1 dimensional 
spacetime and involves fields taking values on the manifold 82- Skyrme and 
baby-Skyrme models are non linear sigma models. 

In our framework, the spatial dimensions of the various models we consid- 
ered were interpreted as the extradimensions. It is therefore natural to consider 
a spacetime in 6 dimensions and to interpret the baby Skyrmion as the local- 
ized brane. The 2 extradimensions are the spatial dimensions where the baby 
Skyrme resides, following the philosophy of the various extended brane mod- 
els considered in this chapter. The case of the static baby-Skyrnie branes has 
already been considered in [JSl [SH] ■ Here we will focus on the case where the 
brane are inflating. 

In this section, we used a modified version of the baby-Skyrme model: 

• we consider the static version of the model 

• we replace the usual space by the extradimensional part of the six 
dimensional metric (12.31) with n ~ 2. 



(fx. 



— — {din) ■ (d^n) ^(^^n x djn){d^n x d^n) — Kf)V{n) 



Modified in this way, the relevant action of the baby-Skyrme model reads 



(2.65) 

where kq,K4,K2 are coupling constants, n is a three dimensional vector field 
taking values on the unit sphere (i.e. = 1), 1^ is the potential, which will be 
detailed later, gij is the six dimensional spacetime metric (|2.3p with n = 2. 

Note that following the modification of the model, n depends only on the 
extradimensional coordinates, i.e. r, 9. It follows that the action can be written 
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Sra = I I ^/^iVhh'' 



Y(a^n) • {djH) - h'^h^^'^{dS x dkn){djn x din) 



- KoV{n)] drdOd^x, (2.66) 

where is the determinant of the 4 dimensional part of the metric and h is 
the determinant of the extradimensional part of the metric, given by 

ds\jj dr'^ + L{rfde^. (2.67) 

It should be noted that the four dimensional part of the metric somehow plays 
a spectator role, since the vector field ft doesn't depend on the four dimensional 
coordinates. However, the warp factor M will contribute non trivially. 

The aim of this final section is to investigate the counterparts of the three 
sorts of extradimensional geometries due to the two cosmological constants, 
in a different context. The first three section considered models where Higgs 
fields were present; here it is not the case. As we will see, the geometry of the 
extradimensions will again be opened, closed or flat. 

2.6.1 The ansatz and boundary conditions 

The baby Skyrmion is usually constructed by using the hedgehog ansatz for 
the vector n: 



ft ~ (sin(/(r)) cos(n6'), sin(/(r)) sin(n6'), cos f{r)) , (2.68) 

where n is an integer and /(r) an arbitrary function. We supplement the ansatz 
with the following boundary conditions: 

/(0) = -(m-l)^, /(oo)=^, (2.69) 

where m is also an integer. The condition on /(O) is necessary for the mapping 
to be well defined at the origin and the condition on /(oo) ensures that all 
points at infinity in the extradimensional manifold are mapped to the south 
pole of the sphere. 

Note that this is a parametrisation of 5*2, where /(r) plays the role of an 
colatitude angle and 9 an azimuthal angle; it is clear from the ansatz that 

= 1 for all values of /(r) and 6. However, it is important to stress that 
the presence of the integer m and n have the consequence that n covers many 
times the sphere: m times in the colatitude direction and n times azimuthally. 

The winding number N of the mapping n is a topological invariant, defined 

by 

N ^ — ( n- {dru X den) drdO = - (1 - (-1)'") G Z. (2.70) 
All J 2 
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This is nothing else than the flux of the vector ft thought the sphere. 
We choose the potential term to be [11] 

V{n) = l-n- n'-°°^ = 1 + cos(/(r)), (2.71) 

where ~ linir->oo ft = (0, 0, —1) is the vacuum configuration of the baby 
Skyrme. 

As usual, the boundary conditions for the metric fields are given by 

M(0) = 1, M'(0) = 0, i(0) = 0, L'{0) = 1. (2.72) 

2.6.2 Equations of the baby-Skyrme brane 

First we introduce the following dimensionless coordinates and quantities 

r = . &, l{r) = .[^L{r), (2.73) 
a = 87rG6'«2, P = — 2", 7 = — \ — — (2.74) 

^2 ^2 

and introduce the functions u{r)^v{r) such that 

9 9 
uir) = 1 + sin2 /(r), z;(r) = 1 - sin^ /(r). (2.75) 

Then, the equations of motions are given by: 

+ yTl + + -)"/'- (1 + /") sin / cos / + sin / = 0, 
M'2 L" M'L' M' 7 

6 — ^ + 4 2^ = a(Tr - 0), 

M2 Mi A/2 ^ 

Af" A/f'2 T 

where 

TO ^ -ir-i^sinV-^^rsinV-Ml+cos/), 

= ir-i^sinV+^^rsinV-^l + cos/), 
re - -ir + i^sinV+^^rsinV-M(l+cos/) (2.77) 
are components of the dimensionless energy-momentum tensor. 
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2.6.3 Asymptotic expansion 
Large r limit 

In the hmit where r becoraes very large, the function / tends to the vacuum 
configuration, / = tt. The asymptotic solutions for the functions M and L are 
then the one described in section 12.3.31 

The decay of the functions /, M, L can be computed easily for /? < and 
/3 = by considering the vacuum configuration plus a correction and keeping 
the dominant terms. Doing so, we find 

6M w M-eV^^^'^ + ]\I^e-V^^^\ (2.78) 

for /3 > and where SM, 5L and 5f denote the correction around M{r), L{r) 
and /(r), i.e. 



M{r) = y/-10j/{3al3)smh{y/-a(3/10r) + dM{r), 
L{r) = LoM\r) +6L{r), (2.79) 
/(r) = 7r + 5/(r). 

The parameters /j~ , , A/j~ , M2 , > ^2 arbitrary constants. Of course, 
we are looking for the modes with = /2~ = = such that 6f, SM and 
6L vanish asymptotically and are indeed fluctuations. 

For /? = 0, however, it is possible to solve for the corrections without ap- 
proximations: 

Sf = flx-iK._{x), 
/VfO 

5M = + (2.80) 
SL = LU^, 



where x = r. — h /i, /°, Af°, Af" , Li, L2 arbitrary constants and Kn{x) 



is the modified Bessel function of second kind. 



L„-r 



Note that for large r, the function / decays as tt + /i ^ for a real 

constant /i. 

Expansion around the origin 

Close to the origin, the expansion of the functions /, L, M depends on m and 
n. It is possible to give the leading terms in a generic form, but the case n = 
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has to be considered separately: 

/(r) = -(?7i- l)7r + /(I V + C'(^)^ (2.81) 

a(-5(/W)^-4(/W)^-2(-l)-M + 2/? + 2M)-87 r3 , , ,4 
L{r) = r + ^ 37+^ W ' 

In the case where ti > 0, we can write the first terms in a general form: 
/(r) = -(m-l)7r + /(")r" + 0(r)"+\ 

L{r) = + (2.82) 

Mir) ^ l + ^l^i^i^±pi^^ + 0(.)3, 
where are real constants. 

2.6.4 Numerical solutions 

We solved the system of ordinary differential equations numerically with the 
solver Colsys |69j for many values of the parameters. Due to the large number 
of parameters, we decided to adopt the following approach: we keep fixed the 
value of the gravitational coupling a, the bulk cosmological constant /3 and 
m,n. Then we vary the Hubble factor 7 for different values of the Skyrme 
coupling ^; the way we treat the problem allows to have a direct visualisation 
of the influence of 7 on the solutions. 

In the analysis, we focused on the energy density of the Skyrmion E and 
its mean square radius MSR |70j : 

/•oo /'OO 

E = 2tt T^L{r)dr , MSR^ / f {r)siv? f{r)dr. (2.83) 

JQ Jo 

The mean square radius allows to characterise the extension of the brane in the 
transverse direction; the more MSR is small, the more the brane is localized. 

Our results are presented in figure [^T71 and [2 . 1 81 for a = 0.1, /3 = 0.1 and 
a — 0.05, P = 0.1 respectively and m = n = 1 for both cases. Note that 
these are preliminary results, a more detailed analysis, including stability is in 
preparation |11| . 

For fixed a, /3, /i, the mean square radius decreases when the value of 
7 increases, while the energy decreases for small values of 7 and increases 
for larger values. In some intermediate values of 7, the energy of the baby- 
Skyrmion passes thought a minimum. The minimum occurs at smaller values 
of 7 when the value of yu is smaller. Note also that increasing values of /i leads 
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Figure 2.18: The value of the energy and of the mean square radius as a function 
of 7 for various values of /i and for a = 0.05, /3 = 0.1. 



50 



CHAPTER 2. EXTENDED BRANES 



to decreasing values of the mean square radius and decreasing values of the 
energy. 

We present typical profiles of the solutions for rn = 2,3 in figure I^TTOl and 
12.201 respectively for non vanishing values of the parameters. These figures 
show that there are three possible geometries depending on the sign of the 
cosmological constant, as in section [2.3.31 namely opened (/3 > 0), flat {(3 = 0) 
and closed (/3 < 0). Again, all three geometries have an angular deficit. 




Figure 2.19: The typical profiles of the solutions for different signs of f3 for 
m = 2. The function / tends quickly to its asymptotic value, so it is not 
possible to distinguish the different profiles for the function /. 




r 



m=3, n = l 

Figure 2.20: The typical profiles of the solutions for different signs of (3 for 
m = 3. Here again, it is not possible to distinguish the different profiles for the 
function /. 
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2.7 Concluding remarks 

In this chapter, we investigated the effect of both a bulk cosmological constant 
and a brane cosmological constant on essentially four classes of warped topolog- 
ical braneworld models: the gauged Einstein-abelian-Higgs model, the gauged 
and global Einstein-non abelian-Higgs models and the Einstein-baby-Skyrme 
model. 

We reviewed solutions available in a wide region of parameters, giving a 
classification for the solutions to the 6-dimensional abelian Higgs brane model, 
and we presented new solutions where the branes are inflating, completing the 
classification. We also considered models with more than 6 dimensions, since 
the case c? = 6 is particular. 

The four dimensional cosmological constant have been chosen to be posi- 
tive in accordance with the current observational status and was modeled by 
infiating branes. Note that the four-dimensional slices of the generic metric 
(j2.3p can be replaced by any Einstein four dimensional spacetime (i.e. a space- 
time satisfying Rat = Cgab, C £ M) with positive curvature, since only the 4 
dimensional scalar curvature appears in the equations of motions. 

It should be stressed that the inclusion of such a cosmological constant 
has heavy consequences on the solutions available in the models. It seems 
to be a generic feature that the geometry of the extradimensional spacetime 
with positive cosmological constant on the brane is either opened, closed or 
flat, with deficit angles. This is comforted by the vacuum solutions presented 
in section 12.21 which are precisely of opened, closed or flat extradimensional 
geometry. Far from the brane core, the matter fields should asymptote their 
vacuum value, so the metric functions should approach the vacuum solutions 
asymptotically. 

This idea has been checked on different models, three of them (Einstein- 
Maxwell-Higgs, Einstein non abelian Higgs, Einstein Yang-Mills Higgs) con- 
tained a Higgs field, the fourth is somehow different in its nature (baby-Skyrme 
model). In every cases, we found the three sorts of extradimensional geometries. 

Note also that the effect of considering branes with an extension in the 
extradimensions somehow regularizes the metric functions, as it was already 
foreseen in the previous chapter. 
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Chapter 3 

Trapping realistic matter 
fields: a toy model 

In this chapter, we consider a model describing an interaction between a boson 
and a fermion field in 1 + 1 dimensions. The procedure will actually consist 
in considering solutions available in the bosonic sector of the theory and then 
coupling the fermionic fields to these background solutions. The bosonic sector 
is chosen in such a way that it admits localized 'bell-shaped' classical solutions. 

This kind of models can be seen as toy models where the scalar fields mimics 
a brane (and/or an antibrane) on which we investigate bound (localized) states 
of fermions. In the first section, we will discuss the fermionic modes while in the 
second section, we will study the stability of three particular bell shaped scalar 
field solutions and show that the stability equation reduces to the Poschl- Teller 
equation (see appendix [C]) , which admits analytic solutions. These two points 
constitute our original results, published in |16| . 

3.1 Fermionic modes on kink and kink- anti-kink 

systems 

Up to now, we considered brane models described by gravitating topological 
defects such as cosmic strings or monopoles. However, realistic matter is of 
fermionic nature. The goal of this section is to consider a toy model of a 
topological defect seen as a brane and to couple fermionic fields to it. The 
simplest way to do so is to consider something which resemble a domain wall 
or a system of wall/anti-wall (see |15| ) in 1 + 1 dimensions. Among other, 
one possibility is to consider the Sine-Gordon model [TB]. It is well known 
that the Sine-Gordon model admits kink and anti-kink solutions, as well as 
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superposition of these solutions. The solution of interest for our purpose is the 
kink-anti-kink system. 

Previous work has been done in this direction, considering a 0** model, 
where kink or anti-kink are available [7T]. The author of [7T] studied coupling 
of fermions to a single kink and to a superposition of the kink and anti-kink 
solution. However, the superposition is valid as long as the kink and anti-kink 
are distant enough. We will briefly review their result before turning to the 
Sine-Gordon model. But first, we present the formalism used throughout this 
section. 

3.1.1 Fermion - Boson coupling in 1 + 1 dimensions 

The model we consider describes a coupling between a bosonic field to a 
fermionic field and an effective mass for the fermion generated by a Yukawa 
interaction of the fermion with the scalar field 0. Using the notations of Chu- 
Vachaspati [71], the lagrangian is given by (with the metric signature (+, — )) 



where ■(/; is a two component spinor and where the two dimensional Dirac ma- 
trices can be chosen in terms of the Pauli matrices 7° = (73, 7^ = iai, g is 
the Yukawa coupling constant and the shift constant C has to be chosen ap- 
propriately according to the vacuum expectation value of the scalar field. V 
is a given potential, defining the model considered for the bosonic part of the 
theory. 

Let (/ic be a classical solution of the bosonic equation derived from (|3.ip . i.e. 
a solution of 



In order to solve the Dirac equation in the background of a classical solution 
(f>c available in the bosonic sector of the model, it is convenient to parameterize 
[7T] the Dirac spinor according to 



where I3± are functions of x. 

Then the system of first order Dirac equations is transformed into two 
decoupled second order equations for /?_ and : 



(3.1) 




(3.2) 




(3.3) 



i'dl + y±(0c))/3± = E^p± , F±(0e) = g^(j>l T gdAc- 



(3.4) 
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3.1.2 Coupling to a 0^ kink 

In the case where the potential V is given by 

T/(0) = i(</)2-l)', (3.5) 
the bosonic part of the lagrangian p.ip admits static locahsed solution given 

by 

(j){x) = zbtanhx. (3.6) 

The solution corresponding to the plus sign is called a kink while the solution 
with the minus sign is an anti-kink. The solutions, kink or anti-kink, extrapo- 
late between the two vacua of the potential (j)± = ±1, between x — ?• — oo and 

X — > +0O. 

The Dirac equation in the background of one kink corresponding to this 
model was studied in detail in |71| . It was found in particular that the equations 
(|3.4p with the solution p.6p are Poschl- Teller equations respectively with N = g 
and N = g — 1 and ^ E'^ — g^ . For any value of g such that n — 1 < g < n 
there exist n normalisable solutions with eigen energy 

Ej = V.]i^9-:i) , j = 0, 1, . . . ,n - 1 (3.7) 

At each integer value of a supplementary bound state emerges from the 
continuum E = g and exist for g > n. 

Next, the Dirac equation was studied in the background of a kink-anti-kink 
configuration say 0^7^ approximated by a superposition of a Kink centred at 
X ~ —L and an anti-Kink centred sX x ~ L with L ^ 1: 

4>xTi{^) = tanh(.T + L) — tanh(a; — L) — I. (3.8) 

Note that this is indeed a good approximation of a kink/anti-kink solution, 
due to the localized properties of the solution. However, if the separation L 
between the kink and anti-kink becomes too small, then the approximation 
does not hold any more. 

Configurations like (|3.8p are "bell-shaped", in the sense that the graph of 
the solutions looks like a bell. 

The main result is that the spectrum of the Dirac equation in the back- 
ground of (jij^j^ deviates only a little from the spectrum available in the back- 
ground of a single Kink (or an anti-Kink) for L ^ 1. The analysis of the 
spectrum in the region i ~ 1 is unreliable since the linear superposition p.Sp 
is an approximate solution of the classical equation only for large values of L. 
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3.1.3 Coupling to the Sine Gordon kink and kink-anti- 
kink 

Since the considerations of the previous section cannot be extended in the 
region where the kink and anti-kink are too close, we reconsidered the spectral 
problem in a different theory, namely the Sine-Gordon model where explicit 
solution are known for the kink, anti-kink as well as for a kink-anti-kink. 
The potential in the Sine-Gordon model is given by 



1 



(3.9) 



and the kink and anti-kink solutions are given by 

(j) = ±4Atan(e^), 



(3.10) 



the -|- (resp. — ) sign describing the kink (resp. the anti-kink). 

We first solved equation p.4p in the background of a single kink (|3.10p . 
On figure [331 a- few fermionic bound states are represented as function of the 
parameter g (for clarity, we omitted the modes n = 4, . . . 9 on the graphic). 
The figure present exactly the same pattern than in the case of the 0^ kink. In 
particular, new bound states emerge regularly from the continuum at critical 
values of g, say g = gn- This is in the case of the single kink solution. 




Figure 3.1: The fermionic bound state corresponding to the single Sine-Gordon 
kink as function of the constant g. The level n = 0, 1,2,3 and 10 are repre- 
sented. 



Because these qualitative properties of the solutions are similar to the one 
of the (f)'^ kink we expect the features discovered in [71] to hold in the case of 
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well separated sine-Gordon kink and anti-kink. In the sine-Gordon model, we 
can take advantage of the fact that an exact form of the kink-anti-kink solution 
is available: 

sinhi ut/B) „ / 

where w is a constant related to the relative velocity of the two lumps. We used 
this solution to study the Dirac equation in the background of a moving wall- 
anti-wall system (approaching or spreading each other, according to the sign 
of t) and were able to study the spectrum of the bound states in the domain of 
the parameter t become small, i.e. when the kink and the anti-kink interact. 

Solving (|3.4p in the background p. lip , we implicitly assume that the mo- 
tion of the kink and of the anti-kink is treated adiabatically. The results are 
summarized on figure \S7I\ where the fermionic eigenvalue E for the ground state 
n ~ and the first few excited states are represented as function of t. 




Figure 3.2: The fermionic bound state corresponding to the sine Gordon kink- 
anti-kink solution as function of time for g = 1 and u ~ 0.5. The level n = 
0, 1, 2, 3, 4 are represented. 

Our numerical integration of the spectral equations p.4p for several values 
of t reveals that only the ground state subsists in the t limit where it 
enters in the continuum (this could be expected since the classical background 
solution vanishes in this limit) and that the excited states join the continuum 
at finite values of t, depending on the level of the excitation. 

It should be pointed out that the fermionic eigenvalue of the fundamental 
solution (line n = 0) approaches E = for f — > oo although staying positive. 
It decays exponentially with time. 
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The fact that the excited energy levels become pairwise degenerate in the 
large t-limit can be understood from the form of the potential. Indeed, for 
t ^ 1, the potential possesses two well separate valleys centred about the 
points X ^ ztut. In the neighbourhood of a; = ztut, the form of the potential is 
given by V± ~ g^cfP'^ g(j)'. Far away from the two regions of the x-line situated 
around x = ±ut (i.e in the region of the origin and in the asymptotic regions), 
the potential is exponentially small. In fact the effective potential under inves- 
tigation results from a superposition of two (suitably shifted) potentials which 
are supersymmetric partners of each other; accordingly they have the same 
spectrum apart from the ground state of K). . It is one of the striking property 
of supersymmetric quantum mechanics (see appendix |D] or j72| for a review) 
that if V±{x) are supersymmetric partner potentials, the fc*'' energy level of V- 
coincides with the {k + 1)*'' energy level of V+. 

As a consequence, in our case, two eigenvectors exists with roughly the 
same eigenvalue when the wall and the anti-wall are well separated, like e.g. 
in figure l373l While t decreases, the two valleys of the potential have tendency 
to merge in the region of the origin, see figure 13.41 and the degeneracy of the 
two energy levels corresponding to each other by the supersymmetry is lifted, 
as shown on figure 13.21 It is tempting to say that the supersymmetry of the 
spectrum occurring for 3> 1 is broken in the — ^ limit. 




Figure 3.3: The profile of the potential and of the first three fermionic bound 
states corresponding to the single sine Gordon kink-anti-kink solution for t = 4, 
(7=1 and u = 0.5 are represented. 
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Eigen states (shifted) for t=1 ,g = 1 
Sine Gordon kink-antikink 

\ I 

\ 




Figure 3.4: The profile of the potential and of the first three fermionic bound 
states corresponding to the single sine Gordon kink-anti-kink solution for t = 1, 
g ~ I and u ~ 0.5 are represented. 



Said in other words, it turns out that when the wall and anti-wall approach 
each other, the system cannot support fermion bound states and their spectrum 
merge into the continuum of the Dirac equation. Seen oppositely, fermion 
bound states (ground state and a number of excited modes, the number of 
them depending on the coupling constant g) can emerge from the continuum 
when the wall and anti-wall separate from each other. 



3.2 Bell-shaped lumps: stability 

In the previous section, several lumps which have the shape of a bell were 
used, mimicking a kink-anti-kink non topological solution. The stability of the 
underlying fundamental solution is well known and can even be formulated 
in terms of a Poschl- Teller potential (see appendix [C|) where some relevant 
stationary solutions can be computed algebraically. However, a more general 
question can be addressed: which potential V{cf>) admit kink like configuration 
as static solution and for which of them is it possible to formulate the stability 
equation as an algebraic equation? We won't answer fully this question, because 
classifying all the possibilities can be a bit ambitious, but we will focus on three 
models, namely the (j)'^, inverted 0^ and 4>'^ models. New lumps solutions to 
these models were presented recently in |73) . 

These models are interesting by themselves but once again, they can be 
seen as toy brane models. The models considered are 1+1 dimensional scalar 
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field theory of the form 

C^^d^^df^cb-Vi^) (3.12) 

with suitable potential V{(j)). 

For static fields, the classical equation has the form d'^cj) = The linear 
stability of a solution, say 4>c{x)^ of this equation, can be studied by diagonal- 
izing the quadratic part of the lagrangian in a perturbation around (j)c- 

Concretely, one has to compute the eigenvalues allowing for normalisable 
eigenfunctions of the spectral equation 

)7? = uj^T] , (f,{t, x) = 0c(.t) + e'"'7j{x) (3.13) 

Negative eigenvectors o;^ < 0, reveal the existence of instabilities of the classical 
solution (pc- 

3.2.1 Case 1: (j)"^ potential 

The case of a quartic potential is well known but we mention it for completeness 
and for the purpose of comparison of the eigenmodes with the cases mentioned 
earlier. 

The potential and the corresponding kink-solution read 

V{(f>) = i(<^2 _ 1)2^ ^ ±tanh(x), (3.14) 

and the stability equation corresponds to the Poschl- Teller equation (see ap- 
pendix [C]) with N — 2 and uj'^ = ujp + 4. The eigenvalues are then = 0, 3, 4. 
It is well know that the kink has no negative mode and is stable. The 
corresponding eigenfunctions are 

1 d0c sinh(x) ■ t,2/ x 1 i r^ 

''0 " TT^ = ^3 TITT' ^4 = tanh a; - -. 3.15) 

cosh(a:)^ dx cosh (a;) 3 

The zero mode corresponds to infinitesimal translations of the classical solution 
(jjc in the space variable. 

3.2.2 Case 2: Inverted 0"^ potential 

This case corresponds to a quartic potential with inverted sign. The potential 
and the corresponding lump are given by 

y(0) = h'ii ~ 4>^), u^) = ±^^. (3.16) 

2 cosh a; 
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In this case also, the stabihty equation corresponds to the Poschl- Teller equa- 
tion with N = 2 but with a different shift of the eigenvalue, for instance 
oj^ ^uj^ + 1; leading to = -3, 0, 1. 

As a consequence the lump is characterized by one negative mode and one 
zero mode. Up to a normalisation they are given by 

1 , sinh(x) d0c - + i 2/ n 1 lv^ 

as usual, the zero mode is associated with the translation invariance of the 
underlying field theory. 



3.2.3 Case 3: Cubic (p^ potential 

The last case investigated corresponds to a potential of third power in the field 
(j}. The potential and the corresponding lump are given by 

y(</)) = 20^(1 - </)), M^) = ± A, (3.18) 

cosh (x) 

The corresponding normal mode equation corresponds to the Poschl- Teller 
equation with iV = 3 with an appropriate shift of the spectrum; cj^ = + 4; 
leading to = -5, 0, 3, 4. 

The lump here is then characterized by one negative mode, one zero mode 
and two positive modes. Up to a normalisation they are given by 

1 sinh(x) 1 dcjjc 

" cosh3(x)' """^ cosh3(x) ="2^' 

1 2 1 2 

7^3 = 2 (tanh (a;) — -), 774 = tanh(x)(tanh (x) 

cosh [X) 5 

In [75] , deformations of the three potentials given above are studied as well 
(e.g. is cases 1 and 3, deformations breaking the reflection symmetry (f> — ?> —</>)• 
The corresponding stability equation is not of the Poschl- Teller type but their 
spectrum could be studied perturbatively from the ones obtained above. It 
should be interesting to see in particular how the unstable mode (or the zero 
mode in the case 1) evolve in terms of the coupling constant parameterizing 
the deformation. 



(3.19) 



5' 



3.3 Concluding remarks 

In this chapter, we studied fermionic localization on kink and kink anti-kink so- 
lutions available in the Sine-Gordon model in 1-t-l dimensions. This kink/anti- 
kink system can be viewed as a toy model for a brane - anti-braneworld solution. 
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Note however that the solution is time dependent, but we treated t as an ex- 
ternal parameter, fixing the distance between the kinks. As mentioned in the 
text, this consists in treating the system adiabatically. The interesting result is 
that the fermionic field is indeed localized on the kinks. Bound states emanate 
from the continuum as the coupling between the kink or kink-anti-kink and the 
fermionic field is increased. 

In the case of the kink/ ant i- kink system, fermionic bound states emanat- 
ing from the continuum appear when the kink and anti-kink start spreading 
apart. Moreover, when the kink and anti-kink are well separated, consecutive 
fermionic levels degenerate, i.e. level 1 and 2, level 3 and 4, etc. 

We also considered the stability of three types of potential, namely 
inverted 0^ and cubic (p^ potential. The stability equation for these three 
models are of the form of the Poschl- Teller equation. 

As a result, lump solutions with the 0^ potential are stable, but solutions 
with the two other potential have a negative mode and are then unstable. This 
seems reasonable since the first potential is positive while the two last have a 
local minima and lim0_j.oo = — oo. Intuitively, such potentials allow a decay of 
localized solutions to the unbounded minima via quantum tunnelling. 
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Higher dimensional black 
holes and strings 

Black holes solutions have always been of considerable interest in general rel- 
ativity. Such objects are solutions to the Einstein equations describing space- 
times having the property that there exists a spatial region whose causal future 
is trapped [55] • The boundary of this spatial region is called the black hole 
horizon. Black holes have been intensively studied in four dimensions. In par- 
ticular, the only possible horizon topology in four dimensions is ^2 - the surface 
of a 2-sphere - and stationary (charged) black hole solutions are classified in a 
three parameters family, namely the Kerr-Neumann family |74| . These three 
parameters are the charge, mass and angular momentum of the black hole. 
Black hole solutions are of considerable interest from an observational point 
of view, since they can result from a collapse of massive stars, but they are 
also fundamental objects from a theoretical point of view. They are actually 
theoretical labs for general relativity. Black holes of the Kerr-Neumann family 
are characterised by a curvature singularity inside the horizon; close to the 
singularity, the curvature blows up and the gravitational interaction becomes 
very strong. 

On the other hand, higher dimensional gravity has been considered long 
time ago, as discussed in chapter [U and have regained interest the last three 
decades, mainly motivated by String theory. 

It seems then natural to investigate the higher dimensional counterparts of 
the black hole solutions available in four dimensional gravity. However, in more 
than four dimensions, there are several drastic changes with respect to what 
is known in four dimensions. First, in higher dimensions, say in d dimensions, 
there are different possible horizon topology, such as Sd-2 (spherical), Sd-3 x 
(cylindrical or toroidal), Sd-2-p x (black branes), etc. There is also the 
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possibility of having several asymptotical spacetime topology, such as x 
Sd-A, Md,--- while in four dimensions, the asymptotic spacetime is A^4, Md 
denoting the d-dimensional Minkowski spacetime. In four dimensions, a black 
hole is uniquely specified once the mass, angular momentum and charge is 
given [75l [76l [77l [78] . It is not the case any more in higher dimensions (see for 
instance the case of the black ring and of the Myers Perry solution, [79l I80|). 

Note that the three parameter family of solutions describing four dimen- 
sional black holes has a counterpart in the presence of a cosmological constant, 
namely the {A)dS- Kerr-Neumann solutions, where the asymptotical spacetime 
is the {A)dS4 spacetime. 

In this chapter, we will first focus on higher dimensional black objects pre- 
senting a spherical horizon topology (5^-2)- We will present such known solu- 
tions in pure Einstein gravity (vacuum) and for the Einstein-Maxwell equations 
(i.e. in the presence of an electromagnetic field). Furthermore, we will present 
the construction of new solutions, namely the de Sitter charged, rotating black 
hole [Him. 

The second part of this chapter is devoted to black strings, black objects 
presenting a 5i x 5^-3 horizon topology, in asymptotically locally flat space- 
time, which have been intensively studied in the literature (stability |19j . non 
uniform black strings [811 182], etc.) and with a positive cosmological constant. 
Black strings with a negative cosmological constant have been studied in pp] 
and will be discussed in more details in chapter [SJ black strings with a positive 
cosmological constant have not been discussed in the literature yet. 

In asymptotically locally flat spacetime, black strings and more generally 
black branes are known to suffer from a long wavelength dynamical instability, 
the so-called Gregory-Laflamme instability [12] . We will briefly review this in- 
stability in the end of the chapter. 

To sum up, this chapter is organised as follows: in section 14.21 we give a 
brief overview of the analytically known higher dimensional black hole solutions 
with a spherical horizon topology and their properties. In section H751 we present 
new black hole solutions, namely the charged rotating de Sitter black holes and 
their properties. Then we review the black string solution and its instabilities 
in section 14.41 finally we discuss the effect of a positive cosmological constant 
on the black string solution in section 14.51 In particular, we give evidences 
that such black strings do not approach an asymptotically locally de Sitter 
spacetime but indeed lead to a singular spacetime. We briefly conclude this 
chapter in section WM 
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4.1 The general framework 

The various solutions presented here are all extrema of the Einstein-Hilbert 
action with boundary term (see appendix minimally coupled matter field 
and a cosmological constant: 

1 r . . . 1 



167rG 



/ ^iR-2A + L,n)d'^x + -^ [ V^Kd'^~^x, (4.1) 

Jm JdM 



where G is the d-dimensional Newton constant, is a d-dimensional manifold, 
dA4 the boundary manifold, g the determinant of the metric, R the scalar 
curvature, A the cosmological constant, h the determinant of the induced metric 
on the boundary manifold, K the trace of the extrinsic curvature Kab (see 
appendix |A| of dAi and Lm is the matter lagrangian. depending on the model 
we will be interested in. 

We will not discuss in detail all the thermodynamical properties of the 
various known solutions apart from the simplest case, namely the Tangherlini 
black hole, but we will present the various solutions instead. The interested 
reader can find more details on the thermodynamical properties in the literature 
[Ml [m [111 [861 [87], between many others (see also for a review). 



4.2 Known higher dimensional black holes 
4.2.1 Tangherlini black hole 

The Tangherlini solution is the generalization of the four dimensional Schwarz- 
schild solution to any number of dimensions. It extremizes the action (|4.ip 
with A = 0, Lm = and is given by 



ds' = -f{r)dt' + ^ + r'dnU, f{r) = 1 - (!^)'"' . (4.2) 



f[r) V r 

The event horizon is located at r = rh. The solution (|4.2p generalizes 
the spatial spherical symmetry of the Schwarzschild black hole and is time 
translationally invariant (it is a static solution) . The vector dt is then a Killing 
vector and the associated conserved quantity (the mass, M) can be computed 
using the standard Komar integral (see Appendix \X\ : 

IdttG 

where Vd-2 is the surface of the unit {d — 2)-sphere. 

The entropy S and temperature Th are computed in the standard way; 
using the Hawking-Bekenstein formula |31l I32| and by demanding regularity in 
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the euclidean section, i.e. by avoiding conical singularity in the near horizon 
Euclideanised geometry of the black hole: 

S^±V,^A-^T„ = '^. (4.4) 

Note that S cx 1/T^~^; this imphes that the specific heat Cp = Th^^ 
of the Tangherlini black hole is negative. This is the typical signature of a 
thermodynamical instability. 

The solution obeys the first law of black hole thermodynamics 

dM = TndS, (4.5) 

and the thermodynamical quantities are related to each other by the Smarr 
relation, 

{d~3)M ^ {d~2)T„S. (4.6) 

These two relations can be checked explicitly for the particular case of the 
Tangherlini solution but it is possible to prove such relation more formally (see 
for example \2G\). 

The Tangherlini solution generalizes to asymptotically {A)dS spacetimes 
with 

= (4-7) 

where £ = ^'^~2|A| {A)dS curvature radius and e is the sign of the 

cosmological constant. 



d-3 



The mass is still given by M ~ '"'^ '^^lelrG^'' — ' ^^^^ horizon radius now 
depends on the cosmological constant. As a consequence, the entropy and 
temperature depend on the cosmological constant as well. The presence of the 
cosmological constant has serious implication on the thermodynamical proper- 
ties of the black hole, for instance, in asymptotically AdS spacetime, the black 
hole develops two phases: small black holes {r+/£ << 1) and big black holes 
r+/i 1, r+ being the largest root of f{r). Small black holes are thermody- 
namically unstable while big black holes are thermodynamically stable; this is 
purely an effect of the AdS curvature [88] . 



4.2.2 Myers-Perry black hole 

The Myers-Perry black hole is the higher dimensional generalization of the 
Kerr solution. However, in more than four dimensions, there are multiple 
independent planes of rotation. It follows that the Myers-Perry solution is 
characterized by multiple angular momenta, associated with the rotation in 
these various planes. 
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In d dimensions, the number of independent planes TV is given by [(d— 1)/2J , 
\_X\ denoting the floor integer value of X. For simplicity, we will consider only 
odd number of dimensions and equal angular momenta. More general case can 
be found in [75] . 

In order to describe rotating black holes, we consider the following ansatz: 



-b{r)dt 



dr^ 



JV-l /i-l 



fir) 

N /fc-1 



i=l \j=0 



def 



h{r) J2 n ^o^' ^i^' ~ w{r)dty 

k=i \i=a / 

Kk=l \i=0 / 
N /fc-1 \ 



(4.8) 



.fe=l \i=0 



In the above formula 9^ := and 7r/2 are assumed; the non trivial angles 

have 9, g [0, tt/2] for i = 1, . . . , A^- 1, while ipu S [0, 27r] for fc = 1, . . . , iV. Note 
that the ansatz (|4.8p imposes p{r) = g{r) ~ h{r) for consistency. 

The Myers-Perry solution |89) with equal angular momenta is then given 
by ([48)) with 



f{r) 



w{r) 



2Ma 



h{r] 



> .9(0 



6(r) = 



r-V(r) 
Mr) ' 



(4.9) 



where M and a are two constants related to the solution's mass and angular 
momentum. 

The symmetries of (|4.8p are given by the following (commuting) Killing 
vectors: dt,dip-, i € [1, [(d — 1)/2J] and once again, conserved quantities can 
be computed using standard Komar integrals leading to the Mass M and the 
angular momenta Ji associated to the N independent rotation planes (here, we 
assume Ji = J, i = 1, . . . ^ N: 



Vd-2d-2^^^ ^ Vd-2d-l^^ 
M = ^^—^ — M, J = — — - Ma, 



(4.10) 



47rG 2 ' 47rG 2 

where Vd^2 is the surface of the unit {d — 2)-sphere. The entropy and Hawking 
temperature are computed in the standard way, leading to 

d+2 



H 



An 



, S 



Vd-2r[ 



(4.11) 
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where r+ is the event horizon, defined as the largest root of /(r). The first 
law becomes dM = TndS + Q^^dJ^, where fi^^ is the angular velocity at the 

ft (i) f) • 

horizon defined by the requirement that k = + f2]j is null on the event 
horizon. In the case under interest, the angular velocity of the horizon reduces 
to ri^'' = w(r^). In the following, we will omit the superscript (i) and the 
subscript H since all angular velocities at the horizon are equal and are always 
to be understood as evaluated at the event horizon, except if stated explicitly. 

The counterpart of the Myers- Perry solution in the presence of a cosmolog- 
ical constant has also an analytic solution |86) . 

2ME 2Ma2 2Ma^ 



2Ma 2 r^jr) 

= ^^rsj^y 9{r)=r, &W = (4.12) 



where M and a are still related to the solution's mass and angular momentum: 

while S = 1 + ea^/i'^. The parameter e is again the sign of the cosmological 
constant. 



4.2.3 Charged higher dimensional black hole 

The generalization of the four dimensional Reissner- Nordstrom black hole with 
(or without) cosmological constant is a solution of the equations of motions 
resulting from (|4.ip with a minimally coupled Maxwell field |90| : 

Lm. = -^FabF^'', Fab = OaAb - ObAa, (4.14) 

where A, B run from to d — 1 and Aa is the Maxwell field. 

The equations resulting from the variation of (|4.ip with respect to the grav- 
itational field and the Maxwell field are then the Einstein-Maxwell equations, 
eventually supplemented by a cosmological constant. 

The d-dimensional {A)dS charged black hole solution is given by 

ds' = -fir)dt' + ^+r'dnj_,, 

^^'^ = '^£^ + ^'H^+ 2(rf-2)(d-3)r^('^-3) ' (4-1^) 



AAdx^ = V{r)dt, V{r) = 



Q 



{d~3)r'i-^'^ 
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where M and Q are the mass and charge of the solution, £ 



(d-l)(d-2) . 
2A 1 



IS 



the (A)^^ radius and e = |A|/A if A 7^ and e = if the cosmological constant 
vanishes. 

4.3 {A)dS charged rotating higher dimensional 
black hole 

As we have seen in the previous sections, known solution are either rotating, 
either charged (with or without a cosmological constant). The missing point is 
the charged and rotating solution. Unfortunately it seems that finding analytic 
solutions is hopeless, so we have to rely on numerical computations. 

Here we consider the dS charged rotating black hole in odd dimensions. 
The charged rotating solution have been constructed in asymptotically flat 
spacetime and asymptotically AdS spacetime [911 US- The difficulty for the 
construction of the dS charged rotating black hole is essentially due to the pres- 
ence of a cosmological horizon. The number of independent angular momenta 
is given by iV = {d— l)/2. We will use the ansatz (|4.8p for the metric field, the 
matter content will be described by the minimally coupled Maxwell action, as 
in the previous section. The ansatz for the Maxwell field reads: 



where V{r) and aip{r) are the electric and magnetic potentials. 

Note that we include a magnetic part in the vector potential. This is related 
to the fact that the solutions are not static (but they are stationary) . Obtaining 
charged rotating black holes indeed requires electric and magnetic potentials. 

4.3.1 Equations and boundary conditions 

The action (|4.ip with minimally coupled Maxwell fields (|4.14p supplemented 
by the ansatz (|4.8p . (|4.16p leads to a system of non-linear coupled differential 
equations: 




(4.16) 




((3d - 5)/i 



(4.17) 
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^ _ 4(5 - 2d)wa,V + -^b h -^^h 



-2{2d - 5)y'2 + ^(3 - 2d)hw'^ - + {2d - 5)w'^)a'^ + 



(4.18) 



-(l^(^-A + ^^'/^' + ^^/^' 

2 f' 2b r 

-^oi^""^ '2 ^Az3!lh' (^-3)(2d-3)fe^ 
+2(_ + 2d - 5)a, - —^6 -j^^ 

I2{d-i)alh {d-3){d-A)h {d~l)h d-2 d-3 



rf.2 ' J ^ £2 „5 ) ) ^' 



Aw ,2 Aa'^V (d - 3)u;' , 1 b\ f\ 3h' 



(4.20) 



for the metric components and 

,, l,2(d-3) 6' /' h', , h, , , 2(d-3)a^h „ 

< + 2 + -b+J-J^<-b ^""^ + ,4; " = 0' (4-22) 

for the matter fields. Here, the gauge choice 5 = is used and prime denotes 
derivative with respect to r. 

It can easily be seen that the equations of motion present the first integral 

r^"-^^ + ('^ - 3)g, (4.23) 

where q is an integration constant related to the charge of the solution. 
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Thus, similarly to the asymptotically flat case [9T] case, the electric potential 
can be eliminated from the equations (|4.17|) - (|4.22p by making use of the first 
integral 

Note also that the cosmological constant parameter can be arbitrarily rescaled 
by a rescaling of the radial variable r and of the fields h and w. In the following, 
we use this arbitrariness to choose rh ~ 1 without loosing generality. 

Constraint of regularity about the horizon 

We are interested in black hole solutions, with an horizon located at r = r^. 
The solutions can be expanded in the neighbourhood of the horizon leading to 

b{r) = 6i(r-rh) + O(r-r,0', /(r) = /i(r - r,,) + 0(r - r,,)', 
h{r) = ha + hi{r - rh) + 0{r - rh)'^,w{r) = Wh +wi{r - rh) + 0{r - rh)"^, 
a^ir) = aa+ai{r-rh)+0{r-rh)'^, V{r) = Vq + Vi{r - rh) + 0{r - rh)^ 

where ao, Wh, Vq, ho, ai,bi, wi,Vi, hi are real constants. For the solutions to be 
regular at the horizon Vh (or at the cosmological horizon Tc), the equation for 
h leads to the condition Ti{r = rh) = with 

Ti{r) = 8b'h^{12a^ + 7h)+4:xb'hh'{12a^ + 5h) 

- 'S2b'h'^x^ + 8xH'h{f'h-4h') (4.24) 



+ 2hx\—b'h^ - fma')^hw^ + Sa'hwV + b'h' + 5h^{w')^ + ih{VT)) 
2A 

+ h'x''[-=^b'h^ - f{4{a')^hw^ + Sa'hwV - b'h' - h\w')^ + ih{V')^)) 

where we posed a^p = a. The value f'{rh) can be extracted from the equation 
for /, giving 



In the same way, the two Maxwell equations (|4.2ip . (|4.22p lead to a single con- 
dition r2 (r = r/i) = with 



24 



.f'irh) 



Abiho Qrl/P + 8r2 - 12ag - 5/io 

rl Sbiho + rh{4:ho{Vi + aiWo)^ - bihi - h^wf) 



(4.25) 



r2(r) = 4ab'h + r^f'{hw'V' + a'hww' - a'b') 



(4.26) 



4.3.2 Thermodynamical quantities 



Horizon quantities 



Although it is not clear if the thermodynamical properties are well defined 
in the presence of the cosmological horizon, one may still define them in the 
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standard way. The Hawking temperature and the event horizon area of these 
solutions are given by 



Th = 4^, Ah = V,^,ri-\ (4.27) 
47r 



Along with [22], we also write the mass and the angular velocity at the horizon. 
They are obtained by using the appropriate Komar integrals: 



(4.28) 

where Vd^i denotes the area of the id — 2)-sphere. Similar quantities, say 
MctJc can be associated with the cosmological horizon. In fact, the masses 
and angular velocities at the two horizons can be related to each other by 
Smarr-like formulae, as we now demonstrate. 

We consider a Killing vector and integrate the two sides of the identity 



V^V^A^, = -R^-.K, = -8^(rt - -^T5i)K, (4.29) 



over the (truncated) hyper- volume S covering the space-like region between the 
two horizons. Then, using both the Einstein equations and the Stokes theorem 
appropriately result into the identities we are interested in. 
Choosing first K = d^p leads to 



Jc-Ki-2y^(^«a^ + 't^')|r=r. = Ji/-Ki-2y^(«;a^ + l/')|r=r, , (4.30) 

expressing the conservation of the total angular momentum. 

The identity associated with the time-translation invariance {K = 9t) is a 
bit more involved but leads to 



Mc + $cQc - Mh - '^hQh = ^ ' \j ^dr, (4.31) 

where Qh (the electric charge at the horizon [53]) and (the electrostatic 
potential) are defined by 



dO{d^)\r=r,, ^ {V + ^—na^)\r=r,. (4.32) 



The same relations hold at the cosmological horizon for Qc, $c, but with the 
quantities evaluated on the cosmological horizon. Similarly to [91j , the mass at 
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the horizon can further be expressed in terms of the temperature and angular 
momentum according to (recall that D, = w{rh)) 

= \thAh + ^OJff. (4.33) 

The various identities above are fulfilled by our numerical solutions, providing 
several nice crosschecks of our method. 



4.3.3 The asymptotic and global charges 

In this section, we present the global charges characterizing the solutions asymp- 
totically. The metric functions have the following asymptotic behaviour in 
terms of three arbitrary constants a, /? and J 

Kr) = + 1 + ^ + fir) . + 1 + + 0(l/r^-i), 

h{r) = r^il + ^ + Oiiyi-^yi u,(r) = -^ + 0(iy''~^). (4.34) 

The asymptotic expression of the gauge potential is similar to the asymptoti- 
cally flat case: 

= -Jz^ + Oil/r'"-'), a^ir) - ^ + 0(1/^^'^-^). (4.35) 
The mass-energy and angular momentum of the solutions are given by 

^^^(/3-(^-l)«), J^lfi. (4.36) 

The above relations can be obtained by using a background subtraction ap- 
proach (see appendix]^ or using a suitable counterterm technique |94| . 

The electric charge and the magnetic moment of the solutions are given by 

m 

^ = 4nG = 4nG ('-^^^ 

The ansatz in /, b, h has the advantage to present a direct connection with 
the closed form vacuum rotating solution. This is not the case in [55] where an 
isotropic coordinate system is used. 



4.3.4 Numerical procedure and results 

We solved the equations (|4.17p - (|4.22p in the case d = 5 and we hope that this 
case catches the qualitative properties of the pattern of the solutions for d > 5. 
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The equations have been integrated numerically with the ODE solver Colsys 

m- 

The positive cosmological constant leads to the occurrence of a cosmological 
horizon aX r — Tc (with rh < rc < oo) where /(re) = h{rc) = 0. This creates 
a difficulty since Vc constitutes an apparent singular point of the equations. In 
order to overcome this difficulty, we have solved the equations in two steps. 

In the first step, we supplemented the system by the trivial equation dP' jdx = 
and imposed the conditions of two regular horizon at r ~ and r = fixing 
fc, fh by hand and solved the equations for r e [r/j, rc]. The appropriate set of 
twelve boundary conditions at the two horizons then reads 

/ 0, 6 = 0, 6' = I, w = Wh, V = 0, a' = a'^, 

Ti = 0, Ta = for r = r^; (4.38) 

/ = 0, & = 0, Ti = 0, Ta = for r = r^. 

where the functions Ti^2 are defined above. Here we take advantage of the 
arbitrary scale of the field b{r) to impose b'{rh) ~ 1 and of the arbitrary 
additive constant of the electric potential to assume V{rh) = 0. The function 
h{r) will be renormalized appropriately after the second step in order to follow 
the asymptotic behaviour (|4.34p . say h — > h/Nb. Then, the functions w and V 
have to be normalized as well, according to w — > w/A^b, V ^ V/Ni,. 

The values wu and a'j-^ are, a priori, arbitrary and control the total angular 
momentum and the (electric and magnetic) charges of the black hole. The 
numerical value of l"^ is determined by the first step, together with the values 
of all the fields at r = rc. The value of these fields at r = rc and of (."^ can then 
be used as a suitable set of Cauchy data to solve the equations for r S [rc, oo]. 
The disadvantage of the method is that we cannot perform a systematic analysis 
of the solution for a fixed value of the cosmological constant, but fortunately, 
the numerical value of £^ depends only a little on w(rh) and a{rh) once rh,rc 
are fixed. 

In the following, we present the results corresponding to r^, = 1, rc = 3, this 
case corresponds to a large value of the cosmological constant but allows one 
to analyse the effect of A on the solution. We hope that the results for this 
case capture the main features of the solutions for generic values of A. 

A generic profile of the solution corresponding to r/j = 1, rc = 3 , Wh = 
0.6, ajj = 0.5 is presented in figure |4J] it corresponds to l/£'^ = 0.0945 and 
q = —0.08538. The smoothness of the profile at r = Tc can be appreciated on 
the plot. It is also worth to point out that the numerical solutions approach 
the asymptotic behaviour (|4.34p although all boundary conditions are imposed 
at r = rc for the second step of our construction. We supplemented figure |4T] 
with the plot of the t — t component of the metric gtt = b — huP' , showing that 
there is a small region around the event horizon where gu is negative, defining 
an ergoradius where gtt{re) = (on the figure, re ~ 1.17). 
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We managed to construct several branches of solutions for different values 
of a'j^ and increasing the parameter Wh- For fixed a'^, the solutions exist only 
for sufficiently large values of Wh', this can be interpreted as an effect of the 
balance between electric attraction and centrifugal repulsion. This is illustrated 
on figure |4j2] where we plot the asymptotic charges M, J, Q, /i as functions of Wh 
for two values of a^; for aj^ = 0.1 (resp. aj^ = 0.5) the branch exist for > 0.1 
(resp. Wh > 0.465). The figure further suggests that, if a second branch exists, 
its asymptotic mass will be larger than the one of the first branch that we 
constructed. The magnetic moment depends only a little in the parameter Wh- 

In fact, our analysis indicates that the second branch indeed exists: for 
ajj fixed, we find solutions for Wh larger that some minimal value of Wh, say 
uih,min- It turns out that a second branch of solution emerges at the minimal 
value Wh,min- When we increase Wh, we have seen that the second branch ends 
in an extremal solution where f'{rh) = 0. This is shown in figure l473l where we 
present a few parameters that are significantly different in the two branches as 
a function of Wh- Figure WM presents the value of /i, /iq and l/£ as a function 
of il, suggesting that the solution exists in a finite range of 51, for fixed value 
of ajj. When we solved the equation for larger values of Tc, corresponding to 
smaller values of A, we observe that the asymptotic quantities plotted on figure 
14.21 depend weakly on rc. 



1.5 




Figure 4.1: The profiles of a generic solution corresponding to rh = 1, rc = 3 , 
Wh = 0.6, a'^ = 0.5 
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4.4. BLACK STRINGS AND INSTABILITIES 



77 




Figure 4.4: The values /i, ho and l/£ are reported as functions of flh (i-e. 
w{rh), properly normalized) for a'^ = 0.1 and a'j^ = 0.5. 



4.4 Black strings and instabilities 

The simplest higher dimensional black object presenting another horizon topol- 
ogy than 5^-2 is the black string, presenting an horizon topology 5*1 x Sd-3. 
The black string metric reads 

ds^ = -f{r)dt^ + -— + r'^dnl. + dz^, z e 0, L. (4.39) 
fir) 

for some length L of the 5*1 circle and where f{r) = 1 — (■^)'^ ^• 

This is simply a {d— 1)-Tangherlini black hole [53] with an extra Ricci-flat 
direction. Note that the asymptotic spacetinie is locally Minkowskian, more 
precisely, it is of topology Aid-i x Si. 

4.4.1 The black string solution and its properties 

The black string solution has obviously ^ and ^ as killing vectors. Using the 
Komar formula (see appendix|X| , we can compute two global charges associated 
to these killing vectors: the mass M and the tension T. Obers and Earmark 
[55] have shown that for the general black objects (e.g. black strings, localised 
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black holes with or without matter) asymptotically approaching Md-3 x 5*1 
spacetime, the metric component gu and gzz are of the form 

9u--l + ^,+o(j:y \ .9.. = l--££^ + o(^iy (4.40) 

They then established the following formula for the mass and tension of these 
objects: 

M = {{d - 3)q _ c,) , r - ^ (-(d - 3)c, + c) . (4.41) 

The black string solution (|4.39p is characterised by ct — rf^'^ and — 0, 
leading to 

M ^ Y±Z^^d-3)rf-\ r^^ri-\ (4.42) 



The entropy S and the temperature Th are computed in the standard way: 

These quantities are related to each others by the Smarr relation 

(d - 3)M -TL^{d~ 2)ThS, (4.44) 

which can be checked explicitly using the definitions above. Note that the 
Smarr relation can be derived in a formal way, see |96) . The first law of thermo- 
dynamic has an additional contribution coming from the gravitational tension: 

dM = TndS + TdL, (4.45) 

which can also be checked explicitly, considering that the varying parameters 
are Vh and L. 



4.4.2 Dynamical stability 

The dynamical stability of the black string (|4.39p has been studied in 1993 by 
R. Gregory and R. Laflamme |19| . This was partly motivated by the thermody- 
namical instability and from the fact that a black hole is more entropic than a 
black string of same mass when the length of the string is too large as a simple 
calculation shows: 



Sth cx r^-'^^Mbh cx 



(4.46) 
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where Sbs,Mbs (resp. Sbh,Mbh) is the black string (resp. black hole) entropy 
and mass. 

Eliminating the horizon radius, we can express the entropy as a fmiction of 
the mass in each cases: 



/ A/' 



L 



Sbh oc M, 



bh ' 



(4.47) 



Then, comparing the entropy assuming the same value of the mass Mbh 
Albs = AI leads to 

1 

d-4 

(4.41 



Si 



bs 



M~ 



Sbh 



Then, it turns out that if the length is too large for a given mass, the entropy 
of the black string is less than the entropy of the black hole. Note that the 
mass is roughly r'f^^ (for the black hole) so the ratio is essentially proportional 
to rh/L. 

Note also the analogy between this thermodynamical instability of the black 
string and the Plateau Rayleigh instability (a fluid cylinder is unstable if the 
cylinder is too long) |97| . 

In order to consider the dynamical stability of the string, one has to consider 
a small perturbation around the metric (|4.39p : 



(4.49) 



where e is a small arbitrary parameter, gAB are the component of the metric 
(|4.39|) . Has are the component of the fluctuation and A, B = 0,1, . . . , d — 1. 

To first order in e, the equation are given in the transverse-traceless gauge 
{V^hAB = 0, g^^'hAB = 0) by 



&Rab = {^L)'AEh. 



1 



CD 



-- [U5^^5'^ + 2i?i/) hcD = 0, 



(4.50) 



where is the Lichnerowitz operator and the d'Alembertian and Riemami 
tensor are evaluated using the background metric gAB (see Appendix [A|) . 
We consider a perturbation of the form 



hAB = 



( .f{r)Htt Htr \ 

Htr Hrr/f{r) GO 

r'^K 

r^Ksin^e 



■■• 

0/ 



^Q,t—ikz 



, (4.51) 
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where Htt, Htr, Hrr, K are functions of r only and k = 2it/L. The values of fl 
and k are not independent as will be shown later in a dispersion diagram. 

It is possible to reduce the linearised Einstein equations given by (|4.50p to 
a single equation, eliminating all but one field, say Htr following [1^1, using the 
gauge conditions and the Lichnerowitz equations. The system then reduces to 
a single second order ordinary differential equation. 

The reduced equation is of the form 

K, + Q(r)i/;, + Sir)Htr = 0, (4.52) 

where of course, Q{r),S{r) depend on k,il,r,d and can be found in [98l US] 
between many others. The asymptotic solution of the equation is given by 

Htr oc e±'-^^^+^ (4.53) 

while close to the horizon, the field Htr behaves like 

i/t, oc (r-rh)-i±*5. (4.54) 

Gregory and Laflamme found that the conditions required in order to describe 
well behaved perturbation [981 [W[ [M] leads to the requirement that close to 
the horizon, Htr behaves like (|4.54p with the + sign and tends to zero asymp- 
totically. They found solutions to (|4.52p with 51 > for small values of k, 
which means that black strings are unstable towards long wavelength (small 
wavenumber) perturbations. However, larger values of k lead to f2 < 0, the 
black strings are then stable toward short wavelength perturbations. 

This is illustrated in figure where the value of as a function of k is plot- 
ted. Note the presence of a zero mode; i.e. a solution with k ^ 0, say fc^, fl = 0. 

4.4.3 Phases of asymptotically locally flat black strings 

The existence of the zero mode implies that the particular perturbation with 
wavenumber k = k^ '^s, static. This is a hint for a phase transition. 

It was widely believed that the string would decay to a localized black hole, 
i.e. a black hole in a spacetime which is asymptotically A4d-i x Si and with 
horizon topology Sd-2- However, it has been shown by Horowitz and Maeda 
that such a phase transition would take an infinite proper time at the horizon 
|100| . This led to postulate an intermediate phase, namely the non uniform 
black string phase. 

The non uniform phase is characterized by the fact that the solution depends 
non trivially on z; as a consequence, ^ is no longer a Killing vector. However, 
^ is still a Killing vector in the asymptotic region of the spacetime, allowing 
the definition of the tension. The horizon topology of the non uniform black 
string is still Sds x 5i, in contrast with the localised black hole characterise 
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by an horizon topology 5*^-2 (the asymptotic spacetime is Md-i x Si in both 
cases) . 

The non uniform black string solution is a static solution to the vacuum 
Einstein equations and have been studied first perturbatively by Gubser |81| 
and in the full non linear regime by Wiseman [82] • The non uniform phase 
is connected to a localised black hole phase; at the connection point (merger 
point), there is a horizon topology change. The localised black hole phase 
have been studied by T. Wiseman and Kudoh [8211101) and the resulting phase 
diagram connecting the localised black hole, the non uniform black hole and 
the uniform string has been studied in [lOlj . The phase diagram is presented 
in figure for d = 5,6. The lines represent the various static solution which 
should be equilibrium configurations since they don't depend on time. 

It should be noted that various 'problems' appear in this context; for ex- 
ample, it is not clear that cosmic censorship (which roughly state that every 
singularities should be hidden behind an horizon) is not violated. At the merger 
point, the string disconnects and form a localized black hole; precisely at this 
point, the areal radius (square root of the coefficient of the angular part of the 
metric) of the horizon vanishes leaving the singularity naked. Another issue is 
that the horizon topology change from Sd-3 x 5*1 to 6*^-2; it is indeed possible 
to study these solutions starting from the black hole side or by the string side, 
but not both in once. This of course makes the picture quite involved |101| . 

The transition from black strings to black hole is also from considerable 
interest from the thermodynamical point of view. For instance, the non uniform 
string is less entropic than the uniform string for small deformations, up to 
certain number of dimensions (13). This is called the critical dimension |102j . 
We will come back to the question of critical dimension in the next chapter; we 
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refer the interested reader to the review from Kol |103| for a deeper discussion 
on phase transition between locahzed black hole and black strings. 




Figure 4.6: The fi — n phase diagram, n being the relative tension TL/M and 
/u = M/L'^~'^ is the rescaled mass. Reprint from [SB]- The straight line is the 
uniform string phase, the line going down and emanating from the uniform 
phase is the non uniform phase while the line emanating from the origin is the 
localized black hole phase. 



4.5 Black string solutions with a positive cosmo- 
logical constant 

Considering (j4.ip with a positive cosmological constant and restricting to pure 
gravity (Lm = 0), we can try to construct the generalization of the black string 
solution presented in the previous section. This was actually already done in 
[20) for a negative cosmological constant, but we will come back to the AdS 
black strings in more details in the next chapter. 

In order to construct a black string solution with a cosmological constant, 
we supplement the action (|4.ip by the following ansatz 

ds^ = -h{r)dt^ + — + r'^dnl_^ + a{r)dz^. (4.55) 

where d^l'^_^ denotes the square line element on the unit {d — 3)-sphere. 

The reason why we add the new field a can be intuitively understood by 
the fact that the cosmological constant gives a curvature to the extra-direction. 
From a mathematical point of view, the ansatz without the function a doesn't 
lead to a consistent set of equations. 
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4.5.1 Equations and boundary conditions 

The resulting equations of motion are given by 

^, ^ 2(d-4) 2(d-l)r 2(d-4)/ /a' ^ b' 
r P r \ a b 

^„ _ {d-3){d-4)b (d-3)(d-4)b {d-l){d-A)b {d - 3)ba' 

j,2 J.2 j- pj 

{d-A)b' {d-4)b' {d-l)rb' a'b' 6'^ 



f ' P2f ■ 9 ■ ^ ' (4.56) 

^ ^ ^b[~£'id-3}{d-A)il- f) + id-l)id-2y]+id-3)r£^fb' 
a 



~re^f[rb' + 2{d-3)b] 

We supplement these equations by the boundary conditions corresponding to 
a black string: 

Hrh) = /(r,0 - 0, b'in,) = 1, a{rk) = 1, (4.57) 

where we set the normalisation of a, b such that a(rh) — 1, b' [rh) = 1 without 
loss of generality since the normalisation can be absorbed in a redefinition of 
the time and z coordinate. The conditions on b{rh), f{fh) are necessary to 
implement a regular horizon at r = r/j. 

It is interesting to notice that the above equations may be written in a 
decoupled form. The decoupled form of the equation is valid for both A > 
and A < 0. In the following, we define e = — A/|A|. 

First, let us introduce 

. , X ra' , , rb' , , 

A{r) = — B{r) = — 4.58) 
a b 

in terms of which the equations for /, a, b rewrite: 

r/' = 2{d - 4) {k - /) - 2{d - 1) J - /(r) (A + B) , (4.59) 

rB' f = {d-l){2- B)—r - {d- ^)kB, (4.60) 

{AB + 2{d-3){A + B) + 2(d-3)(d-4))/ (4.61) 

= 2k{d-i){d-4)-2{d-l)(d-2) — . 

Solving the third equation for A = A{B,f), substituting A{B,f) in the first 
equation and solving for B = B{ f, /') and finally evaluating the second equa- 
tion with A = A{B{f,f'),f), B = B{f,f') leads to the following decoupled 
equations: 
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r'ff 



r//' + (r/T + 2((d-l)-(rf-4) — 



(4.62) 



2(d-4) (d- 1) - (d-4)— /-3 fc(d- 1)_ (d-4)— rf 



B = -B (4.63) 

± ^i.--2(.-3,((.-4)(i-i)-if)-. ^--'^-;y-j Z 

where B := ((rf - 4) + \^ + (d- 4)i - (d- l)i^) and 



^ -2(d - l)(d - 2)^ + 2(d - 3)(d - 4)/ + 2(d - i)fB 

2{d-3)f + fB ■ ^^-^^^ 

Note that the decoupUng is still valid in the more general case considered in [20] 
where topological black holes are investigated as well (i.e. space-times where 
the spherical part dfl^_-^ in ()4.55p is replaced by the metric of an hyperbolic 
or flat manifold with the same dimensions). 



4.5.2 Asymptotic solution 

We have not been able to solve ()4.63p explicitly, but we can study the asymp- 
totic behaviour of / with this equation. In fact, to the leading order, it appears 
that 

/(r) « For"* (4.65) 

for a constant Fq and for any G {2}U]4, oo[, so we cannot conclude on the 
asymptotic behaviour of / at this stage. In fact, it is not possible to fix ana- 
lytically 0, but it is possible to give a description of the asymptotic behaviour 
of the metric functions in terms of only one parameter. 

We have to consider two cases : = 2 and > 4. In the case e = +1, 
where the cosmological constant is negative, the asymptotic solution admits 
= 2 [20]. In the case e = — 1, there are strong numerical evidences that </) = 2 
leads to a solution where the metric functions diverge at the event horizon |10) . 
Since this section is precisely devoted to the case e = — 1, we will then assume 
(t)>4. 

Assuming a{r) oc r", b{r) oc and / = F^r'^ for large values of r, it 
is possible to express a and as a function of /3 only, i.e. to describe the 
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asymptotic behaviour of the metric functions with only one parameter: 

, /32 + 2(d-4)/3 + 2(d-3)(d-4) 



/ -> For"' , 
A{r) a(^) = - 



/3 + 2(d-3) 
2(d-3)(d-4) + 2(d-3)/3 



/3 + 2(d-3) 

B{r) P, (4.66) 

for large values of r. 

Note the particular relation between the exponents: 

a + /3 + + 2(d-4) (4.67) 

Since these behaviours are compatible with the equations of motion for each 
value of P such that ^ > 4, it is not possible to fix /3 without more assumptions, 
but there are strong numerical evidences for a(r) = 6(r) at the leading order 
in the asymptotic region. This means that the value of /3 chosen by the system 
is a fixed point of a(/3) and is given by 



P = -2{d-i)- (4.68) 

We have constructed the higher order corrections for the function /(r) and 
obtained 

/(^) - - + j4tIu- '« - 

with (j) discussed previously, cf) = ^^^'^^p+2{d-t}^'''^'^^'^'' ■ 

Note the subleading term — r^/^^; it is the term we expected to the leading 
order to describe an asymptotically locally de Sitter spacetime. As a conse- 
quence, it is obvious that the solution does indeed not describe an asymptot- 
ically locally de Sitter spacetime. In the next sections, we will evaluate the 
energy and try to determine whether this solution can be physically relevant. 

4.5.3 Energy of the solution 

The energy of the solution is given by the following quantity |104j : 

E = -^f N{Kd-2-Ko.d^2) (4.70) 

where is the lapse function, is the boundary of the fix time sliced manifold 
Ej, Kd~2 is the trace of the extrinsic curvature of 5^ in Et and Ko.d-2 is the 
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trace of the extrinsic curvature of the fixed time boundary of some reference 
background spacetime (see appendix |A|). 

The trace of the extrinsic curvature of the r = Const, hypersurfaces, say 
Kd-2 is given by: 

This quantity, once integrated according to (|4.70p . is not divergent for r — > c» 
despite of the unusual asymptotic form of the solution. It follows that we 
don't need a reference background in order to have a finite energy. Anyway, it 
wouldn't have been not clear which background to use with such an asymptotic. 

We can now compute the energy of the solutions which exhibits the asymp- 
totic behaviour as described in the previous section. Using (|4.67p . the energy 
is simply given by 

E = -l-Vd-3L^2Foid - 2){d - 3) (4.72) 

where L is the length of the extradimension in the z direction and V^s is 
the surface of the unit d — 3 sphere. This energy is finite and depend only 
on _Fo which depends essentially on the dimension of space-time and of the 
cosmological constant. Note that _Fo is determined numerically. 



4.5.4 Geodesic equations and curvature invariant 

We will now compute some curvature invariants in order to check whether the 
solution with the asymptotic behaviour described above have a chance to be 
regular. The Kretschmann curvature invariant is given by: 

K ^ R.,.dR^-- - '^-f : - '^f' + ^'-f" . (4.73) 

Plugging the asymptotic behaviour of the function / in this expression shows 
that K is diverging for large values of r: 

K ^ {d- 3)i^o (2(rf - 4) + 0^) r^"^"'' ^ oo for r ^ oo, (4.74) 

since > 4 by assumption. This shows that the solution with the asymptotic 
behaviour of previous section is asymptotically singular. 

Worst, the asymptotical singularity can be reached in a finite proper time 
for a massive or massless observer (say with mass m) starting from outside the 
horizon. This can be seen by integrating the geodesic equation at fixed 9i (the 
angular sector of the metric). 
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Using the fact that k constant along geodesies, provided fc is a KiUing 

vectoiQ, the geodesic equation can be written in the form jl7| : 

where £, Z are conserved quantities along the geodesic, related to the observer's 
energy and momentum along the z direction and dot denotes derivative with 
respect to an affine parameter along the geodesic, say A. This equation can be 
integrated numerically with particular solutions admitting the singular asymp- 
totic behaviour. One can also assume that the observer is in a region well 
described by the asymptotic solution and integrate the geodesic equation using 
a = b = bor^ and / = For"^ to first order in r, with (3 defined above and obtain: 

It turns out that 1 + < for all values of > 5. Consequently, in the limit 
r — >■ oo, A is finite. In other words, an observer far from the horizon can reach 
the asymptotic singularity in a finite affine time. A numerical computation for 
observers that start close to the horizon reveals that the singularity can still 
be reached in a finite time for a massive or massless observer. 



4.5.5 Numerical results 

We solved numerically equations ()4.59p . ()4.60p and (|4.6ip with the boundary 
conditions (|4.57|) for many values of the horizon r^. We present the evolution 
of the ratio rf / f for A > black string solutions and for several values of d 
and = 0.5 in figure HTTl The figure clearly demonstrates that the power law 
configuration (|4.66p is approached. The exact values of the exponents (|4.68|) 
coincide with our numerical values within the numerical accuracy required for 
the numerics, i.e. typically 10~*. 

We also investigated the solution in the interior region, i.e. for r < r/i, to 
check whether there exists a second horizon which would 'hide the asymptotic 
singularity', but it turned out that this was not the case. We have a naked 
singularities both, at r = and at r = oo. We are tempted to interpret this 
solution as an hypercylindre with an horizon at some equator by matching 
the origin with the point at infinity by means of an appropriate system of 
coordinates, still to be found, but that we expect to exist. 

^In order to prove that, one must use the geodesic equation: i^Vai'' = 0, dot denoting 
derivation with respect to an affine parameter parametrising the geodesic, and the Kilhng 
equation: Va^t + V(,fca = 0. 
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Figure 4.7: The ratio rf / f is given for several values of d 

4.6 Concluding remarks 

In this chapter, we quickly reviewed black hole and black string solutions avail- 
able in higher dimensional general relativity. The analytically known black 
hole solutions are either rotating and charged, either rotating with a cosmolog- 
ical constant, either charged with a cosmological constant. The missing part 
was the charged rotating with cosmological constant solutions. Such solutions 
have been constructed numerically for a negative cosmological constant in [87] 
and we constructed their counterpart for a positive cosmological constant. The 
main difficulty in the case of a positive cosmological constant is the fact that 
there exists a cosmological horizon where some metric coefficients become zero. 
We imposed regularity constraints at the cosmological horizon and added an 
equation, expressing the fact the cosmological constant is indeed a constant in 
order to have sufficient freedom to impose the constraints. 

We also reviewed the black string solution, in flat spacetime and with a pos- 
itive cosmological constant. We left the AdS black string for the next chapter. 
We briefly introduced the Gregory-Laflamme instability and the phase diagram 
of the black string in flat spacetime. We presented black string solutions with a 
positive cosmological constant; an interesting result is that there doesn't exist 
an asymptotically de Sitter black string, the spacetime becoming asymptoti- 
cally singular. This can be interpreted as the fact that the compact direction 
cannot support the pressure induced by the cosmological constant. 



Chapter 5 



Phases and thermodynamics 
of AdS black strings and AdS 
non uniform black strings 

In this chapter, we consider the black string solution in asymptotically lo- 
cally AdS spacetime. First, we review the AdS uniform solution and revisit 
some particular thermodynamical properties of the latter. These uniform black 
strings in AdS have been considered by several authors in different contexts 
[201 11051 \n\ . Then we investigate the stability of the AdS uniform black string 
by looking for the existence of the counterparts of the Gregory-Laflamme zero 
modes. This approach was introduced by S. S. Gubser [ST] for asymptotically 
locally flat black strings. It consists in considering the non uniform solution as 
a perturbation of the uniform solution and in solving the first order equations 
in the perturbation. 

In a third step, we present solutions to the non uniform problem, first in a 
perturbative approach then in the full non linear regime. Whereas the uniform 
black string does not depend on the extra-coordinate, the non uniform black 
string is characterized by a non trivial dependence on this extra-coordinate. 
Up to our knowledge, this is the first time this technique is applied to asymp- 
totically locally AdS spacetimes. 

The second order in the perturbative approach is already sufficient to see 
that the cosmological constant provides drastic changes in the phase diagram 
of AdS black strings. In particular, it suggests the existence of a thermo- 
dynamically stable phase of non uniform black string in a regime where the 
corresponding uniform black string with same length is thermodynamically 
unstable. This is one of the unexpected features of the asymptotically locally 
AdS spacetime. 
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Finally, the solution to the system of partial differential equation provides 
hints for the existence of a localized black hole phase and for the possibility 
that the phase transition between the uniform and non uniform phase can be 
of second order in every number of dimensions, provided the horizon radius is 
large enough. The AdS localized black hole solution has not been constructed 
yet, up to our knowledge. 

This chapter presents the counterpart of the black string (uniform and non 
uniform) which we briefly discussed in the previous chapter. Our original results 
have been presented in [H [Ml 123 [TUB] . 

5.1 The generic model 

We consider the d-dimensional Einstein Hilbert action with the presence of a 
(negative) cosmological constant and the Hawking-Gibbons boundary terms. 



where G is the d-dimensional Newton constant, g the determinant of the metric, 
R the scalar curvature computed with the metric gAB, A,B = 0, . . . ,d — 1 
associated to the manifold Ai, A the cosmological constant, h the determinant 
of the induced metric on the boundary dM of the manifold and K the extrinsic 
curvature of the boundary manifold. 

In order to investigate the various phases of asymptotically locally AdS 
black strings, we consider the equations resulting from the variation of (|5.ip 
with respect to the metric components. We then study the equations for metrics 
of the form 



ds^ = -e2^('--^)6(r)dt2 + e2^('''^) ^ + air)dz' + e^'^^^^^'^rH^ll,^, (5.2) 



for some functions a, &, /, A, _B, C and where dfi^-s is the line element of the 
(d — 3)-unit sphere. We impose an horizon radius at a given value of the 
radial coordinate, say rh- The variable r then ranges from rh to oo while 
the coordinate z ranges from to L G R+. We furthermore define the AdS 



The metric ansatz (|5.2p is a generalization of the uniform black string metric 
(|4.55p of chapter IH we recover (|4.55p ior A = B = C = Q. The functions 
A{r^ z), i?(r, z), C(r, z) parametrise the non uniformity and we will refer to these 
as the correction functions (although they are more than 'just' corrections in 
the non perturbative approach presented in the last sections of this chapter) . 




(5.1) 
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In the presence of a non vanishing cosmological constant, the field equations 
imposes 

bir)^'^ + 0{rf , /(r) = ^ + 0(r)° , a(r) = ^ + O (r)° , (5.3) 

to the leading order for ?' — > cx) for any value of the z coordinate. This asymp- 
totic behaviour describes an asymptotically locally AdS space (Appendix [B|) . 
The correction functions vanish at infinity; if it was not the case, the asymp- 
totic spacetime wouldn't be of the form (|5.3p . It is worth noting that the decay 
of the correction functions cannot be arbitrary, in particular it shouldn't spoil 
the asymptotic form of the metric components. This will be discussed in due 
course. 

The ansatz (j5.2p is a general axisymmetric parametrisation (in a certain 
gauge); the functions A, B, C parametrise the z component dependence of the 
line element whereas a,b, f are relevant for a uniform setup. The need of the 
function a(r) was already discussed in the previous chapter, in the black string 
with positive cosmological constant; this is in contrast with the asymptotically 
locally flat case where the extra-coordinate z was Ricci-flat. 

In the next sections, we will consider the uniform string phase, its sta- 
bility, the non-uniform string phase in a perturbative approach and in a non 
perturbative approach. 

From now on, we will omit the coordinate dependence of the various metric 
functions. 



5.2 Uniform string phase 

The uniform phase is characterized by 

A = B = C = 0, (5.4) 

since the functions A, B, C induce a non trivial dependence on the z coordinate 
when non vanishing, making the string non uniform (along the z direction). The 
line element of the uniform black string then reduces to 

ds^- = ^b(r)dt^ + -— + r^dnjo + a(r)dz^. (5.5) 
/('') 

Unfortunately, there is no analytic solution describing AdS c?-dimensional black 
strings, so one has to rely on a numerical resolution of the Einstein equations, 
except in some specific case, where other fields are present (for example with an 
additional 2-form and particular relations between the coupling constants, see 
[105j ). Solutions for the asymptotically locally AdS black string have been first 
constructed in |20| ; we will review the main properties of the uniform string in 
the following section. 
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5.2.1 The equations and boundary conditions 



The variation of the action (j5.ip with respect to the metric components, sup- 
plemented by the ansatz (|5.5p and (|5.4p leads to a coupled system of three non 
linear ordinary differential equations : 

^, ^ 2(d-4) ^ 2(d-l)r 2(d-4)/ ^ /a^ ^ b' 
r P r \ a b 



b" 



{d-'d){d-A)b _ (d-3)(d-4)& _ (d-l)(d-4)6 ^ (rf - 3)6a^ 

j,2 J.2 J £2 J- j,^ 

(d-4)&' {d-4)b' {d-l)rb' a'b' b'^ 
b[e{d - 3){d - 4)(1 - f) + {d- l){d - 2)r2] - (d - 3)r£2/&' 



a 



rPf[rb' + 2{d-3)b] 



(5.6) 



Boundary conditions 



In order to describe black objects, we impose an event horizon on a given value 
of the radial coordinate, say by setting 6(r/,,) = 0. The regularity of the 
equation then imposes f{rh) = 0. Since the system of differential equations 
has a total degree 4, we still need two boundary conditions. It is important 
to notice that the equations (|5.6p are invariant under independent rescaling 
of the functions a and b. So we choose to arbitrarily fix a{rh) = b'{rh) = 1 
and globally rescale the a and b functions a posteriori in order to obtain the 
correct asymptotic a = r"^ / f- + O (r)" , b ^ r'^ j f- + O {rf. Note that the 
equations (|5.6p are not invariant under a rescaling of the / function, moreover, 
once f{rfi) = is fixed, the asymptotic behaviour of the / function is naturally 
to the leading order. 
To sum up, the boundary conditions read 

/(r,0 = &(nO=0 , a{rh) = b'{rh) = l. (5.7) 

for a given value of the horizon radius r/j. 

Note that there exists a solution of (|5.6p regular on r e [0, oo]. However, in 
this case, regularity conditions and thus boundary are different: 

/(O) = 1, b{0) = bo, a(0) = ao, b'{0) = 0, (5.8) 

for real values of op , bo chosen such that the solution is asymptotically locally 
AdS. 

The black string solution obeying (|5.7p approach the regular solution for 
rh — > 0; the convergence is point-like for r g]0,oo]. Further details on the 
globally regular solution can be found in |20| . 
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5.2.2 Near horizon and asymptotic expansion 

In the region where r r^, the solution can be expended as a Taylor series 
according to 



2ah {d - 1) Th 



\2 



2 



2 



0[r) On{r i ^ {d - ^) Pr^ + {d - 1)tI 2 ^ ^ ' 



where ah = a(rh), hu = b'{rh) are real constants. 

The asymptotic behaviour of the background functions is given in [20] : 
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f] \ 2i / /} \ d—3 / n\ ti— 3 



+ (ct + + Co) ( - j + ^ ^d^afc+iC log 



(5. 



where 5 is the Kronecker delta, [XJ stands for the floor integer value of X , 
Ct,Cz are parameters to be determined numerically and ai,fi,co,S. depend on 
d and are given in |20) . 

For fixed values of and £, the black string solution extrapolating between 
and (|5.10p has very specific values of ah,bh. 



5.2.3 Thermodynamical properties of the uniform black 
string 

The Hawking temperature of the string can be found by computing the surface 
gravity or by demanding regularity in the Euclidean section near the horizon. 
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Both procedure give the same result: 




(5.11) 



The entropy of the black string is computed using the Hawking-Bekenstein 
formula and is given by one quarter of the horizon area: 



where Vds is the surface of the unit {d — 3)-sphere. 

The quantities defined in the asymptotic region, namely the mass and ten- 
sion require more attention. Standard approaches like Komar integrals (see 
appendix [X]) are not applicable because the action is diverging asymptotically. 
However, it is possible to use a regularising procedure; such a procedure has 
been first proposed by Balasubramanian and Krauss |107) and was inspired 
by the AdS/CFT correspondence. This procedure consists on adding to the 
action ()5.ip appropriate boundary terms Id, which are functional only of cur- 
vature invariants of the induced metric on the boundary. Such terms will not 
interfere with the equations of motion because they are intrinsic invariants of 
the boundary metric. By choosing appropriate counterterms, which cancel the 
divergences, one can then obtain well-defined expressions for the action and 
the energy momentum of the spacetime. This procedure is intrinsic to the 
spacetime of interest and it is unambiguous once the counterterm action is 
specified. 

Thus we have to supplement the action (|5.ip with |107] : 



where R and R"'' are the curvature and the Ricci tensor associated with the 
induced metric 7 on dAi and Q{x) is the step function such that (a;) = 1 
provided x > 0, and vanishes otherwise. 

The series truncates for any fixed dimension, with new terms entering at 
every new even value of d, as denoted by the step-function. 



(5.12) 




(5.13) 
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However, given the presence of \og{r/£) terms in the asymptotic expansions 
(|5.10p (for odd d), the counterterms ()5.13p regularise the action for even di- 
mensions only. For odd values of d, we have to add the following extra terms 
to (EH) [TUg] : 

^ct = fd'^-'xV^logi-jL.J-^il-R'-Ra.R'^'^) (5.14) 



SttG J ^ ' { ' 8 '3 

/S / o -1 

"128 \ 95 ^ ^'^ r<ac6(i — Y^r< VaVf,K 

where □ denotes the d'Alembertian computed with the boundary metric ^ab- 
Using these counterterms in odd and even dimensions, one can construct a 
boundary stress tensor free of divergence from the total action I = S + 1^^.+ /^^ 
by defining a boundary stress-tensor: 



(5.15) 



Thus a conserved charge 



f d'^"'^" ^'Tab, (5.16) 

can be associated with a closed surface E whose normal unit vector is noted n", 
provided the boundary geometry has an isometry generated by a Killing vector 

If ^ = d/dt then is the conserved mass/energy M. Similar to the A = 
case |109) . there is also a second charge associated with d/dz, corresponding 
to the solution's tension T. 

For example, in the case d = 6, the components of the boundary stress 
tensor reduce to 

Ttt = TTr7^(-4ct + c.)4> (5-17) 



1 £■ 
(-Ct -I- 4c;j) ■ 



■2 



1 



where 9, 0i, 02 are angle of the angular sector. 

After some algebra, the mass (associated with dt and with E chosen as 
the asymptotic cylinder Sd-s x Si) and tension (associated with dz and the 



96 



CHAPTER 5. PHASES OF ADS BLACK STRINGS 



asymptotic sphere Sds) are given by 

M = Mo + Me, Mo = (c, - (d ~ 2)ct) , 

IottG 

r = To + Tc, To = ^'^g J^'' (id - 2)c, - , (5.18) 

and where Afc = ~LTc = ^ i^^G^ ^ (12^^,5 — i^^'^d,? + . . .) are Casimir-Uke 
terms in the dual theory defined on the boundary spacetime |20| ; however, 
they are not relevant for our purpose. 

Note that the thermodynamical quantities defined above obey a Smarr re- 
lation [10] 

M + LT ^ThS. (5.19) 

This relation can be obtained by evaluating the action with suitable boundary 
terms, using the equations of motions: i? — 2A = 2i?( or i? — 2A = 2i?|. But 
since dt and dz are Killing vectors, Rf and i?| can be expressed as divergences. 
For example, in d = 6, we find: 

^ -I i^'"'^)'- ^ -I ('•'"'y'l) ■ 

Then, using the suitable counterterms in d = 6, we find that the total action 
I evaluated using i?* is given by 



{-Act + cz)t' + riJan^ [{d ~ 4)rl + {d - 1)£2] \=M- ThS, 



V3L , __3 /_ bh , ,, 2 

-r l^u, — _L;t J 1 — 1,^ - 

(5.21) 



while the total action evaluated with R^ is 

/-^(-4c. + c.) = -rX. (5.22) 

It then follows that M + TL = ThS, as stated previously. The computations 
are similar in higher number of dimensions. 

The solutions also obey the first law of thermodynamic: 

dAI = TndS + TdL. (5.23) 



5.2.4 Integration of the first law and the string tension 

As described above, the mass and tension are defined in terms of quantities 
appearing in sub- leading orders of the asymptotic solution (|5.10p . In practice, 
these quantities are difficult to extract with a sufficient precision, especially 
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when the number of dimensions is large. An alternative is to integrate the 
first law of thermodynamic (|5.23p and to use the Smarr relation (|5.19p . This 
procedure has been used previously in the context of non uniform black strings 
in asymptotically locally flat spacetime by Kudoh and Wiseman |101| : it is 
however the first time it is used in asymptotically locally AdS black strings 
and with numerical solutions. 

Since the length L doesn't enter explicitly in the equations of motions, we 
can assume working at fixed length, without loss of generality. Moreover, the 
length appears as a multiplying factor in the mass and the entropy and doesn't 
enter the definition of the tension nor of the temperature. So we choose to 
work with mass and entropy per unit length. 

The first law then reduces to dAI/L = TudS/L. In practice, for a given 
value of the cosniological constant, Th and S/ L are thermodynamical quantities 
defined on the horizon and we compute them as functions of the horizon radius. 
So one can construct Tnirh)-, S{rh)/ L and integrate the first law in order to 
find the mass as a function of the horizon radius: 

M{rh) _ M{rl) 1 , ,dS{x 



L TH{x)^dx. (5.24) 



Of course, we still need one value for the mass, M{r'^)/L, which we can be 
computed in a region of the parameters where the extraction of the numerical 
coefficient is still manageable. Moreover, the effort is to be done only once for 
each value of the number of dimensions. 

Once the mass is computed, we use the Smarr relation (j5.19p to construct 
the tension: 

T-/ X TH{rh)S{rh) - M{rh) 

T{rh) = . (5.25) 

From the thermodynamical quantities, one can compute whether a given 
uniform black string is thermodynamically stable or not. The object needed to 
deal with thermodynamical stability is the specific heat, defined as 

dS 

Cp = Th^, (5.26) 

where the length L is kept fixed. If Cp > 0, we have thermodynamically stable 
solutions while if Cp < 0, we have thermodynamically imstable solutions. Pos- 
itive specific heat means that increasing the temperature increases the entropy 
and conversely for negative specific heat. Since Th is always positive, it is 
sufficient to look at the sign of 
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5.2.5 Scale invariance 

Note that the equations (|5.6p are invariant under the following scaling trans- 
formation 

£ = X£^ f = Ar, (5.27) 

for a real constant A. 

This transformation has an effect on the various thermodynamical quanti- 
ties as follows 

S = X'^-^S, Th = Th/\. (5.28) 

We take advantage of the scaling relations (|5.27p to vary the AdS radius 
and work with a fixed horizon radius, since it is technically easier to keep the 
horizon radius fixed. We set arbitrarily r?i = 1 and varied the AdS radius 
i. Afterwards, we rescaled all the thermodynamical quantities in order to 
investigate the physics of the black strings at fixed AdS curvature, which seems 
to make more sense physically. 

5.2.6 Numerical solution 

The profiles of the a, 6, / functions are generically the same for every number 
of dimensions and are presented in figure 15.11 for ^ = 1 . 

70 
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Figure 5.1: Profile of the metric functions a,b,f for d = 6,rh = 0.5, £ = 1 
(reprinted from [20]). The profiles are qualitatively the same for every num- 
ber of dimensions, difference arising at sub-leading terms in the asymptotic 
development. 
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Thermodynamical quantities defined in the previous section are computed 
as functions of rh and are presented in figure [5^ for d = 6,i = 1. Once again, 
the general case with d arbitrary presents the same main features than in 6 
dimensions. 




Figure 5.2: Temperature Th, mass M, tension T and entropy S for d = 6,£ = 1 
as functions of the horizon radius. Here again, the figure is reprinted from |20) . 

It should be noted that the entropy is an increasing function of the horizon 
radius while the temperature first decreases for small values of rh and start 
increasing about = rf^ {rf^ ~ 0.5 for d = 6,^ = 1). In other words, there is a 
phase of uniform black strings where the entropy increases when the tempera- 
ture decreases, i.e. Cp < and another phase where Cp > 0; these phases are 
thermodynamically unstable (resp. stable) and are referred to as small (resp. 
big) AdS black strings. This phenomenon is similar to the AdS black hole case 
[88]. 

Moreover, there exists a regime where the tension is negative. This is an 
unexpected feature of the AdS uniform black string and has been pointed out 
in |22| . as we explain now. 

A regular solution in the asymptotically locally flat spacetime would be 
the {d — l)-dimensional Minkowski spacetime times a circle, Aid-i x S^. Of 
course, this solution has zero mass and zero tension. The d dimensional uniform 
black string in the asymptotically locally flat case is characterise by a positive 
mass and tension. However, the AdS counterpart of the regular solution is 
characterized by a positive non vanishing mass |20| , a temperature and entropy 
such that ThS = 0. Using the Smarr relation, it follows that the tension of the 
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regular solution is negative: 

M 

M + TL = O^T= — -. (5.29) 
1^ 

It turns out that the black strings with small values of also have negative 
tension since these black strings are continuously connected to the AdS regular 
solution. This is confirmed by a numerical analysis as shown in figure 15.31 and 
15.41 for d ^ 6,i = 1 and d = 7,i = 1. These values were obtained using the 
definition (|5.18p of the mass and tension as functions of ct , Cz (in units where 
the Newton constant is set to one) and are in accordance with the first law 
integration procedure proposed in section [5. 2. 41 Note however that the regular 
solution {rfi ~ 0) has a vanishing total energy since E ^ M + TL = 0, where 
E denotes the total energy [5S] (intuitively, M is the energy due to the mass 
of the black string and TL is the binding energy of the string) . 

It should be stressed that uniform black string solutions in asymptotically 
locally AdS spacetime have been constructed previously in different contexts 
|106i 11051 [2T] (e.g. with or without gauge fields). Our solution and thermo- 
dynamical quantities are compatible with these previous results; however, the 
negative tension phenomenon was not noticed in these references. This phe- 
nomenon should hold for every number of dimensions, provided the mass in 
the r/i — > limit does not vanish. 




Figure 5.3: The mass and tension per unit length for the uniform black string 
in d = 6,£ = 1. Large values of exhibit a positive tension while smaller 
values of rh have negative tension. The figure on the right is a zoom of the 
region where the tension is negative. 

The negative tension is an effect of the AdS curvature since this phe- 
nomenon doesn't appear in the asymptotically locally flat uniform string prob- 
lem. This might be interpreted as follows: on the one hand, a negative cos- 
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Figure 5.4: This figure is the same as 15.31 for d = 7. 

mological is characterized by a positive constant pressure. On the other hand, 
tension on strings typically acts oppositely to the pressure. In this picture, 
when the gravitational tension is too small, the pressure fi:om AdS dominates, 
resulting in a negative tension. If the string is massive enough, the gravita- 
tional tension is higher and then dominates the AdS pressure. The effect on 
the sign of the tension then results from a balance between the AdS pressure 
and the gravitational tension. 

We checked the sign of the tension numerically for d = 6, 7 but we expect 
this property to hold for every value of the number of dimension. 

5.3 Stability of the AdS uniform string 

From the black string instability in asymptotically locally flat spacetime (see 
section IT4| . it is expected that the AdS uniform black strings are also unstable, 
at least for small values of the cosmological constant. In order to investigate 
the stability of the AdS black strings, we consider a static perturbation around 
the string by choosing the functions A, B, C in (|5.2p as: 

A ^ eAi(r) cosfcz + 0(e)^, 

B = eBi(r)cosfcz + 0(e)^ (5.30) 
C = eCi{r)coskz + 0{e)'^, 

where e is a small real parameter and fc = ^ is the critical wavelength. We 
expand the Einstein equation in a series in e and keep the equations appearing 
at order e. 

Note that all the wavelength should not lead to a static perturbation, so k 
will have to assume a precise value depending on the parameters of the model 
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in order to lead to a static perturbation. This precise value of k will be the 
AdS counterpart of the critical Gregory-Laflamme wavenuniber. 

Note also that once the length L is fixed, one can say whether the black 
string is stable or not: for a given value of L and a given perturbations with 
wavenumber K, the perturbation cannot fit in the compact extra-direction if 
if < 27r/L; by analogy with the asymptotically locally flat case, smaller value 
of K are unstable perturbations. As a consequence, black strings are unstable 
if they are longer than the length associated with the critical wavenumber. 



5.3.1 The equations and boundary conditions 

We consider the Einstein equations with cosmological constant supplemented 
by the ansatz (|5.2p where the correction functions are developed as a Series in 
e according to (|5.30|1 . The relevant equations for the stability problem are the 
equations appearing at order e. 

Following [81] , we express Bi as a function of Ai , Ci and their first derivative 
using the equations of motions: 

_ a {2b (r {A[ + {d ~ 3)C[) + {d ~ 3)Ci) + rAjb') - rba'jAi + (d - 3)Ci) 
^ ~ a(r6' + 2(d-3)fe) 

(5.31) 

The resulting equations of motion for the fields Ai , Ci are given by 
A'l = aiAi+a2A[+a3Ci+aiC[, C" = (^Mi + (^2^'i + <^3Ci + (^4C^ , (5.32) 
where prime denotes the derivative with respect to the r coordinate and where 
26 ((-3 + d) k'^e + (r - dr) a') + r {k'^P + 2 (-1 + d) a) b' 



Oil = 



a2 



Paf{2{-3 + d)b + rb') 



1 b' 1 



r 2b Prf 



i-l + d) r"^ + {-A + d)e - 



2 4(-l + d)r^6 



2 (-3 + rf) 5 + rb' 



as 



a4 



2 (-3 + d) (-1 + d)b (2a - ra') 
£2a/ (2 (-3 + d)b + rb') ' 
{-'i + d)b' A{-d + d){-l + d)rb 



2b {2{-3 + d)b + rb')' 
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V2 



^3 



ip4 



2 ((-1 + d) + (-4 + d) P) {-ba' + ab') 
Praf (2 (-3 + d)b + rV) ' 

4&((-l + d) r2 + (-4 + d) £2) 
^ Pf{2{-?, + d)b + rb') 



(5.33) 



1 



- [2b{d - 3) {k'^fr + 2{d - l)ra - {{d - l)r^ 



Praf {2{d-i)b + rb') 
+ (d - 4)^2) a') + f [k^r^ - 2 (-4 + d) a) b'] , 

- £Va/(2(d-3)5 + .50 [^^-3 + ^) (^^^V + 2 (-1 + .) .a 

- ((-1 + d)r^ + (-4 + d) f) a') + £^ [k^r^ - 2 (-4 + d) a) b'] 



This system of equations requires five boundary conditions; the appropriate 
boundary conditions are given by 



^i(c») = Ci(oo) = 0, Ci{rh) = l, 
A'i(r,) = An , C[{rn) = Cii. 



(5.34) 



The first two conditions are required in order to describe asymptotically van- 
ishing perturbations, the third condition is imposed using the linearity of the 
equations and the last two conditions are required for regularity of the equa- 
tions at the horizon, with An, Cn given later by (|5.39p . 

Note that we have 5 boundary conditions for a system of total order 4; fc^ 
will have to assume very specific values in order that all the boundary condi- 
tions are fulfilled. The system of equations is in fact similar to an eigenvalue 
problem given by a system of two second order coupled differential equations, 
the eigenvalue being the square critical wavenumber fc^. The technical way to 
tackle this problem is to supplement the system of equations with the trivial 
equation 



dr 







(5.35) 



and to hope that the numerical algorithm will be able to find the appropriate 
solution and the corresponding value of k^ . It turns out to be the case. 

Let us also stress that the system of equations still admits the same scaling 
relations that the uniform string one (|5.27p , supplemented by the scaling of the 
critical wavenumber: 



k = fc/A. 



(5.36) 
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5.3.2 Near horizon and asymptotic expansion 

In the limit where r is close to rh, one can perform a Taylor series expansion 
of the solution according to 

Ai{r) = Aio + An{r-rh)+0{r~rhf, 

Ci{r) - Cio + Cn{r-rh)+0{r-rhf, (5.37) 



where 



, _ 2a. (Aio-Cio)(d-4) id-3)f 

- 3ahr,Aid~4)P + id-l)rl) ^^'^^^ 
(2aUio id - 5) (d - 1) + (-2Aio + Cip (d - 3)) Pi^) rl 
^ 3ahrhiid-4)e + id-l)rl) 

_ 2(Aio-Cio) Cio(2a,,(rf-l) + fc'^')^/. 
" " +a,(d-4)£2 + a,(d-l)r2' ^"'^^^ 



and ah was defined previously, Aio, Cio are real constants. Note that due to the 
linearity of the first order equation, either Aiq or Cip can be fixed arbitrarily, 
only the ratio Aiq/Ciq is non arbitrary. We used this freedom to impose 
Ci(r/i) — Cio = 1 in the boundary conditions 

In the asymptotic region, where r — !• oo, the fields Ai{r) and Ci(r) decay 
according to 

d-l d+2 

r 



Ai(r) = _(d-3)7i(^-j +0 

Ciir) = jJ-) +o(-)'' , (5.40) 



where 71 is a real constant to be determined numerically. The influence of the 
logarithmic terms appearing in the background expansion for odd values of d 
shows up in higher order terms in (|5.40p . In fact, this is not the most general 
asymptotic expansion for the first order fields. The most general expansion is 
of the form 

Mr) = 7o(l + ^j-(d-3)7i(-j +0 
Cir) = 72(l + ^)+7i(^)' (5.41) 



for real values of 70 , Ci , C2 • These additional terms all cancel once the bound- 
ary conditions (j5.34ll are imposed. 
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5.3.3 Numerical results and the Gubser-Mitra conjecture 



We used the ODE solver Colsys in order to integrate the equations (j5.32p . 
Once given a starting profile sufficiently close to the solution, the solver was 
able to produce an accurate solution for both the functions Ai , Ci and the 
corresponding eigenvalue fc^ for a given value of the number of dimension, the 
horizon radius and of the cosmological constant. 

A positive value of k'^ indicates the presence of an unstable mode of the 
uniform black string, as in the asymptotically locally flat case. A negative 
value of fc^ leads to the absence of an instability, since in this case, the critical 
wavelength is imaginary. 

It should be stressed that the thermodynamical quantities are not affected 
by the perturbations. Indeed the perturbations are linear in cos kz and the 
thermodynamical quantities involve an integration over z from to L = 27r/fc. 

As mentioned before, we worked at fixed horizon radius and rescaled the 
solutions (in particular the critical wavenumber) in order to bring the stability 
analysis in a framework where the AdS radius is fixed. Our results suggest 
that the Gregory-Laflanime instability persists in asymptotically locally AdS 
spacetimes, at least when the ratio rii/£ is small (compared to one). The 
dynamically unstable uniform black strings (with rh/£ ^ 1) are referred to 
as small black strings while the dynamically stable black strings are referred 
to big black strings. The same phenomenon occurs for AdS black holes, small 
AdS black holes are unstable while big AdS black holes are stable. As stated in 
[88] , the AdS radius acts like a confining box; if the horizon radius is small, the 
black objects doesn't 'feel' the curvature, while for larger values of the horizon 
radius, the black objects feel the 'box' and are stabilized. 

Note that the dynamical instability matches the thermodynamical insta- 
bility, in accordance with the Gubser-Mitra conjecture [23] , stating essentially 
that there is a correlation between thermodynamically stability and dynamical 
stability. This is shown in figure [5751 Note that the Gubser-Mitra conjecture is 
not always satisfied |110) . but seems to hold for simpler models where matter 
fields are absent. 
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Figure 5.5: Plots of the entropy and as functions of Th for d = 5 and 
A = — 1. The red hne corresponds to the temperature where the thermody- 
namically stable phase > 0) meets the thermodynamically unstable phase 
(resp. < 0). This temperature is precisely the temperature separating the 
dynamically stable (A:^ < 0) and unstable (fc^ > 0) phase. This confirms the 
Gubser-Mitra conjecture in the case of AdS black strings. 




1.74 1.75 1.76 1.77 1.78 1.79 1.E 

T„ 

Figure 5.6: The same figure as 15.51 for d = 8. 
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In asymptotically locally flat spacetimes, we can consider that there is only 
one unstable black string, since the only dimensionless parameter is rh/L. So 
one can fix arbitrary and compute the critical wavenumber, fixing the critical 
length Lc = 27r/fcc- However, in asymptotically locally AdS spacetimes, we 
have a new dimensionful parameter, namely the AdS radius £. So now, we can 
construct two dimensionless quantities: the rescaled horizon radius rh/i and the 
rescaled length L/£. In other words, for each values of rh under consideration, 
we will find a critical length for a given value of £ (figure [577]) . This of course 
complicates the picture since now we find a continuum of unstable uniform black 
strings characterized by G [0, r^], depending on the number of dimension. 
This critical value defines the boundary between small uniform black strings 
and large uniform black strings. Note also that for every number of dimensions 
we considered, the value of increases when — > 0. This implies that the 
globally regular solution also presents a Gregory-Laflamme instability. 
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Figure 5.7: The critical wavenumber as a function of the horizon radius for 
different values of d and for A = —1. 

In analogy with the flat case, we expect a non uniform string phase to 
emerge from the unstable black strings. In the next sections, we will present 
numerical arguments for the existence of non uniform AdS black strings. We 
use first a perturbative approach, following the lines of Gubser |81| . then we 
perform a non-perturbative approach, solving the appropriate system of cou- 
pled non linear partial differential equations. 

The non uniform string phase has A,B,C ^ and emerges from the un- 
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stable black strings. There are two approaches in order to construct the non 
uniform branch. The first involves a perturbative approach, following the lines 
introduced by Gubser [ST] for asymptotically flat black strings and leads to 
a system of ordinary differential equations. In this case, however, only an 
approximated solution can be produced. 

The second approach consists in solving the non linear partial differential 
equation problem (see [55] for the asymptotically locally flat non uniform black 
strings). We will present the perturbative analysis in the next section and the 
solution to the system of non linear partial differential equations, using results 
based on the perturbative results in section 15.51 

5.4 Perturbative non uniform strings 

In the perturbative approach, we develop the functions A, B,C as a Fourier 
series and in term of a small parameter e |81j : 

A = eAi(r)cosfcz + (Ao(r) + A2(r) cos2fcz) + 0(e)^, 

B = eBi{r)coskz + e^{Ba{r)+B2{r)cos2kz) + 0{ef, (5.42) 

C = eCi{r)coskz + e^{Ca{r)+C2{r)cos2kz)+0{ef, 

where k is the critical wavenumber. The perturbative equations are given by 
the Einstein equations developed as a series in e. 

As already noted, order e*^ is the uniform black string problem, while order 
is the stability problem treated in section [5731 Order involves the back- 
reacting modes Xq of the metric, induced by the first order perturbation and 
the massive modes X2, X denoting generically A, B, C . 

The backreacting and massive modes are decoupled to the order since 
they appear linearly to this order and with linearly independent Fourier modes. 
In this thesis, we will consider only the backreacting modes since they are the 
only modes that contributes to the thermodynamical corrections at order e^, 
as we will see later. 

5.4.1 Equations of the backreacting modes 

Like for the first order equations, it is still possible to solve algebraically Bq 
and to reduce the problem to a system of two second order coupled differential 
equations. However, the expression for Bq is quite long but is straightforward 
to compute. We will not give the expression for Bq but only mention that it is 
of the form 

So = L(Ao,Co) + Q(Ai,Ci), (5.43) 
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where L is a linear diflFerential operator and Q is a quadratic expression involving 
Ai , Ci and their derivatives, with coefficients depending on the background 
functions, on the number of dimensions and on the AdS radius. 
The resulting equations are of the form 

Aq = ixoAq + oiCq + terms quadratic in Ai,Ci, A'l, C[, 

Co = 7oCo + 71^0 + terms quadratic in Ai,Ci,A[,C[, (5.44) 

where 

ao = -|2r-&6' (^(d-2)(2d-5)/+^^^— ^^^-^i^ + (d-4)(d-2)^ 
+ (d - 2yfb" + 2b' [{d - 3)(d - 2)\f + (d - 4) ^ (d~2){d~iy 
+ (d - 3)(d - 2)] ] I / {2{d - 2)rbf {rb' + 2{d - 3)6) } 
ai = [(d-3)(-4(d-3)r6/6'-rV6'^ + 26^(-(d-4)(d-3)/ 

/ [2rbf {rb' + 2{d - 3)b)] 

70 = {-r66'(^-(d-2)(2d-9)/+^^— ^^g-^i^ + (d-4)(d-2)^ 

+ {d - 2)r'fb'' + b' (^{d - 8)(d ~ 3)(d - 2)/ - (d - 4) (^(d-2)id-iy 
+ (d - 3)(d - 2))) |/|(d - 2)rbf [rb' + 2(d - 3)b) } 

b ((d - 3)(d - 2)/ + + (d _ 4)(d - 3)) 

" (d-3)r/ (rfc' + 2(d - 3)fe) ' 

The quadratic terms in (|5.44p are really long and we prefer not writing them 
explicitly. They are anyway straightforward to compute. 

We note that only the derivatives of Aq, Cq appear in equations (|5.44p . As 
a consequence, the equations are invariant under 

Ao ^ Ao + Const. , Co ^ Co + Const. (5.46) 

This will be useful to integrate the equations since it allows to shift the solutions 
down to zero asymptotically and thus not to care whether the functions Aq, Cq 
go to zero at infinity as long as they go to constants. 

5.4.2 Near horizon and asymptotic solution 

The second order fields have the following near horizon expansion: 

Ao{r) = AQQ+AQi{r-rh)+0{r-rh)'^, 

Coir) = CQQ + CQi{r^rh)+0{r-rhf, (5.47) 



no 
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where Aqo , Coo are to be fixed using the invariance ()5.46|) such that the fields 
Aq, Co decay to at infinity, Coi is an arbitrary real constant and 

^ _ ahC'oi id - 4) (d - 3) £2 + (Alo + 2AwCio (d - 3) - C% (d ~ 3)) kH^rh 
Aqi — 



anCoi [d - 3) (d - 1) rl 



^au [d-A)P + Zah {d-l)rl 



(5.48) 



This expression for Aqi is obtained by imposing regularity at the horizon; if A^i 
is not chosen to be the above expression, there will exist terms proportional 
to [r — rh)~^ in the near the horizon development. Such terms are clearly 
diverging near the horizon. 

The backreacting fields A^ , Cq follow the same asymptotic form as the first 
order fields to the leading order: 

Ao(r) = -(d- 3)70 (J)' \o[^-^'^\ 

Coir) = 70 (^)' (5.49) 

where 70 is a real constant, also to be determined numerically. Once again, 
the most general expansion involves lower order terms which would give infi- 
nite contribution to the mass and thus are non- physical [21]. We will impose 
boundary conditions in such a way that the non physical terms vanish. 

5.4.3 Boundary conditions and numerical technique 

The numerical technique is the following; we first integrate the background 
fields then we integrate the first order fields with the solver Colsys. The back- 
reacting fields are then integrated with the following initial conditions: 

Aoirh) = ao, C(ni) = Xo, C"(r,i) = xi 

A'oirh) = -[ahXiid-4)id^3)e^ + {Ala + 2AioCioid-3)^ 

C\^ id - 3)) + a,,xi id - 3) (rf - 1) rl}/{3ah id - 4) f 

+3ahid-l)rl'j, (5.50) 

where q;o,Xo are arbitrary constant to be fixed a posteriori using the shift 
invariance (|5.46p . while xi is a constant which is tuned such that the fields 
Ao, Co follow the decay ([05]) . 

In practice, we integrate the second order equations using a Runge-Kutta 
algorithm at order 4. The integration is carried out from the horizon radius, 
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to some R'> Th. The fields Ao,Co follow the asymptotic ()5.49p if RC'^{R) + 
{d — l)Co{R) = 0. The problem here is that Co can always be shifted by an 
arbitrary constant so we impose the decay (|5.49p transposed to the derivative 
of Co: 

RCq{R) + dC'a{R) = 0, (5.51) 

the value of Cq{R) being obtained using the field equations. 

The value of RCq{R) + dC'Q{R) is a function of xi, say C(xi)- The backre- 
acting fields will follow the asymptotic decay (|5.49p if xi is chosen to be a root 
of C{xi). We used a Newton algorithm in order to find the value of xi such 
that C(xi) = 0. 



5.4.4 Thermodynamical corrections 

As mentioned before, the first corrections to the thermodynamical quantities 
arise from the backreacting modes, at order 2 in e. This is due to the fact 
that there is an integration over z € [0,L] which discards contributions from 
massive fields at this order in e. Of course, the corrections at order depends 
on the first order fields since they come squared in the development of the 
metric component. 

The correction at order on the entropy is given by |25) 

ah {-Alo + 2AioCio(d - 3) + (4Coo + C^d - 3)) {d - 3)) {d - 1) 

Oj/bo = 



4a,,(d- 1) 
CMAw+Cw)id-3)ee 
Aahid~l) 



(5.52) 



where the background quantities are evaluated with the length L = 27r/fc and 
So is the background entropy. The total entropy is given by 5*0 + e^SS. 

The correction to the temperature at the horizon is obtained by demanding 
regularity in the Euclidean sections and developing the condition as a series in 
e. Here again, the first correction arise at order and is given by: 

,^ 2a,(d - 1) (2^00 + Afo) +{d- 3)(Aio - Cio)Ciofc^£^ 

^^^/T"- i[^;;id^) ' (5-^3) 

where Th is the background temperature and STh the correction from the 
backreaction, the temperature of the non uniform phase being denoted by Th + 
e'STH. 

The mass and tension can be computed using the counterterm procedure 
|107| defined in section 15.2.31 Note that an application of the counterterm 
procedure to perturbative non uniform black strings in 6 dimensions can be 
found in [Ml. 
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However, the thermodynamical quantities involved in these procedures are 
difficult to extract from the numerical solution in the general c?-dimensional 
case because of the backreacting fields asymptotical decay (|5.49p . Instead, it is 
in principle possible to compute the mass by integrating the first law with fixed 
asymptotical length (L fixed) and to use the Smarr formula [20] to extract the 
tension. However, we will postpone discussion of the mass and tension to the 
next section where the full non perturbative solution will be presented. Instead, 
we will investigate the thermodynamical stability properties of the perturbative 
non uniform black string. 

In order to investigate the thermodynamical stability of the non uniform 
phase, we consider the specific heat, defined as the derivative of the entropy 
with respect to the temperature. In the case of perturbative non uniform black 
strings at order e^, the specific heat reduces to 

C^^=T^^, (5.54) 

with L fixed, 6S and 6Th as defined above. 

This quantity gives the variation of entropy with respect to the Hawking 
temperature in the non uniform phase at the emanating point (where the non 
uniform phase emerges from the uniform phase); non uniform solutions with 
Cp < Q are thermodynamically unstable, while solutions with Cp > are 
thermodynamically stable. Since the temperature Th of the background is 
always positive, it is sufficient to compare the sign of SS and STh in order to 
investigate the existence of a stable phase. 



5.4.5 Numerical results and discussion 

We integrated the field equations for values of d from 5 to 15. We considered 
perturbative non uniform black strings emanating from the uniform phase with 
length fixed by the critical length L = 27r/fc found in the stability analysis. 
Our results show that for every number of dimensions considered, there exists 
a phase of non uniform black string presenting a positive specific heat (figure 

This stable phase exists for small values of the critical wave number, i.e. 
large value of the length (with £ fixed). This is the analogue of the small/big 
black strings or black hole in AdS where the relevant parameter for the ther- 
modynamical stability was rh/£. Here, the relevant parameter is fii = L/£; 
short non uniform black strings solutions with /ii <C 1 are thermodynamically 
unstable while long non uniform black strings with /ii w 1 are thermodynami- 
cally stable (see table lS.ll where we present the value of /ii separating the stable 
phase from the unstable phase, say nf). 
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Figure 5.8: The correction on the entropy 5S as a function of the correction 
on the Hawking temperature 5Th for various number of dimensions. There are 
points on this figure with STh > 0,(55 > corresponding to a thermodynam- 
icaUy stable phase. The box in the upper left corner is a zoom of the region 
with 6Th > 0. 



d 


6 


7 8 9 


10 


11 


12 


Ml 


8.66409 


8.76464 9.02511 9.3052 


9.56981 


9.80617 


10.1825 



Table 5.1: The value of the ratio ~ L/£ where the stable phase occurs for 
various d. 



Let us stress the fact that these short/long black strings are present in 
the small AdS phase, since it is the phase we are dealing with. We expect 
this property to be a generic feature of non uniform black strings in arbitrary 
number of dimensions. 

Figure 15.91 shows the direction of the phase transition in a 5* — Th diagram 
for d = 9. In the short black string phase (i.e. with L ^ £), the entropy 
increases and the temperature decreases along the non uniform phase while 
in the long black string phase (L « £) the entropy increases and so does the 
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temperature. Let us emphasize once again that this new phase of non uniform 
black strings occurs for small black strings (ro ^ £) since the length is not well 
defined for big black strings (rg ~ £). 




1 1 1 1 

0.7 0.8 0.9 1.0 1.1 1.2 



Figure 5.9: The phase diagram in the S — Th plane for d = 9. The background 
and the direction of the corrections for some value of the rescaled length are 
shown. The same pattern holds for all the number of dimension we considered. 



5.5 Non perturbative solutions 

In this section, we present the construction of the non perturbative solutions. 
The problem is given by a system of three coupled partial differential equations 
and a set of boundary conditions. 

5.5.1 The coordinate system 

The solutions are a priori independent of the normalisation of the functions a, b; 
this means that the functions a, b determined by (|5.7p as boundary conditions 
have to be multiplied by a suitable factor in order to follow the asymptotic 
behaviour (|5.3p a posteriori. This invariance is related to the arbitrary rescaling 
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of z and t. In order to avoid this scaling ambiguity, we perform the following 
change of functions 

2 2 2 

air) = ^a(r), b{r) = ^6(r), f{r) = ^j^f{r). (5.55) 

The new functions a,b,f then obey a{r — > oo) = l,b{r — > oo) = 1 which are 
imposed as boundary conditions independently of £ and d. Consequently, the 
functions reconstructed from d and b are normalized properly. 

Note that the z coordinate ranges from to some length L. The value of L 
is fixed to be the critical length where the linear perturbations are static. This 
precise value of L is again provided by the stability analysis coming from the 
first order in perturbation theory |106j , see section 15.31 

It must be stressed that the equations resulting from the ansatz (j5.2p re- 
quire complicated regularity conditions for the functions A, B, C at the horizon. 
Going in a 'conformal-like' gauge, 

+g(f)e2^('--^)df7L3, (5.56) 

where g{f) = f^^+r^, with the horizon radius, leads to much simpler regular- 
ity conditions. Note that we rescaled the z coordinate such that z G [0, 1], z = 
z/L hy factorising the length L. 

The relation between (|5.2p with the change of functions (|5.55p and (|5.56p 
is given by 

f2 + rl = r\ ~a{f) = a(r), &(r) = 6(r), />) = /(r). (5.57) 

Note that a subtlety arises when one wants to consider a conformal gauge, 
starting from the line element (j5.5p . Naively, one would write 

g{r)=f^+rl, (5.58) 

following the modification from (|5.5p to (|5.56p . This doesn't work at all; the 
near horizon expansion suffers from divergences of order {r — rh)~'^- The correct 
way to enter the conformal gauge is to replace the diverging r^/f^ by g{f)/£'^ 
since the change of gauge is such that — > g{f)- To the leading order, it 
doesn't change anything asymptotically, but it appears that the near horizon 
expansion doesn't suffer any more from unmanageable divergences. So after 
some algebra, the conformal gauge with factorized AdS asymptotic is given by 
(|536| . 

From now on, we omit the tilde over the r and z coordinates and the hat 
over the various functions, except if this spoils clarity. 
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5.5.2 The equations and boundary conditions 

The equations are the Einstein equations with the ansatz (|5.56p , leading to the 
following set of coupled partial differential equations: 

(-l + d)e2Sr2\ b' a'h' h'^ b' f g' a'g' {2 + d)b'g' 



+ 



fg^ J 2rb Aab 46^ 46/ 2rg Aag Abg 

fg' i-l + d) g'^ b^ g^ ^^^^^("'i)' {-3 + d) ^V^(04)c(°-i) 
4/^ ^ V + 26 ^ 2^ ^ L^afg^ ^ L^afg^ 
^V^(0,2) ^(1,0) ^,^(1,0) ^,^(1,0) f>^{l,Q) (2 + d)g'A(i>o) 



L^afg^ r 2a b 2f 2g 

26 2g 



i-i + d) (-3 + d) e^^-2g^V (-4 + rf)(-l + d)e^^r^ a' g'^ 
2/^ ^ 27^^2 2^ ^ 4^ 

_^«T _ ^ ^ ^ _ (-3 + rf)&V _^ /V _ (-4 + rf)(-l + rf)g-^ 
4a/ 2r.9 ag Abg Afg S.g^ 

a" g" (-3 + d)^Vyl(0'i)C(0'i) (-4 + rf) (-3 + d) 



2a 2g L^afg^ ^L^afg^ 

^4^2^(0,2) (_3 + rf)g^^(l,0) ^(1,0) g'gd.O) /^(l-O) 3ff^5(l-°) 

L2Q/g3 2.9 r ^ 2a 2/ ^ 2g 

{-3 + d)b'C(^'"-> M + d)Vc(i-") _ ^ ^(1.0)^(1,0) 

^(-4 + .)(-3 + d)Ca^o)-^^,,,^^ ^^^^^^ 



_ / (-4 + rf) e^^-^'='<^V^ \ _ j-l + d) e^^r^ ,g' a'ff' 6'ff' fg' 
V 7? / 7? 2^ ^ 4^ ^ 46^ ^ 4/^ 

+ g'' ^4^2^(0,1)^(04) (_3 + rf) £Vc(°-i)' £^^^(^(0,2) 

V ^2^^ L^afg^ ^ L^a/gS ^ L^afg^ 

/^(l.O) ^(1,0) „,^(l,0) 5^(^(1.0) j,^(l,0) (_i + rf)g'c(i.o) 

1 1 h h -^^ — 

2g r 2a 26 2/ g 

^^(1,0)^(1,0) ^ (_3 ^ ^) ^(1,0)2 ^ ^(2,0) ^ Q (5 g^) 
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where the exponent (m, n) denotes the m*'' derivative with respect to r and n*^ 
derivative with respect to z and the primes denote the derivative with respect 
to r for functions depending on r only. 

It should be noted that the terms containing second derivatives of the func- 
tions A, B, C are all of the form 

X denoting A, B and C. 

The combination of the Einstein equations leading to equations (|5.59p . 
([5:TO)1 and (i53T|l are El +Ei:, El and E^ where El = Gl- Ml and 9 denotes 
the angular sector. Note that there are two other linearly independent non van- 
ishing combinations of the Einstein equations; these are constraints: El — E'!^ 
and £'*. It has been argued in [51| that the constraints satisfy the Cauchy- 
Riemann equations, as a consequence, if they are fulfilled on the boundary of 
the domain of integration, they will be fulfilled in the interior of the domain as 
well. Here, although there is a cosmological constant, the argument is the same 
since the cosmological constant cancels for the particular combination leading 
to the constraints. 

The system of partial differential equations above is supplemented by the 
following boundary conditions: 

A{r, z) = B{r, z) = C(r, z) = for r oo, 
dzA{r,z) = dzB{r,z)^dzC{r,z)={)iov z = Q,l, (5.63) 
drA{r, z) = drC{r, z) = 0, B(r, z) = A(r, z) + do for r 0, 

where do is a real parameter. The first condition ensures the fact that A, B, C 
are corrections and thus vanish at infinity, the second conditions imposes peri- 
odicity along the z coordinate while the last three conditions are the regularity 
conditions and ensure that the constraints are fulfilled. The parameter do is 
related to the variation of the temperature along the non uniform branch: the 
variation of temperature along the non uniform branch is a function of to B — A 
evaluated at the horizon (i.e. of do)] we will come back to the relation between 
do and the temperature later. Let us stress anyway that do plays a crucial role 
in the construction of non uniform solutions: it is the parameter which forces 
the solution to deviate from the uniform case; if do = 0, the solution to the 
system of partial differential equations is the trivial solution A = B = C = 0. 

The domain of integration is z e [0,1], r G [0, oo] (recall that r is the 
radial coordinate in the conformal gauge). We used the solver Fidisol (see 
appendix [E| in order to solve this system. Fidisol is sensitive to the mesh 
and boundary conditions (it needs exact boundary conditions, i.e. imposed at 
infinity in this case); in order to impose appropriately the boundary conditions 
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(|5.63p especially the one at r = cx), we have chosen to work with a compactified 
variable x ranging from to 1: 

X = "—II^, (5.64) 

r 

We have also tried different meshes in the x direction, as well as in the z 
direction in order to test the robustness of the solution. 

Fidisol is based on a Newton-Raphson algorithm and needs a 'good initial 
profile' (see appendix [E|to efficiently provide a proper solution. We tried to 
start with the uniform string solution, A = B = C = and to increase the 
value of do but this approach failed to provide convincing solutions. Starting 
from A = eyli(7~) cos27rz, B = eBi{r) cos27tz, C = eC'i(f) cos 27rz. where 
Ai{r) — ^i(-\/r^-»^) and similarly for Bi,Ci { Ai, Bi, Ci being the solution 
from first order in perturbation theory, section 15. 3p , for a small real value of 
e and a small value of do works much better and provides robust solutions. 
Then, we increased the value of do starting with the previous solution as an 
initial profile. 



5.5.3 Asymptotic behaviour 

The asymptotic behaviour of the background functions was reminded in the 
first sections of this chapter; since we factorised the divergent term, the 
asymptotic development in [5D] is of course to be divided by g{r)/£'^, leading 
to 

1^1 ,4-1 "i-l 'i-l 

L 2 J /^2y+i F g fe 



+0 



\9 



d + 1 

2 \ — 



1+ E «Mt +^*hr +E'^'i^2fe+4iog^ - 



I I / ^a,zh,-^i n /J9 I 



d + 1 

2 \ — 



+0 



(f) 



d-1 



/ = 1 + E ( ~ ) + (ct + + Co) ( - ) + 



£2 \ »+l / p2\'-2 



d+1 

l?2\ — 



Ew4log|(^) , (5.65) 
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where g = r'^ + 



Note that 1/g in the above formula could be further 



expanded, leading to a much longer expression. 

It has been argued in |24[[^ that the modes of A, B, C in the perturbative 

approach decay as (I)'' ^ in all order of the perturbations. Since the non 
perturbative solution is a combination of these modes, the leading order should 
also decay as 



A(r,z) 
C(r,z) 



ai - 



7i - 



d-1 



r 



oi- 

r 



d+l 



, B{r,z) 



o - 



d+l 



d+l 



(5.66) 



for some ai , /3i , 71 G K. The decay (|5.66p is confirmed by direct computations 
on the numerical solution. 

We were able to obtain solutions to the equations (|5.59p . (|5.60p . (|5.6ip 
with the boundary conditions (j5.63p for various values of d, rh/i, do. For a 
given value of the number of dimensions, the non uniform solutions are given 
as families of solutions, labelled by rh/i and do. The value of rn/f- fixes the 
critical length L of the corresponding uniform string and da allows to slide 
along the non uniform branch for the given value r^/i- Before presenting 
the thermodynamical properties of the non uniform phase, we present generic 
profiles of the functions A, B, C in figure [5TTU1 





1 



d=7 , =1, A=-0.1 , da =0.0175 

Figure 5.10: Typical profiles of the correction functions A, B, C for d = 5, = 
1, A = -0.1, da = 0.00375 and d = 7, r,, = 1, A = -0.1, do = 0.0175. 
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5.5.4 Properties of the solution 
Thermo dynamical properties 

The thermodynamical properties of the background have been presented above. 
The correction functions A, B, C provide corrections on the thermodynamical 
quantities computed from the background (as the name 'correction function' 
suggests). The quantities computed on the background solution are denoted 
with a subscript J7, while the quantities referring to the full non uniform solu- 
tion are denoted with a subscript NU . 

The corrections on the thermodynamical quantities are computed using the 
standard procedure (regularity in the euclidean section and one quarter of the 
horizon area) and are given by 

Tnu - ^ V&'(nO/'(nOe^('--^)-^('---) = e-'''>Tu, 
47r 



L Jo 

where Tnu ^nd Snu a-re the Hawking temperature (resp. the entropy) of the 
non uniform string. Note that the first relation makes the role of do clearer: 
once the background is set, the parameter do allows to slide along the non 
uniform branch in a temperature-entropy phase diagram. 

These quantities are defined at the horizon and are relatively easy to com- 
pute from a technical point of view. The asymptotic quantities can be parametrized 
according to 

Mnu = + 

Tnu = Tu(^l + ^^y (5.68) 

where SM, ST are the corrections to the mass and tension due to the non 
uniformity and are functions of the parameters ai,/3i,7i appearing in (j5.66p 
and of ct,Cz,co (in (|5.65p '). 

Note that at the boundary spacetime, dz and dt are both Killing vectors, so 
the definition of the mass and tension as functions of the asymptotic parameters 
|20j remains the same. Note also that the decay (j5.66p combined with the 
asymptotic form of the background solution gives contribution of order r~^'^~^^ 
to guidzz and leads to correction on ct,Cz. As a consequence, the mass and 
tension are also modified; it follows also that another decay of the corrections 
would then be physically non acceptable. 

The quantities Ct,Cz and worse, ai,/3i,7i are difficult to compute, first 
because of the sharpness of the decay and second because the determination of 
these quantities is plagued by numerical noise. 
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In order to avoid this problem, we were able to integrate the first law of 
black holes and then to compute the masses (of the uniform solution and of the 
non uniform solution). Assuming that the Smarr relation (|5.19p holds for the 
non uniform solution, we were able to compute the corrections on the tension 
for the non uniform solution. 

The procedure to integrate the first law is performed in two steps, first at 
the background level, then for the non uniform solution. The procedure for the 
background level has been described in section 15.2.41 

The integration of the first law for the non uniform solution is more straight- 
forward than for the background itself since we the functions A, B, C are cor- 
rections to the background and indeed provide corrections on the thermody- 
namical quantities, which naturally appear as functions of the parameter do. 
These thermodynamical corrections are related to each other by the first law of 
black holes. It follows that the procedure to integrate the first law seems more 
natural in this case; given a uniform black string with horizon radius r/j, the 
non uniform black string with a given value of dg emanating from this uniform 
black string has a mass given by 

MMuirh) = A%(r,)+ r"T#^(r,„do)^^ddo 

= AIu{rh) + TS{rH) e-'^^ddo, (5.69) 

where SS/Sdo has to be evaluated numerically. Terms of the form TdL don't 
appear since the length of the non uniform black string is fixed along the non 
uniform branch. 

Once the mass of the non uniform string for a given value of do is computed, 
it is straightforward to obtain the tension for the corresponding value of do- 
using the Smarr relation; 

r^^ ^ _ (5.70) 

Geometric properties 

The order parameter for the transition from the uniform string phase to the 
non uniform string is defined by |81j 

X=^l^-1, (5.71) 

where Rmui (resp. Rnax) is the largest (resp. smallest) value of the areal 
horizon radius, defined as R = r^e'-^^'''"^^ in Schwarzschild-like coordinates. 

The parameter A measures the deformation of the solution; the uniform 
string is characterized by A = 0, while the hypothetical localized black hole 
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phase has an infinite deformation parameter. In that sense, it measures how 
'far' the non uniform solution is from the uniform solution; i.e. it characterizes 
the non uniformity. This is the reason why it is used as the order parameter in 
many references dealing with phases of black strings (Wiseman (55] , Wiseman- 
Kudoh Uni], Gubser EI] etc). 

Relation between deformation and thermodynamics 

It has been argued in |109j that using the mass and tension is more relevant 
in the construction of the phase diagram than another choice of parameters 
(including A); for instance, in the asymptotically locally flat case, the mass and 
tension remain finite in the three phases (uniform, non uniform and localised 
black hole/string), while the entropy diverges for small temperature. However, 
the mass and tension are defined at the boundary spacetime and are difficult 
to compute from a the numerical data, as already stated. Moreover, in the case 
of asymptotically locally AdS strings, the fact that the background tension is 
negative in some region of the parameters makes the picture more difficult to 
interpret. 

An alternative is to study the phase diagram in the Th, S plane; these 
quantities are defined at the horizon of the black objects, and are easier to 
compute. Moreover, the entropy does not diverge in the AdS case since there 
exists a minimum temperature for the uniform strings where the stable branch 
connects the unstable one. 

Two remarks are in order concerning the deformation as a function of the 
parameter do, as pointed in |25| : 

• A increases with do, 

• the slope of A as a function of do around do = is larger for larger absolute 
values of the cosmological constant. 

These remarks can be formulated by the following statement: 

The deformation is related to the variation of entropy and do is related to 
the variation of the temperature. Moreover the rate of variation of A with do 
is proportional to the variation of the entropy with respect to the temperature 
along the non uniform phase. 

The above statement can be argued in three steps: 

1. For small values of A, the non uniform string is well described by the first 
order of perturbation theory (where C(r, z) « eCi(r) cos27rz, e being the 
small parameter of the perturbation). In this approximation, one can fix 
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the value of e for a given value of A: 



From the boundary conditions of the first order in perturbation theory, 
Ci{rh) ~ 1, leading to e = A/2. 



2. The variation of the temperature along the non uniform branch is given 
by T^^ = Tije~'^°; a variation 5do in do gives rise to a variation of the 
temperature 5T^^ in the non uniform phase according to 

ST^v = -T^e-'^°5do. (5.73) 



Moreover, in the perturbative approach, the correction to the entropy 
appears in order e^; it is then proportional to A^ from the first point: 

,55^^ « A2. (5.74) 

In every cases we considered, the correction on the entropy was always 
positive. 



3. So, increasing do from to do <C 1 implies a variation of A of to A ^ 1 
and 

In other words, the value of A corresponding to a (small) value of do is a 
monotonic negative function of |^jvir, related to C^^ defined previously: 



-C^^do, for do < 1. (5.76) 



To sum up, for small values of do, if gj^Nu is large and negative, A will be 
large and positive and conversely. This is illustrated in figure 15.111 for d = 7 
and for two values of A. 
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Figure 5.11: For small values of do, the deformation increases faster for larger 
absolute values of the cosmological constant, this is to be related to the increase 
of 6S^^ /6T with A (large absolute values of A at fixed r/^ correspond to large 
values of r/j at fixed A from the scaling relations (|5.27p . |20|). 



Entropy difference at same mass 

In the case of asymptotically locally flat spacetime, it is possible to evaluate the 
entropy difference at same mass between the non uniform and uniform string 
using the relation given in (|5.79p below (see for example |102( I81|). To this end, 
the dimensionless entropy s, temperature 9, mass fi and tension f as functions 
of the deformation must be introduced: 

^ e -.^ThL, fi:^ f -.^ (5.77) 



in units where G = 1. 

Then, we define the variation of the rescaled quantities under a variation of 
the thermodynamical quantities: 

X T 

(5.78) 
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The entropy difference at same mass (the variation due the non uniformity 
is compensated in the uniform phase by a variation of the length) is the given 
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by ^mm 

Snu _ 
Su 

where cti , a2 are given by |102) 



l + cri\^ + <J2X* + (5.79) 



d-2 1\ , , 

-1 = ^1 - dT^^'i' -2 - -^(^ (^^1 + ^^ij (5-80) 

In order to derive this (|5.79p . the relation between the entropy and the mass 
of the (d ~ 1)-Tangherhni black hole, which is a foliation at constant z of the 
uniform black string, is used at some stage (see |81j). The computation of (T2 
(ui vanishes |81[|102| ) for various values of the number of dimensions permitted 
to establish the existence of a critical dimension above which the order of the 
phase transition between the uniform black string and the non uniform black 
string changes from one to higher |102) . 

However, in the asymptotically locally AdS black string case, the relation 
between the mass and the entropy for the constant z foliation of the uniform 
phase is not clear. So one cannot derive a similar relation and numerical com- 
putations must be used. 

From a technical point of view, the mass and entropy in the uniform phase 
are functions of r^, say M(rii) and S{rfi) while the corrections on the mass 
(resp. on the entropy) along the non uniform string are functions of and do-, 
say SM{rji, do) (resp. 6S{rh, do)), for a fixed value of the AdS length. One can 
compute the variation of the mass along the non uniform branch and find the 
value of r/i, say r^, such that the uniform string has the same mass: 

rl is such that M (r/J + 6M (r^, do) = M (rl). (5.81) 

Once r,* is known, it is trivial to compute the difference between the entropy 
of the non uniform phase emerging from with a given value of do and the 
entropy of the uniform string with horizon radius r^, as a function of do'- 

^ ^ Snu - Su ^ {Sjrh) + SSirh,do)) - Sjrl) 
Su S{rl) 

We will use the relation (|5.82p in order to investigate the order of the phase 
transition. If the entropy of the uniform string is higher than the entropy of the 
non uniform string of same mass in the small deformation regime, the system 
will first choose the uniform phase. At some stage the non uniform string will 
be more entropic (the non uniform black string somehow interpolate between 
a localized black hole, which is more entropic than the uniform string as the 
argument of chapter |4] suggest), and this phase will be entropically preferred; 
the phase is discontinuous in A and the order of the phase transition is 1 |lllj . 
Conversely, if the entropy of the non uniform string is larger at same mass, the 
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transition is continuous and the order is higher than one. In other words if 
A < for small do, the order of the phase transition is one, if A > 0, the order 
is higher. 

We will present the result of this investigation in the next section. 
5.5.5 Results and discussion 

We first considered the model presented above in 6 dimensions in |24| . with- 
out factorising the asymptotic behaviour of the background functions. Many 
difficulties were encountered during the procedure, at various levels, but we 
were able to obtain some preliminary results indicating the existence of the 
non uniform asymptotically locally AdS black string at least in 6 dimensions. 
The results from the non perturbative regime were in agreement with the per- 
turbative approach. 

The existence of non uniform asymptotically locally AdS black strings was 
confirmed in a more careful analysis, using the 'asymptotic factorisation' and for 
an arbitrary number of dimensions in |22| . However, the numerical procedure 
is still plagued by many difficulties: 

1. The background functions are numerical and it is not clear how the nu- 
merical errors in the background propagates in the corrections. However, 
the background solution is nevertheless computed quite accurately, with 
relative errors typically of order 10^^. 

2. If the parameters e (the small parameter in the perturbative approach) 
and do Sive not chosen in a suitable manner, the solver produces a uniform 
solution which seems to be a 'correction' on the uniform string leading 
to a uniform string with a different temperature. Of course, we are not 
interested in this kind of solutions. Moreover, it happened that the so- 
lution is a kind of combination of this uniform correction and some non 
uniform correction, leading to a decreasing of the deformation along the 
non uniform branch. 

3. The region of larger horizon radius (but still in the unstable phase of the 
uniform black string) is even more difficult to investigate: we expect a 
transition from thermodynamically unstable non uniform black strings to 
thermodynamically stable non uniform black strings following the argu- 
ments presented in section [5.41 

These thermodynamically stable (resp. unstable) non uniform black string 
were referred to as small-short (resp. small-long) non uniform black strings 
since they are characterized by a small horizon radius and a short (resp. large) 
critical length (compared to the AdS curvature). The boundary separating 
these two phases, small-short and small-long black string, is characterized by 
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6S/6Th oo. As we argued in section [5.5.41 the deformation would become 
extremely high even for very small values of do (thus of 5Th)- Of course, this 
extreme behaviour is not manageable from a technical point of view and sets 
the limit of our numerical investigation. However, we were able to find some 
hints for the existence of non uniform solutions in a neighbourhood of this 
region, but the corresponding solutions have large numerical errors and are not 
reliable enough to perform quantitative predictions. 

Despite all these difhculties, we were able to obtain accurate solutions in the 
small-short non uniform black string region, with relative errors of the order of 
2-3%. Of course, in the highly deformed region, the solver starts failing giving 
accurate solutions. 



Hint for a localised black hole phase 



For the reasons described above, we restrict the study of the AdS non uniform 
black string to the region of the parameter space where rh <^ i and L £. 

We were able to construct several solutions for various values of d, £ (or 
equivalently r^) and do- We focused on d = 7 and constructed a number of 
solutions for d = 5,8. We typically considered three or four different values 
of the cosmological constant, ranging from —0.05 to —0.3 and rh = 1. Then 
we rescaled the solutions in order to have the results for fixed cosmological 
constant. For every values of the number of dimensions and cosmological con- 
stant, we constructed many solutions with increasing values of do, typically 
ranging from 0.0001 to 0.01. We were able to construct solutions characterized 
by A w 2.5 in most cases and higher for d — 7. It should be mentioned that 
we did not perform a systematic analysis for every number of dimensions; a 
deeper analysis, confirming the existence of the solutions and providing a more 
detailed analysis should be carried in order to complete the picture and to 
confirm the existence of the solutions. 

For all the solutions we found, the deformation is an increasing function of 
do and we systematically found extremely deformed solution, providing strong 
indications for the existence of an AdS localised black hole phase. In figure 
I5.12( we present embeddings of the event horizon in a 3 dimensional euclidean 
space i.e. of the areal radius R{r, z) = re'~^^'^'^^ evaluated at r = rh\ the strings 
is almost uniform for small values of do and looks like a (deformed) sphere 
touching the borders for larger values of do. 
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Figure 5.12: Embedding of the horizon of non uniform black strings with d ~ 7, 
A = —0.25 in a 3 dimensional euclidean space. The equation of the surface is 
given by + = r^e^'-^^'"'"^-'. The horizon looks like a black hole horizon as 
the deformation increases. 



Horizon radius depending critical dimension 



We have evaluated the entropy difference between the non uniform black string 
and the uniform black string with same mass for d ~ 7 and various values of 
Vfi- Above some critical value of the horizon radius, the uniform black string 
is more entropic than the non uniform black string, while for larger values of 
rh, the non uniform black string becomes more entropic. In other words, the 
phase transition can be of first order or of higher order for the same value of 
the number of dimension, depending on the value of the horizon radius. We 
present a plot of the variation of the entropy difference as a function of do for 
d = 7,£ = 1 in figure [5T3l for various values of the horizon radius. Small values 
have a negative slope while larger value acquire a positive slope, changing the 
order of the phase transition. 
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Figure 5.13: The relative entropy difference between the non uniform phase and 
uniform phase at same mass as a function of do for d ^ 7,£ ^ 1 and various 
value of r/i. 

It should be noted that the value of do in figure [^331 is very small, and thus 
the non uniform solutions are weakly deformed. This regime of small deforma- 
tion is the regime where the solutions to the set of partial non linear differential 
equations are well approximated by the perturbative solutions. The variation 
of entropy and mass are computed according to the procedure described in 
section 15.5.41 involving a comparison between the entropy of the background 
and of the non uniform solution. So on the one hand, we have the background 
uniform solution which is computed very accurately (the relative error is of the 
order 10~^) and on the other hand, the non uniform solution is in the best 
controlled region. Moreover, the quantities involved in the entropy difference 
computation are extracted from the horizon data, which is the most reliable 
region. 

Figure [5.131 thus provides strong indications for the fact that the order of 
the phase transition can change for the same number of spacetime dimensions. 
Although we only investigated the case of d = 7, we expect this property to 
hold for all values of the dimension. Said differently, the immber of dimensions 
above which the order of the phase transition changes should depend on the 
value of rh- 

5.6 Concluding remarks 

In this chapter, we studied a part of the phase diagram of the AdS black 
strings. First, after reviewing the uniform string solution and its thermody- 
namical properties, we studied the dynamical stability of the solution. We 



130 



CHAPTER 5. PHASES OF ADS BLACK STRINGS 



found that the uniform string is unstable for small AdS black strings, i.e. 
when the horizon radius is small compared to the AdS radius. This is in per- 
fect correspondence with the thermodynamically stable region, along with the 
Gubser-Mitra conjecture [53] • Inspired by the asymptotically locally flat case, 
we expected non uniform AdS black string to exist. We computed solutions 
for such objects, first in a perturbative treatment up to order 2 then in a non 
perturbative approach. 

The results from the perturbative solutions gave us strong indications for 
the existence of the AdS non uniform black string phase and permitted to 
compute the first corrections on the thermodynamical quantities in the non 
uniform branch. It also allowed to predict a new thermodynamically stable 
phase of non uniform black string, which we referred to as small-long non 
uniform black string. This is purely an effect of the cosmological constant and 
is somehow the counterpart of the big AdS uniform black string, except that 
in this case, L/£ ^ 1 instead of r}i/£ ^ 1. Moreover, we showed that there is 
not only one non uniform black string (in the sense that in the flat case, the 
solution depends only on one dimensionless parameter Vh/i) but a continuity 
of non uniform black strings emanating from the uniform strings with different 
critical length. This is also an effect of the cosmological constant, since the 
AdS radius provides a new length-scale in the theory. 

The non perturbative approach confirmed the existence of the non uniform 
solutions, and gave hints that the counterpart of the critical dimension in the 
flat case |102| should depend on the value of the AdS radius (or equivalently on 
the horizon radius). Moreover, from the embeddings of the horizon in euclidean 
spacetimes, the non uniform phase is likely to connect to an AdS localized black 
hole phase, still to be constructed. 

Note that we revisited some of the thermodynamical properties of the uni- 
form phase, arguing that it has a negative tension for some region of the pa- 
rameters. 

All these results - the new non uniform phase, the negative tension, etc. 
- should have counterparts in a conformal field theory defined on Sd-3 x S'l, 
thanks to the AdS /C FT correspondence, but this aspect hasn't been studied 
in this thesis. 



Chapter 6 

Conclusions and perspectives 



This thesis, entitled 'New features of black strings and branes in higher dimen- 
sional gravity due to a cosmological constant' was essentially divided in two 
parts. The first part dealt with braneworld models while the second part dealt 
with black objects. 

In both cases, we investigated the effect of considering a non vanishing 
cosmological constant in the model. It turned out that a cosmological constant 
induces drastic changes on the pattern solutions or on their properties. In this 
concluding chapter, we will first review our main results, stressing the effect of 
the cosmological constant; then, we will give some possible perspectives to this 
thesis. 

6.1 Braneworld models 

We studied four classes of topological brane models. In these models, the branes 
consist in a soliton extending in the extradimensions. More precisely, they 
are topological defects available in the theories under consideration; namely 
the Einstein abelian Higgs model, the Einstein non abelian Higgs model, the 
Einstein- Yang-Mills-Higgs model and the baby Skyrme model. The geometry 
of these models were chosen to be non factorisable. 

We considered a cosmological constant in the entire spacetime (a bulk cos- 
mological constant) and modelled the effect of a cosmological constant on the 
brane by considering inflating branes. The four dimensional spacetime con- 
tributes to the equations associated to these models only through the Ricci 
scalar associated to the four dimensional subspace, allowing the replacement 
the four dimensional part by any spacetime with constant curvature. Note that 
inflating brane models are also relevant in a cosmological context, as braneworld 
inflationary models. 
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The consequence of the inclusion of the cosmological constant on the brane 
and in the entire spacetime is that the solutions admit essentially three type 
of extradimensional geometry, depending on the sign of the bulk cosmological 
constant: opened, closed or flat extradimensional geometry. All these geome- 
tries have an angular deficit, related to the brane cosmological constant. A 
common feature of the different models is the occurrence of periodic solutions, 
providing a natural compactification of the extradimensions. In the case of 6 
dimensions, we have shown that the vacuum system is integrable by quadrature 
and we were able to express the effective Planck mass in terms of the metric 
functions. 

It should be stressed that the branes in these models have an extension, 
due to the extension of the soliton describing the branes. This is in contrast 
with the Randall-Sundrum model where the branes are infinitely thin, having 
the consequence that the metric function are not smooth (the derivative of the 
metric functions is not continuous). 

We tried to classify the solutions available in the 6 dimensional Einstein 
abelian Higgs model according to the values of the bulk cosmological constant, 
the brane cosmological constant and the self coupling of the Higgs field. The 
pattern of solutions is very rich but quite involved, due to the number of pa- 
rameters in the model. We also constructed solutions presenting a 'mirror' 
symmetry in the framework of non abelian Einstein- Yang-Mills in 7 dimen- 
sions and non abelian Einstein-Goldstone in d dimensions, allowing once again 
periodic solutions, i.e. naturally compactified spacetimes. 

Finally, we turned back to six dimension with the baby Skyrme model; 
we took advantage of a natural definition of a root mean square radius to 
characterise the extension of the brane. The baby Skyrme model is a toy 
model of the Skyrme model, originally designed to describe nucleons. In that 
sense, the 6-dimensional baby Skyrme model can be thought as a toy model 
containing effective hadronic matter, i.e. realistic matter fields. In this case, 
we also systematically found periodic solutions. 

We believe that the presence of periodic solutions (closed geometry) is a 
generic feature of inflating brane models, independently of the matter content, 
as long as the matter fields are localized, which is a necessary condition in order 
to describe localized branes. This idea is comforted by the fact that periodic 
solutions are indeed vacuum solutions; adding localized matter field should not 
influence too much the asymptotic form of the spacetime. 

6.2 Fermion localization on a brane 

We studied a model in 1 + 1 dimensions for the localization of fermionic fields on 
a kink/anti-kink system. Previous works have been made in this direction, for 
example, the coupling of fermions to kink solution available in the 0^ model 
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and to a superposition of the latter has been studied in [71]. However, the 
superposition is an approximate solution only if the kink and anti-kink are 
separated enough. 

We took advantage of the existence of a kink/anti-kink solution in the Sine- 
Gordon model in order to study the bound states of fermions in such a back- 
ground, without approximation. First, we studied the coupling of fermions to 
kink or anti-kink solution in the Sine-Gordon model. The kink solution in the 
Sine-Gordon model is qualitatively the same as in the c/)^ model; as a conse- 
quence, the structure of fermionic bound states in both theories is nearly the 
same. 

More interestingly, there is an exact solution in the Sine-Gordon model 
describing a kink-antikink system. This exact solution is still valid in the limit 
where the kink and antikink are close. We studied the coupling of fermions 
to the kink/ anti-kink solution available in the Sine-Gordon model and found 
that when the kink and anti-kink are too close, only a few fermionic bound 
states exist. When the distance between the kink and anti-kink is increased, 
new bound state emanate from the continuum. 

Moreover, the eigenvalues degenerate by pair in the limit where the sep- 
aration between the kink and anti-kink is large enough. This is because the 
potential seen by the fermion around the kink and around the anti-kink are 
superpartner potentials. If the kink and antikink are not distant enough, the 
supersymmetry is somehow broken. 

Finally, we studied the stability equation of three particular kink models. 
The stability equation took the form of a Poschl- Teller equation, which admits 
exact solutions. 

The kink and/or anti-kink can be thought as modelling brane/anti-brane. 
As a consequence, the 1 -f 1 dimensional model studied here mimics a five 
dimensional brane and/or anti-brane. 

6.3 Higher dimensional black objects 

Black holes in four dimensions are rather well understood. In particular, they 
are unique, in the sense that they belong to a three parameter family of so- 
lutions, the Kerr-Neumann solutions [74]. The horizon topology of the four 
dimensional black holes can only be spherical. However, in higher dimensions, 
there exists black solutions with other horizon topologies. Among the zoo of 
solutions, there are black hole solutions presenting a spherical horizon topology, 
but there are also solutions with cylindrical horizon topology (black strings), 
toroidal horizon topology (black rings jSOl), solutions consisting in a black hole 
surrounded by a black ring (black saturns |112) ). multi-black rings, etc. 

In this thesis, we focused on spherical horizon topology and cylindrical 
horizon topology. 



134 



CHAPTER 6. CONCLUSIONS AND PERSPECTIVES 



6.3.1 Charged-rotating de Sitter black hole 

We constructed the charged-rotating de Sitter black hole. Unlike the four di- 
mensional case, such solutions are not known analytically in higher dimensions. 
The exact black hole solutions available in higher dimensions are the {A)dS 
rotating higher dimensional black hole (Myers and Perry [84] and Gibbons-Lu- 
Page-Pope [85]), the {A)dS Tangherlini solution, the higher dimensional {A)dS 
charged black hole solutions. The missing part was the {A)dS charged rotating 
black hole. Such a solution has been constructed numerically for a negative 
cosmological constant in [57], but the case of a positive cosmological constant 
is more involved, due to the occurrence of a cosmological horizon. 

We were able to integrate the equation of the charged-rotating de Sitter 
higher dimensional black hole numerically, adding an equation expressing the 
fact that the cosmological constant is indeed a constant. This 'trick' allowed 
sufRcient freedom to impose regularity conditions on the cosmological horizon. 
The solution is asymptotically de Sitter, as expected. We focused particularly 
in the case of five dimensions, hoping that it catches the main features of the 
generic solutions. 

It must be stressed that the rotating solutions behaves quite differently in 
d = 5 and d > 5, d being the number of dimensions. In five dimensions, the 
horizon area goes to zero for some maximum angular momentum, whereas in 
six dimensions and higher, there is no upper bound on the angular momentum. 
Although there are differences in the thermodynamical properties, the solution 
behave qualitatively in the same way at least in some region of the parameter 
space. This is why we expect the de Sitter charged-rotating black holes to exist 
for any value of d. 

Note also the existence of two phases for the {A)dS charged rotating black 
hole; this is to be put in relation with the existence of two phases in the AdS 
Tangherlini black hole. This is typically one of the main difference between the 
solutions with and without a cosmological constant. 

6.3.2 Black string with positive cosmological constant 

Considering a (d — 1) dimensional black hole and adding an extra compact 
direction leads to a black string, the horizon topology is then cylindrical. The 
construction of a black string in asymptotically locally flat spacetime is straight- 
forward. 

We tried to construct a generalization of black string solutions to asymptot- 
ically locally de Sitter spacetimes, but it turned out that such solutions do not 
exist. There are actually black string solutions to the Einstein equations with 
a positive cosmological constant, but they are not asymptotically (locally) de 
Sitter. Instead, the asymptotic spacetime is singular and can be reached by a 
massive or massless observer in a finite proper time. Moreover, the singularity 
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is located at a finite proper distance from outside the event horizon. This is a 
naked singularity; there is no second horizon that would 'hide' the latter. 

We interpret this non existence as resulting from the fact that the compact 
extra direction cannot support the negative pressure induced by the positive 
cosmological constant. 

Note that the AdS counterparts of black string solutions indeed exist and 
have been constructed in [50] . We revisited the thermodynamical properties of 
such black strings and as well as their stability and phase diagram. We will 
review the results in the next section. 

6.4 Phases of AdS black strings 

Black strings in asymptotically locally flat spacetimes are unstable, thermody- 
namically as well as dynamically. This was the result of Gregory and Laflamme. 
The existence of a static perturbation (zero mode) suggested a non uniform 
black string phase; i.e. a black object with cylindrical horizon topology but 
which depends non trivially on the compact direction. Note that the existence 
of a zero mode is the signature of a Gregory-Laflamme instability. 

As we said in the previous section, black string solutions have been extended 
to asymptotically locally AdS spacetimes. There exists two phases of AdS black 
string solutions; one is thermodynamically stable (big black string), the other 
is thermodynamically unstable (small black strings). 

We studied the dynamical stability of such objects by looking for a static 
perturbation, following the lines introduced by Gubser |81j . Our results sug- 
gest that there is a well defined Gregory-Laflamme zero mode for the small 
black strings, while there is no static perturbation for the big black strings. 
Moreover, the point where the dynamically stable and unstable branch meet 
coincides perfectly with the point where the thermodynamical phases meet. 
This matching between the thermodynamical and dynamical stability was con- 
jectured by Gubser and Mitra |23| . 

Let us stress that in the asymptotically locally flat case, there is only one 
static mode, in the sense that once one uses rescaled variables, the value of the 
static mode is fixed. In the AdS case, even with rescaled variables, there is 
a continuity of zero modes, labelled by the horizon radius (or equivalently by 
the cosmological constant). This is due to the fact that the AdS radius is an 
additional length scale in the theory. 

We constructed the non uniform phase emanating from the unstable black 
strings. The picture is quite different than in asymptotically fiat space: we 
provided evidences that there should exist a thermodynamically stable phase 
of non uniform black string in AdS, characterized by a long compact direction 
(relative to the cosmological radius). This new phase seems to be the coun- 
terpart of the small/big AdS uniform phase, but where the relevant parameter 
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is the ratio between the length of the compact direction and the AdS radius 
instead of the ratio between the horizon radius and the AdS radius. We called 
these solutions short and long black strings for obvious reasons. Note that this 
phase emanates from the small uniform black strings, since the big uniform 
black string is stable. 

Our results also provide qualitative evidences for the existence of a localized 
black hole phase in AdS. This can be seen by looking at the embeddings of 
the horizon in an euclidean space; the more we slide along the non uniform 
branch, the more the deformation increases. We therefore imagine that the 
deformation can increase until the horizon disconnects and adopts a spherical 
topology. 

Finally, we argued that the order of the phase transition between the uni- 
form and non uniform black string changes depending on the value of the 
horizon radius (for a given value of the cosmological constant). This is also 
in contrast with the asymptotically locally flat black string, where the order 
in the phase transition changes only above a critical value of the number of 
dimensions (sa 13). 

6.5 Outlook and perspectives 

Many extensions of this thesis are possible... For instance, stability of the baby 
Skyrme model is currently under investigation [11]; one can also consider the 
Skyrme model in 7 dimensions. Another possibility is to consider other type 
of matter fields to describe the brane. We may also consider the interaction 
between two solitonic branes in various models. 

Among all possibilities, an interesting outlook for brane models is to couple 
fermionic fields to the solutions presented in this thesis. We already studied a 
toy model in 1 + 1 dimensions; it would be interesting to study the full n + 4 
dimensional problem. 

Concerning the black strings, our results suggest the existence of an AdS 
localized black hole phase. One obvious continuation of this thesis would be to 
construct this phase. This is however technically much more involved than the 
non uniform black string. Another issue of the black strings in AdS is their 
connection with the dual conformal field theory in the AdS/ CFT context. Note 
that the uniform AdS black string solution would provide new background for 
the dual CFT. 

Finally, the techniques developed to construct non uniform black strings 
might be applied to black rings, producing an object which is not uniform 
along the ring. 
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Some topics 
relativity 



in general 



A. 1 Lagrangian formulation of Einstein equations 
and boundary terms 

It is possible to derive the Einstein equations from a variational principle. The 
action giving rise to the Einstein equations (in vacuum) is called the Einstein- 
Hilbert action: 



where G is the Newton constant, M. is the d-dimensional manifold under con- 
sideration, g is the determinant of the metric on the manifold and R is the 
scalar curvature. 

The scalar curvature is constructed from the metric components in the 
following way: first, the Christoffel symbols are defined as 



The indices range from to d — 1. Then, the Riemann tensor is constructed: 
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Finally, the Ricci tensor and Ricci scalar are given by 
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If we want to describe matter fields interacting under gravity, the simplest 
way is to consider an action of the form 

^ = TTTT^ / + (A.5) 

IdttG J 

where Lm is the lagrangian density describing the matter fields. This way to 
couple matter fields to gravity is called 'minimal coupling'. 

Extremising the action (|A.5p with respect to the variation of the metric 
components leads to the Einstein equations. 

In order to obtain the equations, we first use the identity 

SiV^R + WttG^L,,,) = (Rab - \Rgab - SirGTab^ 5g''\ (A.6) 

where Tab = (ffsr ~ ^9abLm^ is the stress-energy tensor. 

Let us compute the variation of each terms in order to set this relation: 

Sg = ggabSg"' => 5^g = -\v~ggab5g''\ (A.7) 
5R = SiRabg""') = g'^'SRab + RabSg^K (A.8) 

But 5Rab is divergence and contribute as boundary terms, so it doesn't enter 
the equations of motions: 

5Rab - V, {5Vl, - . (A.9) 

Putting everything together, we find 

Gab = SnGTab, (A.IO) 

where Gab = Rab ~~ ^Rgab is the Einstein tensor. 

As we have seen, the variational principle gives rise to a divergence term, 
leading to a boundary term (jA.9p . In order to take this divergence term into 
account and obtain a well defined variational principle, one can add an extra 
term in the action. This extra term was introduced by Hawking and Gibbons: 

^/ V^iR+167rGL„,) + -^ f V^K, (A.ll) 

TTLr J »7rG Jqm 
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where A4 is the manifold under consideration, dA4 is the boundary manifold, h 
is the determinant of the induced metric on the boundary and K is the trace of 
the extrinsic curvature of the boundary, Kab = Va?^b, being a vector normal 
to the boundary. 
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A. 2 Symmetries and conserved charges 
A. 2.1 Komar integrals 

In this subsection, we present a technique to compute conserved quantities, 
namely the Komar integrals. We follow the approach of reference Let k 
be a killing vector. Consider the quantity 

abmi...md-2\7, c ,\7Ci,d _ mi...md-2a.b ,V7,V7<^1-'^ 

t V f,fcmi...md_2C<i ^ — *^ '^mi...md-2cd ^ b ^ 

= {d-2y.2\S^^6''JVbV''k'^ (A.12) 
= (d-2)!2!V6V["fc''l 
^ {d^2)\2\Rlk^ 

Contracting equation (jA.12p with eani...nd-i leads to 

V[„je„2...„^_j]c(;V'^^'* = ^_ ^Rbk^iani...nd-ii (A. 13) 



smce 



abmi...md^2c 



(A.14) 



In the vacuum region (typically in asymptotic region), the Ricci tensor 
vanishes. It follows that the left hand side of equation (|A.13|) vanishes as well. 

In other words, considering the Killing form K. = k^dx™, equation (jA.lSP 
states that 

d^dJC^Q, (A.15) 

in regions where Rab — 0. 

Integrating equation (jA.15|) over a volume V bounded by any two [d — 2)- 
spheres S and S' in the exterior vacuum region and applying Stokes theorem 
results in the fact that the integral 



*dlC (A.16) 



s 



is independent on the choice of S and thus is a constant. Expressing the 
form (jA.lSp in components, this integral yields the following expression for a 
conserved charge associated to a Killing vector fc, known as the Komar integral: 

Qk = ^ e™l...™,_.afcVU''d^™-■■™-^ (A.17) 

where (is™i---™d-2 ^]^g surface element on S. 
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A. 2. 2 Background substraction 

The background substraction method has been introduced by Hawking and 
Horowitz in |104| . The main idea is to first consider the Hamiltonian formula- 
tion of general relativity, 

H = — i- / {Nn + N„,nnVld''-^x (A.18) 

where is the sliced manifold with respect to some time coordinate, Pmn is 
the momentum associate to the induced metric on such that Pmn = g^m" 
{Leh being the Einstein- Hilbert Lagrangian and 7„j„ is defined in the next 
equation), N and TV" are respectively the lapse function and the shift vector 
(see [511 1211) associated to the t coordinate, such that the line element takes 
the form 

ds^ = -N^dt^ + {dx"^ - N''dt){dx'' - N''dt)"fab, for some metric jab on Sj, 

(A.19) 

= St n where is the boundary manifold, n™ is the normal unit 
vector to while r™ is the normal unit vector to E°°. (^^"2) jg f^i^^ trace of the 
extrinsic curvature of embedded in E^ and H, W" are constraints associated 
to the diffeomorphism invariance of the theory |26j . 700 is the determinant of 
the induced metric on . 

A solution to the Einstein equations solves the constraints, so the hamilto- 
nian reduces on-shell to 




(A.20) 



Unfortunately, this Hamiltonian evaluated with a solution to the Einstein equa- 
tions is sometimes diverging. The simplest example is the Minkowski spacetime 
in four dimensions expressed in spherical coordinates. The boundary manifold 
is simply the surface of a sphere of radius i? — >■ 00 and the extrinsic curvature 
is cx 1 /R. The lapse function in this case is just 1 and the shift vector is the 
null vector. But the determinant of the induced metric is i?'^sin^ 9, leading to 
a diverging Hamiltonian H ~ J R sin 9d9d(j) for large values of R. 

The idea of Hawking and Horowitz is to consider a background reference 
metric and to subtract the contribution to the Hamiltonian of this reference 
solution. The background has to be static and the solution which we want to 
compute the Hamiltonian should asymptotically be sufficiently close to the 
reference background (see |104j for the precise meaning of sufficiently close). 
Then, since the reference background is static, the conjugate momentum to the 
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induced boundary metric vanishes and the only contribution left to the physical 
Hamiltonian is given by 

Hphys = f ^{n (('^-2)/^ -('^-2) Ko) - iV>™„r") Vl^d^'-^x, 

(A.21) 

where ^'^ ^^Kq is defined in the same way as ^^i^ but for the background ref- 
erence metric. For example, the appropriate reference background for Tangher- 
lini solutions is the c?-dimensional Minkowski space in spherical coordinates. 

A. 3 Linearised gravity 

Let us consider a line element 

ds^ ^ gabdx^dx^, (A.22) 

where gab is a solution of the Einstein equations, i.e. such that Gab = ^T^GTab 
for certain matter fields. 

We can consider a perturbation hab around this solution; the perturbed 
metric then reads 

gab = gab + hab, hab < gab- (A. 23) 

Of course the matter fields are perturbed also, but here, we will focus on the 
gravitational perturbation since the matter perturbation will depend on the 
model under consideration. 

Taking the trace of the Einstein equation and solving for the scalar curvature 
yields the following form of the Einstein equation: 

Rab = SttG (Xab - J^gabT^ , (A.24) 

where T is the trace of the stress tensor and Rab is the Ricci tensor constructed 
out of (1X231) . 

To first order in hab, (|A.24p becomes 

- i {nS^X + '2Ra\') hcd = SnGS (xab - ^~gabT^ ■ (A.25) 

The right hand side is computed assuming that the matter fields are perturbed 
also; 6 then represents the variation of Tab — jz^gabT limited to first order in hab 
and in the matter field perturbation. Equation (|A.25P holds only in the gauge 
such that VaTi"'' = 0, /lab.g"'' = 0; if we don't impose a gauge choice, there are 
terms of the form Va^ch^, Va^bhcdg^"^- Note that the d'Alembertian and the 
Riemann tensor are evaluated using the background metric gab- The derivation 
of equation (|A.25P is straightforward. 
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The left hand side of (|A.25P is the variation of the Ricci tensor and de- 
fines the Lichnerowitz operator (Al)^^'^'* = OS^S^ + 'i.R^^ such that SRab = 



Appendix B 

AdS spacetime 



B.l (A)dS spacetime 



From a geometric point of view, the d-dimensional AdS spacetime is the surface 
of a (d + l)-dimensional hyperboloid embedded in flat spacetime of signature 
(— , —,+,..., +). This hyperboloid and the embedding spacetime are invariant 
under the action of 50(2, d—2). Let us denote the coordinates of the embedding 
spacetime hy U,V, Xi, i ~ 1, . . . ,d — 2; then the equation of the hyperboloid is 
H = 0, the metric and H being given by 

ds^ = -dX^- dXj+6^jdX'dX^, i,j ^l,...,d-l, (B.l) 

d-l 

4=1 

for some i? e M and where Sij is the Kronecker symbol. This is represented in 
figure [BTI in a 3 dimensional euclidean spacetime. 
Defining 

-'^o — Xd Xq + Xd , , 

= t=—, 
Ru ' Ru ' 

the line element on the hyperboloid becomes 

ds^ ^ ^ (~dt^ + Siidx'dx^ + dz^) , i,j = l,...,d-2, (B.3) 

where z = 1/u. These coordinates for AdS spacetime are referred to as the 
Poincare coordinates. 
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Figure B.l: Embedding of the AdS hyperboloid in a 3 dimensional Euclidean 
spacetime. The AdS hyperboloid is embedded in a flat spacetime of signature 
(-,-, + ...,+). 



Note that these coordinates define two patches: z > 0, z < 0, corresponding 
to Xo — Xd > (resp. < 0). These two patches are the Poincare patches. It 
is clear from (jB.3|l that the conformal boundary of AdSd is Aid-i (located at 
z = 0). 

AdS spacetimes has a constant negative curvature; the Riemann tensor is 
given by 

Rabcd~—-^{gacgbd^9ad9bc)- (B.4) 

Note that the Minkowski spacetime obeys the same relation, in the limit 
where R-^ oo. The counterpart of AdS spacetime with a positive curvature is 
the de Sitter spacetime, having also such a symmetric Riemann tensor, except 
that the sign is + instead of — . 

Such spacetimes are said to be maximally symmetric. 

Another choice of parametrisation of the hyperboloid, 



^0 

X,, 



i?^, Xd = 
smp 

RCli tan p, i 



smp 
I d 



(B.5) 



where X^iLi = 1; defining the surface of a unit {d — l)-ball, gives the fol- 



lowing metric form of the AdS line element: 

i?2 



ds^ = 



cos^ p 



-dr'^ + dp^ + sin^ pdill^2) 



(B.6) 
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where (iO^_2 is the square hne element on the (d— 2)-sphere. These coordinates 
cover the whole hyperboloid and are referred to as 'global' AdS coordinates. 
Setting r = i?tanp, the line AdS element becomes 

d-s' = - (l + ^) dt' + ^ + r'dnl^,. (B.7) 

which will be a useful parametrisation for our purpose. 

Using (jB.4l) , the AdS spacetime is obviously a solution to 



_ (rf-l)(rf-2) 

'-'ab — 9ab, [ti.ti) 

Gab being the component of the Einstein tensor, gab the metric components. 
In other words, AdS spacetime can be seen as a universe containing a positive 
energy density Gq and a negative pressure. 

Note that solutions to Einstein equations such that the Einstein tensor is 
proportional to the metric tensor are called Einstein spaces. 



B.2 Asymptotically AdS spacetime 

We will refer to an asymptotically AdS solution as a solution such that the 
asymptotical form of the solution satisfies globally relation (|B.4p . However it 
is not always clear that these asymptotical spacetimes can be recovered as the 
surface of a hyperboloid embedded in a d + 1 dimensional fiat spacetime. 
For example, the d-dimensional AdS black hole: 

ds' = -f{r)dt' + ^^+r'dnl_,, fir)='^ + l-[^y'\ (B.9) 

is clearly asymptotically of the form of (|B.7P ; relation (|B.4[) is immediately 
satisfied and it is obvious that it describes the surface of the AdS hyperboloid 
for large values of r (when rg/r <C 1). 

Another example is the AdS black string, 

ds^ = -b{r)df + -jj-^+r^dnl_^ + a[r)dz^, z € [0,L], (B.IO) 

where a, 6, / are functions to be computed numerically and where the asymp- 
totic form of these function is 

a(r)(x^ + ao + o(^-j , cx ^ + ao + (^-j , /(r) cx - + /0+O (^-j , 

(B.ll) 
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for some values of ao,/o (see e.g. [SO])- This spacetime satisfies (|B.4p asymp- 
totically, but the z coordinate has a finite range. However, it is not clear 
whether this spacetime describes the surface of the AdS hyperboloid or not, 
since the asymptotic form of the black string is the AdS space (|B.7p only to 
the leading order. We will refer to such spacetimes as asymptotically locally 
AdS spacetimes. 



Appendix C 

The Poschl- Teller equation 



The Poschl- Teller equation is known as the following one-dimensional eigen- 
value Schrodinger equation for the corresponding Poschl- Tellel potential. 

NiN + 1) 2 rr^,, 

ax'^ cosh X ^ 

It is considered on an appropriate domain of the Hilbert space of square inte- 
grable functions on the real line. It is a standard result that, for integer values 
of N , the above equation admits -|~ 1 eigenvalues and eigenvectors which can 
be computed algebraically. 

For the first few values of N the eigenvalues are given by 



N ^1:ujI = iV = 2 : = -4, -1, 0, iV = 3 : -9, -4, -1, 0. 



(C.2) 



The corresponding eigenfunctions are given by 



1 



, rjQ = tanhx 



(C.3) 



■n-i = 



cosh a; 



for the case iV = 1, 



?7_4 = tanh^ — 1, ?7-i 



sinhx 



2 — , rjQ = dtanh a; — 1, 



(C.4) 



cosh X 



for the case N — 2 and so on. 
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Appendix D 

Supersymmetric quantum 
mechanics 



In this appendix, we review the basic concepts of supersymmetric quantum 
mechanics, following reference |72| . 

D.l Superpartner potentials and superpotential 

First, consider a Hamiltonian equation for quantum mechanics in 1 dimensions 
of the form ^ ^ 

Assuming that the ground state has zero energy, Eq = 0, it is possible to 
reconstruct the potential (up to a constant) from the zero energy ground state: 

where the prime denotes the derivative with respect to the only coordinate. 
Now, let Hi = A^A, where 

A=^^ + W, A^ = -^^ + W, (D.3) 
V2m ox \j2m ox 

for some function W , called the superpotential. Note that if AifjQ = 0, HiijjQ = 
is automatically satisfied. 

Then, expanding A^ A acting on a test function provides the relation be- 
tween Vi and W: 

Vi = W^- -^W. (DA) 



'2m 
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Note that this is consistent with Aipo = 0; if this is so, we can write W 
as a function of the ground state, leading to W = — 7^ 4^- Inserting this 
expression of W in (|D.4p leads to (|D.2|) . 

Now, let us define another Hamiltonian, H2 ~ AA^ . This Hamiltonian 
corresponds to a quantum theory with a new potential V2 given by 

V2^W^ + -^=W'. (D.5) 

These two potentials Vi and V2 are said to be superpartner potentials. 



D.2 Energy eigenvalues 

The striking feature of superpartner potentials is the link between the energy 
eigenvalues and eigenvectors of each potentials. It turns actually out that the 
operator A and A^ have the effect of taking an eigenvector of Hi to H2 (resp. 
H2 to Hi). 

The eigenvectors of Hi and H2 obey the following equations 

These two relations imply that 

HlA^i^^^ = A^AA^i^l^^ =E^^^AH\^^ 
ilsM^' = AAtA^,W =41)^^(1). (D.7) 



of H2- From the fact that that E^j^^ = and that the energies are positive, we 



which means that Aipif^ is an eigenvector of Hi and A^tp'n^ is an eigenvector 
of H2 ■ From the fact that tha 
obtain the following relations: 
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Numerical method 



In this appendix, we introduce the basic ideas used in the two main algorithms 
used to solve the various differential equations of this thesis, namely Colsys 
and Fidisol. 



E.l Partial differential equations 

As we said in the text, we used the solver Fidisol in order to solve partial 
differential equations. The Fidisol package has been developed by the Karlsruhe 
Institut fiir Technologic and has been checked on a large number of problems. 

The algorithm needs several inputs. First, the effect of introducing a com- 
pact coordinate which maps the semi- infinite range r € [r/i,oo] to x € [0,1], 
such that X = {r — rh)/r leads to the following substitutions in the differential 
equations 

r—^{l-x)— r^—^{l-xf—~2{l~x)— (El) 
dr dx ' dr^ dx^ dx 

for any function Fi. 

The equations for Fi are then discretized on a non-equidistant grid in x 
and z (where z is defined in chapter [5]). Typical grids used have sizes 60 x 40, 
covering the integration region < a: < 1 and < 2: < 1 . 

All numerical calculations are performed by using the program Fidisol, 
which uses a Newton-Raphson method. A detailed presentation of the this 
code is presented in |113| . 

This code requests the system of nonlinear partial differential equations to 
be written in the form P{x, z, u, u^, Uz, Uxx, u^z) = 0, (where u denotes the set 
of unknown functions and , Uz are the derivative of u with respect to x and 
z) subject to a set of boundary conditions on a rectangular domain. The user 
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must deliver to Fidisol the equations, the boundary conditions, the Jacobian 
matrices for the equations and the boundary conditions, and some initial guess 
functions. 

The numerical procedure works as follows: for an approximate solution u*^^-* , 
P(m(^'') does not vanish. The next step is to consider an improved solution 
= 4- supposing that P{u^^'' + Au) = 0. The expansion in the 
small parameter Am gives in the first order = P{u'-^^ + Au) « P(u(i)) + 

This equation can be used to determine the correction Au^^'> = Au. Re- 
peating the calculations iteratively (u^'^^ = u^^^^ + Au^^^^ etc), the approximate 
solutions will converge, provided the initial guess solution is close enough to the 
exact solution. The iteration stops after i steps if the Newton residual 
is smaller than a prescribed tolerance. Therefore it is essential to have a good 
first guess, to start the iteration procedure. 

In each iteration step a correction to the initial guess configuration is com- 
puted. The maximum of the relative defect decreases by a factor of 20 from one 
iteration step to another. However, in some range of the parameters, the con- 
vergence is slower. In this case, we can re-iterated the solution until the defect 
is small enough (about 10^"*). Note, that this defect concerns the discretized 
equations. The estimates of the relative error of the solution (truncation er- 
ror) are computed separately. They can be of the order 0.001. The errors 
also depend on the order of consistency of the method, i.e. on the order of the 
discretization of derivatives. 

E.2 Ordinary differential equations 

The procedure used to integrate systems of ordinary differential equations is 
very similar to the one used for partial differential equations. One starts with 
an initial mesh and an initial profile and provides the equations in the following 
form: 

u\"''^ = F,{x,u,uf^),k = l,...,mj- - 1, (E.2) 

where the superscripts denotes the order of the derivative with respect to the 
variable x. The boundary conditions are provided in the following form: 

G^ia^,UJ,uf)...)=0, fc = l,...,mj -1, (E.3) 

where is the point where the boundary condition is imposed. 

The solver then solves the discretized equations and boundary conditions 
using a Newton-Raphson method; to this end, it requires the derivatives of the 
equations and of the boundary conditions with respect to the functions Uj and 
their derivatives itj*', i = 1, . . . , nij — 1. 
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Here again, it is important to provide an initial guess sufficiently close to 
the real solution. It must be noted that Colsys rearranges the mesh in order 
to minimize the numerical errors. More details can be found in |69j . 

The typical grid we used consisted in 200 points and the relative errors were 
of the order of 10~^ for every cases we considered. 
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of signature (—,—, + ..., +) \lA4 
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